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The  objective  of  this  program  has  been  to  elucidate  the  underlying  mechanisms  by  which 
porous  matrices  give  rise  to  composite  toughness  in  the  absence  of  interphases  between 
the  fibers  and  the  matrix,  the  role  of  microstructure  and  processing  in  attaining  this 
behavior,  and  the  mechanics  associated  with  the  deformation,  cracking  and  failure  of  the 
composites.  The  broader  goal  has  been  to  develop  the  knowledge  base  necessary  to 
exploit  the  nonlinear  phenomena  occurring  within  the  matrix.  Properties  of  prime  interest 
include  in-plane  tensile  strength  and  notch  sensitivity,  subject  to  both  fiber-dominated 
and  matrix-dominated  loadings.  Recognizing  that  changes  in  the  structure  of  the  porous 
matrix  have  opposing  effects  on  the  properties  of  interest  (enhancing  some  while 
compromising  others),  the  research  has  sought  to  identify  the  microstructural  parameters 
that  provide  a  desirable  balance  of  composite  properties.  The  key  variables  explored  were 
those  pertaining  to  the  characteristics  of  the  porosity,  notably  its  content,  as  well  as  the 
nature  and  extent  of  the  bonds  which  determine  both  the  strength  of  the  particle  network 
within  the  matrix  and  the  strength  of  the  fiber-matrix  interface. 

The  highlights  of  the  program  are  summarized  below.  Additional  details  are  in  the 
appended  papers. 


I.  Mullite/Alumina  Mixtures  for  Use  as  Porous  Matrices 

Weakly  bonded  particle  mixtures  of  mullite  and  alumina  have  been  assessed  as  candidate 
matrices  for  use  in  porous  matrix  ceramic  composites.  Conditions  for  the  deflection  of  a 
matrix  crack  at  a  fiber-matrix  have  been  used  to  identify  the  combinations  of  modulus  and 
toughness  of  the  fibers  and  the  matrix  for  which  damage  tolerant  behavior  is  expected  to 
occur  in  the  composite.  Accordingly,  an  experimental  study  on  modulus  and  toughness  of 
particle  mixtures  as  well  as  the  changes  in  these  properties  following  aging  at  elevated 
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temperature  was  designed  and  conducted.  Models  for  these  properties  based  on  bonded 
particle  aggregates  have  been  calibrated.  Upon  combining  the  experimental  and  modeling 
results,  predictions  are  made  of  the  critical  aging  times  at  which  damage  tolerance  is  lost 
because  of  sintering  at  the  particle  junctions  and  the  associated  changes  in  mechanical 
properties.  For  an  aging  temperature  of  1200°C,  the  critical  time  exceeds  10,000  hours 
for  the  mullite-rich  mixtures. 

II.  Controlling  Mechanical  Properties  via  Precursor-Derived  Alumina 

An  assessment  has  been  made  of  the  utility  of  a  precursor-derived  alumina  in  controlling 
both  the  modulus  and  the  toughness  of  mullite-rich  porous  matrices.  Property  changes 
have  been  probed  through  models  of  the  mechanical  behavior  of  bonded  particle 
aggregates.  The  upper  allowable  limit  on  the  concentration  of  precursor-derived  alumina 
to  ensure  damage  tolerance  in  a  fiber  composite  lies  in  the  range  of  7-9%,  depending  on 
the  mullite/alumina  ratio  in  the  particle  slurry  and  the  subsequent  aging  treatment.  Values 
slightly  below  this  limit  yield  composites  that  exhibit  a  desirable  balance  between  fiber- 
and  matrix-dominated  properties. 


III.  Notch  Sensitivity  of  Ceramic  Composites  with  Rising  Fracture  Resistance 

A  theoretical  framework  has  been  developed  for  the  notched  strength  of  ceramic 
composites  that  exhibit  rising  fracture  resistance.  It  is  based  on  established  concepts  of 
crack  stability  under  stress-controlled  loadings.  Using  a  linear  representation  of  the 
resistance  curve  (expressed  in  terms  of  an  energy  release  rate),  analytical  solutions  are 
obtained  for  the  strength  as  well  the  amount  of  stable  crack  growth  preceding  fracture  and 
the  associated  fracture  resistance.  Calculations  have  been  performed  for  several  test 
configurations  commonly  used  for  material  characterization,  including  single-  and  double¬ 
edge  notched  tension,  center-notched  tension,  and  single  edge-notched  bending.  The 
results  reveal  salient  trends  in  strength  with  notch  length  and  specimen  geometry.  An 
assessment  of  the  theory  has  been  made  through  comparison  with  experimental 
measurements  on  an  all-oxide  fiber  composite.  The  utility  of  the  theoretical  results  both 
for  rationalizing  the  trends  in  measured  notched  strength  and  for  guiding  experimental 
studies  of  notch  sensitivity  have  been  demonstrated. 

IV.  Strength  and  Notch  Sensitivity  of  Porous-Matrix  Oxide  Composites 

Effects  of  matrix  strength  on  the  notched  and  unnotched  tensile  properties  of  a  family  of 
porous-matrix  oxide  composites  have  been  examined  both  experimentally  and 
theoretically.  Experiments  were  performed  on  three  composites,  distinguished  from  one 
another  by  the  amount  of  binding  alumina  within  the  matrix.  Increases  in  alumina 
concentration  produce  elevations  in  unnotched  tensile  and  shear  strengths,  but  the 
benefits  are  offset  by  an  increase  in  notch  sensitivity.  The  degree  of  notch  sensitivity  can 
be  rationalized  on  the  basis  of  a  model  that  accounts  for  interactions  between  notch  tip 
tensile  and  shear  bands.  The  model  predictions  are  cast  in  terms  of  the  ratio  of  the  notch 
length  to  a  characteristic  bridging  length  scale.  These  results,  in  turn,  form  the  basis  for  an 
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analytical  formula  for  notched  strength,  accounting  for  effects  of  elastic  anisotropy  and 
finite  sample  size.  The  utility  of  this  formula  in  predicting  notched  strength  has  been 
illustrated.  Issues  associated  with  bridging  law  shapes  and  bridging  length  scales  have 
been  addressed.  The  effect  of  alumina  concentration  on  notch  sensitivity  are  rationalized 
in  terms  of  its  influence  on  the  bridging  length  scale,  dictated  by  the  interplay  between  the 
unnotched  tensile  strength,  the  longitudinal  Young’s  modulus,  the  degree  of  in-plane 
elastic  anisotropy  and  the  fracture  energy.  The  net  result  is  a  decreasing  bridging  length 
scale  and  hence  increasing  notch  sensitivity  as  the  matrix  is  strengthened  with  alumina. 


V.  Effects  of  Combustor  Rig  Exposure  on  a  Porous-Matrix  Oxide  Composite 

Effects  of  exposure  in  a  laboratory  combustor  on  microstructural  stability  and  property 
retention  of  an  all-oxide  fiber-reinforced  ceramic  composite  have  been  explored.  The 
material  consists  of  a  porous  mullite-alumina  matrix  and  Nextel  720  fibers  in  an  8-hamess 
satin  weave.  To  assess  the  effects  of  matrix  strength,  two  matrix  conditions  have  been 
used,  distinguished  from  one  another  by  the  amount  of  alumina  added  through  precursor 
impregnation  and  pyrolysis  (1.8%  and  4.8%).  In  both  cases,  the  dominant  damage  mode 
upon  exposure  involves  interply  delamination  along  the  panel  mid-plane.  However, 
significant  reductions  in  the  rate  and  extent  of  cracking  are  obtained  in  the  material  with 
higher  alumina  content:  a  result  of  the  higher  delamination  resistance.  Mechanical  tests 
performed  on  exposed  specimens  reveal  a  slight  (10-20%)  reduction  in  tensile  strength 
along  the  fiber  direction  and  a  comparable  increase  in  shear  strength.  These  trends  suggest 
some  sintering  of  the  matrix  upon  exposure.  Examinations  of  fracture  surfaces  provide 
additional  supporting  evidence.  Implications  for  long-term  performance  and  strategies  for 
imparting  improvement  in  microstructural  stability  and  delamination  resistance  have  been 
elucidated. 
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A  METHOD  FOR  DETERMINING  THE  STRESS  INTENSITY  FACTOR  OF  A 
SINGLE  EDGE-NOTCHED  TENSILE  SPECIMEN 

M.A.  Mattoni  and  F.  W.  Zok 
Materials  Department,  University  of  California 
Santa  Barbara,  CA  93106 
e-mail:  zok  @  engineering,  ucsb.  edu 

Abstract  A  straightforward  method  for  determining  the  stress  intensity  factor  of 
a  single  edge-notched  tensile  specimen  with  unknown  end  conditions  is  presented.  It 
involves  the  use  of  measured  crack  mouth  opening  displacements,  coupled  with 
existing  solutions  for  stress  intensity  factors  and  crack  opening  displacements  for 
uniform  tensile  loading  and  for  pure  bending.  The  use  of  the  method  is  demonstrated 
through  experiments  on  a  porous  matrix  ceramic  composite. 

1.  Introduction.  The  single  edge-notched  tensile  (SENT)  specimen  is  commonly 
used  for  measuring  the  fracture  resistance  of  engineering  materials.  The  pertinent 
existing  solutions  for  stress  intensity  factors  are  based  on  two  limiting  idealizations  of 
the  conditions  at  the  specimen  ends  where  the  force  is  transmitted  from  the  testing 
machine  to  the  specimen.  In  the  first,  the  loading  fixture  is  assumed  to  have  no 
rotational  stiffness  and  the  stress  is  distributed  uniformly  across  the  specimen  ends. 
To  achieve  this  condition,  the  specimen  length,  H,  must  be  significantly  greater  than 
the  width,  W,  and  the  fixture  must  have  a  mechanism  for  allowing  rotation  with 
minimal  constraint.  The  latter  condition  is  usually  achieved  by  pin  loading.  The  stress 
intensity  factor  for  this  case  is  given  by  standard  handbook  solutions,  e.g.,  Tada  et  al. 
(1985).  In  the  second,  the  loading  fixture  is  assumed  to  be  rigid  and  the  axial 
displacements  are  uniform  across  the  specimen  end.  In  this  case,  a  bending  moment 
is  induced  at  the  specimen  ends,  in  addition  to  the  applied  tension.  To  achieve  these 
conditions,  a  significant  length  of  the  specimen  at  the  ends  must  be  clamped  using 
stiff  grips.  Pertinent  solutions  for  the  stress  intensity  factor  for  the  uniform 
displacement  condition  have  been  developed  by  Marchand  el  al.  (1986),  Blatt  et  al. 
(1994),  John  et  al.  (1995),  and  John  and  Rigling  (1998). 

The  differences  in  the  stress  intensity  factors  in  the  two  cases  are  usually 
significant.  For  example,  for  a  test  specimen  with  aspect  ratio  H/W=2  and  a  crack 
length  a/W=0.5,  the  normalized  stress  intensity  factors  for  the  uniform  stress  and  the 
uniform  displacement  conditions  are  K/cWra  =  2.83  and  1.29,  respectively:  a 
difference  of  more  than  a  factor  of  2.  It  is  expected  that  stress  intensity  factors  for  test 
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conditions  that  yield  intermediate  levels  of  constraint  will  not  be  represented 
adequately  by  either  of  the  two  limiting  cases. 

Here  we  present  a  straightforward  method  for  determining  stress  intensity  factors 
for  the  SENT  geometry  under  conditions  in  which  neither  of  the  two  limiting 
conditions  are  obtained.  The  development  of  the  method  has  been  motivated  by  the 
need  to  assess  fracture  resistance  of  ceramic  composites  that  are  currently  being 
fabricated  in  small  quantities  and  for  which  the  test  specimens  are  necessarily  small. 
Often  this  precludes  the  use  of  long  tensile  specimens  for  the  purpose  of  gripping 
long  sections  near  the  ends  and  achieving  uniform  displacement  conditions. 
Furthermore,  compliant  end  tabs  are  usually  used  to  ensure  a  uniform  pressure 
distribution  along  the  gripped  surface  and  to  prevent  crushing.  Such  tabs  reduce  the 
grip  stiffness  and  compromise  the  goal  of  a  uniform  displacement  condition.  Pin 
loading  configurations  are  also  impractical  because  of  the  potential  for  bearing 
failure.  This  is  due  to  the  high  stresses  needed  to  initiate  crack  tip  failure  coupled 
with  the  low  crushing  strength  of  the  matrix  in  the  composites  of  interest.  As  a  result, 
the  uniform  stress  condition  is  also  difficult  to  achieve. 


Figure  1.  A  schematic  showing  the  uniform  tensile  stress  and  the  bending  moment  at  the  ends  of  an 
SENT  specimen.  Note  that  the  stress  intensity  factor  from  the  moment  is  negative. 

2.  Analysis.  When  neither  of  the  limiting  end  conditions  is  achieved,  the  stresses 
at  the  specimen  ends  can  be  partitioned  into  two  components:  (i)  a  uniform  tensile 
stress.  Or,  equal  to  the  average  (measured)  stress,  and  (ii)  a  bending  stress,  Ob, 
proportional  to  Or,  but  otherwise  of  unknown  magnitude  (Fig.  1).  The  stress  intensity 
factor  associated  with  these  two  stress  components  is: 

K  =  Fj-(a  /  W)  aTViia  -  FB(a  /  W)  a  Bi/jta  =  FE(a  /  W)  Oj-Vrca  (1) 

where  Fr(a/W)  and  FB(a/W)  are  the  usual  non-dimensional  functions  for  tension  and 
bending  given  in  handbooks  (e.g.  Tada  el  al.,  1986),  and  Fn(a/W)  is  the  function  that 
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characterizes  the  normalized  stress  intensity  factor  in  the  present  test  configuration. 
The  corresponding  crack  mouth  opening  displacement,  8,  is: 


5  =  VT(a/W)^|±-VB(a/W)i^ 


VE(a  /  W) 


4  O'j’  a 
E 


(2) 


where  E  is  the  plane  strain  Young’s  modulus,  the  functions  VT(a/W)  and  VB(a/W)  are 
obtained  from  handbooks,  and  V£(a/W)  characterizes  the  normalized  crack  opening 
in  the  test  configuration.  Rearranging  Eqn.  2  yields  an  expression  for  the  bending 
stress: 


£b.  _  Vt  i  Tse" 

VB  VB  vaCf-* 


(3) 


Combining  Eqn.  (1)  and  (3)  gives  the  stress  intensity  factor  as: 


This  result  indicates  that  the  stress  intensity  factor  can  be  obtained  from  a 
combination  of  the  known  functions  and  the  normalized  crack  mouth  opening 
displacement,  VE  =  8  E/aoT ,  the  latter  being  obtained  experimentally. 


3.  Experiments.  The  analysis  was  assessed  through  a  series  of  SENT  tests  on  a 
porous  matrix  ceramic  composite.  The  composite  consisted  of  Nextel  720  fiber  cloth 
in  an  8  harness  satin  weave  and  a  porous  mullite-alumina  matrix.  The  matrix  was 
introduced  in  two  steps:  (i)  pressure-assisted  infiltration  of  a  mullite-alumina  slurry 
through  the  fiber  preform,  and  (ii)  impregnation  and  pyrolysis  of  an  AI2CI  (OH)s 
solution  (an  alumina  precursor).  Processing  details  were  essentially  identical  to  those 
described  by  Mattoni  el  al.  (2001),  with  two  exceptions.  First,  the  concentration  of 
the  precursor  solution  was  increased  to  give  a  volumetric  yield  of  ~  8.3%  (compared 
with  3%  in  previous  studies).  The  more  concentrated  solution  allows  a  prescribed 
matrix  density  to  be  achieved  with  fewer  impregnation  and  pyrolysis  cycles.  In  the 
present  study,  4  such  cycles  were  used,  yielding  a  final  matrix  porosity  of  31.9%. 
Second,  to  prevent  segregation  of  the  precursor  during  drying,  the  panels  were  gelled 
immediately  after  precursor  impregnation  by  immersing  them  in  a  NH4OH 
environment  for  6  hours.  The  absence  of  segregation  was  confirmed  through  a  series 
of  matrix  hardness  measurements  at  various  points  through  the  plate  thickness. 


L6 


The  aspect  ratio  of  the  gauge  section  of  the  SENT  specimens  was  fixed  at 
H/W=2.1.  The  width  of  all  specimens  was  W=24  mm  whereas  the  notch  length  was 
varied  over  the  range  a/W=0.35-0.65.  To  facilitate  gripping,  thin  fiberglass  tabs 
were  adhesively  attached  to  the  ends  over  a  length  of  approximately  4  mm.  The 
specimens  were  instrumented  with  a  displacement  gauge  at  the  crack  mouth  and  a 
crack  gauge  at  the  notch  tip.  Corroborating  measurements  of  crack  length  were 
obtained  through  direct  optical  observations  of  the  notch  tip  region  using  a 
stereomicroscope. 


Figure  2.  Effects  of  end  condition  on  (a)  the  CMOD  response  and  (b)  the  bending  stress. 


Table  I  Functions  for  Single  Edge-Notched  Specimen  under  Uniform  Tension  and  Bending 
(from  Tada  et  al.  1985). 
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The  variation  in  CMOD  with  applied  stress  for  one  specimen  configuration  is 
plotted  in  Fig.  2(a).  Combining  measurements  of  this  type  with  Eqn.  (3)  and  the 
functions  in  Table  I  yields  the  normalized  bending  stresses  plotted  in  Fig.  2(b).  Also 
shown  for  comparison  in  these  figures  are  the  predicted  values  for  uniform  stress  and 
uniform  displacement  conditions,  using  solutions  in  Tada  et  al.  (1985)  and  John  and 
Rigling  (1985).  The  experimental  results  lie  between  the  predictions  for  the  two 
limiting  cases,  as  expected.  Evidently  the  bending  stress  is  small  for  a/W=0.35  but 
subsequently  increases  rapidly  with  increasing  a/W. 

The  preceding  results  were  combined  with  the  load  and  crack  length 
measurements  and  the  result  in  Eqn.  (4)  to  obtain  the  initiation  fracture  toughness,  Ko, 
as  well  as  the  fracture  resistance,  Kr,  subsequent  to  the  onset  of  cracking.  The  results 
are  plotted  in  Fig.  3.  The  initiation  toughness  is  independent  of  the  notch  length,  at  a 
value  Ko=7.8±0.7  MPaVm.  The  fracture  resistance  curves  increase  linearly  with 
crack  extension,  at  a  rate  independent  of  the  initial  notch  length.  The  consistencies  in 
Ko  and  Kr  between  specimens  of  varying  notch  length  provide  confidence  that  the 
present  method  of  determining  the  pertinent  stress  intensity  factors  is  indeed  correct. 


Figure  3.  Effects  of  notch  length  on  (a)  the  initiation  fracture  toughness  and  (b)  the  fracture  resistance 
curve. 

4.  Concluding  Remarks.  The  present  method  for  determining  the  stress  intensity 
factor  provides  flexibility  in  the  test  configurations  used  for  determining  fracture 
resistance.  Specifically,  it  provides  latitude  in  the  degree  of  constraint  exerted  by  the 
grips  without  compromising  the  integrity  of  the  fracture  resistance  measurements.  In 
addition  to  its  use  in  test  configurations  that  are  known  to  experience  intermediate 
end  conditions,  the  method  would  be  useful  for  checking  whether  a  particular 
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targeted  end  condition,  either  uniform  stress  or  uniform  displacement,  has  indeed 
been  achieved. 
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Weakly  bonded  particle  mixtures  of  mullite  and  alumina  are 
assessed  as  candidate  matrixes  for  use  in  porous  matrix 
ceramic  composites.  Conditions  for  the  deflection  of  a  matrix 
crack  at  a  fiber-matrix  interface  are  used  to  identify  the 
combinations  of  modulus  and  toughness  of  the  fibers  and  the 
matrix  for  which  damage-tolerant  behavior  is  expected  to 
occur  in  the  composite.  Accordingly,  the  present  study  focuses 
on  the  modulus  and  toughness  of  the  particle  mixtures,  as  well 
as  the  changes  in  these  properties  following  aging  at  elevated 
temperature  comparable  to  the  targeted  upper-use  tempera¬ 
ture  for  oxide  composites.  Models  based  on  bonded  particle 
aggregates  are  presented,  assessed,  and  calibrated.  The  exper¬ 
imental  and  modeling  results  are  combined  to  predict  the 
critical  aging  times  at  which  damage  tolerance  is  lost  because 
of  sintering  at  the  particle  junctions  and  the  associated  changes 
in  mechanical  properties.  For  an  aging  temperature  of  1200°C, 
the  critical  time  exceeds  10  000  h  for  the  mullite-rich  mixtures. 


1.  Introduction 

The  use  of  porous  matrixes  to  enable  damage  tolerance  in 
continuous-fiber  ceramic  composites  (CFCCs)  has  emerged  as 
a  new  paradigm  in  high-performance  materials.1  The  concept 
obviates  the  need  for  fiber  coatings  to  affect  crack  deflection, 
thereby  providing  opportunities  for  lower-cost  manufacturing 
relative  to  that  of  conventional  coated-fiber  systems.  Furthermore, 
on  selection  of  all-oxide  constituents,  the  long-term  durability 
requirements  in  the  targeted  high-temperature  applications  can  be 
met. 

Although  the  porous  matrix  concept  offers  new  opportunities 
for  the  development  of  damage-tolerant  CFCCs,  it  also  presents 
challenges  in  the  design  and  synthesis  of  microstructures  that  meet 
the  opposing  property  requirements  placed  on  the  matrix.  That  is, 
in  the  absence  of  fiber  coatings,  the  matrix  must  be  sufficiently 
weak  to  enable  damage  tolerance  under  fiber-dominated  loadings, 
yet  retain  adequate  strength  to  ensure  acceptable  off-axis  proper¬ 
ties.  In  principle,  the  combination  of  properties  can  be  tailored 
through  changes  in  the  state  of  the  matrix.  For  instance,  improve¬ 
ments  in  the  interlaminar  strength  and  off-axis  in-plane  strength 
can  be  obtained  by  reducing  the  matrix  porosity,  however,  these 
improvements  come  at  the  expense  of  a  reduction  in  the  damage 
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tolerance  under  fiber-dominated  loadings.  These  offsetting  effects 
suggest  the  existence  of  an  optimum  in  matrix  properties  at  which 
a  prescribed  balance  of  properties  is  attained.  However,  ihc 
relationships  between  matrix  structure  and  composite  performance 
are  understood  presently  at  only  a  rudimentary  level.  Conse¬ 
quently,  the  pathway  to  optimization  remains  ill  defined. 

The  objective  of  the  present  article  is  to  elucidate  some  of  the 
connections  between  composition,  microstructure,  and  mechanical 
properties  of  a  candidate  family  of  porous  matrices,  with  emphasis 
on  the  conditions  needed  to  enable  damage  tolerance  in  liber- 
dominated  loadings.  Effects  of  the  matrix  on  off-axis  properties  are 
addressed  elsewhere.2  The  work  steins  from  concurrent  activities 
on  the  development  of  all-oxide  CFCCs  for  use  in  future  gas- 
turbine  systems.1-6  Because  of  the  interest  in  mullite/alumina 
fibers  for  use  in  these  applications,  the  present  article  focuses 
specifically  on  mullite/alumina  matrixes.  These  matrices  are  both 
chemically  compatible  with  the  fibers  and  exhibit  mechanical 
characteristics  that  make  them  attractive  for  use  in  CFCCs.  A 
perspective  based  on  the  mechanics  of  crack  deflection  is  used  to 
guide  the  experimental  measurements  and  provide  a  framework  for 
interpreting  the  property  values.  Consideration  is  given  to  the 
stability  of  the  properties  following  long-term  exposure  at  elevated 
temperature,  comparable  to  the  targeted  upper  use  temperature  for 
oxide  CFCCs. 

We  present  a  review  of  the  mechanics  of  crack  deflection  tit  a 
fiber-matrix  interface.  The  results  are  used  to  identify  the  critical 
combination  of  matrix  modulus  and  matrix  toughness  needed  for 
deflection  and  to  motivate  the  subsequent  mechanical  measure¬ 
ments.  The  nature  of  the  candidate  matrix  system,  the  processing 
route,  and  the  measurement  procedures  arc  described.  Results  for 
modulus  and  toughness  are  presented,  along  with  predictions  of 
these  properties  from  models  of  bonded-particle  aggregates.  Fi¬ 
nally,  we  address  the  implications  on  crack  deflection  in  CFCCs. 
For  this  purpose,  calibrated  models  arc  used  to  establish  the  critical 
aging  time  at  which  the  crack  deflection  condition  is  no  longer 
satisfied  and  thus  the  damage  tolerance  is  lost. 


II.  Mechanics  of  Crack  Deflection 

As  a  minimum  requirement  for  damage  tolerance  in  CFCCs. 
cracks  in  the  matrix  must  either  arrest  at  or  deflect  into  the 
fiber/matrix  interface  rather  than  penetrate  into  the  fibers.  The 
conditions  that  satisfy  this  requirement  are  obtained  from  the  He 
and  Hutchinson7  diagram  shown  in  Fig.  I.  Crack  arrest  or 
deflection  is  predicted  to  occur  when  the  ratio  of  interlace 
toughness.  13,  to  fiber  toughness.  I  ,,  falls  below  the  corresponding 
ratio  of  energy  release  rates.  CJG  .  associated  with  deflection  into 
the  boundary  and  penetration  into  the  liber.  The  critical  ratio  is 
controlled  by  the  elastic  mismatch  parameter: 
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Fig.  1.  Conditions  for  crack  deflection  at  a  fiber-matrix  interface  (solid 
line).  Experimental  measurements  from  this  study  are  superimposed 
(symbols).  The  latter  assume  a  toughness  ratio  u>  m  r/r„,  =  1  and  fiber 
properties  Tr  =  15  S/m2  and  Ef  =  260  GPa. 


where  E  is  the  plane  strain  modulus,  and  the  subscripts  /  and  m 
refer  to  fiber  and  matrix,  respectively.  For  porous  matrix  systems, 
a  takes  on  relatively  high  values  (>0.5),  and  hence  the  allowable 
toughness  ratio  is  also  relatively  high.  In  the  systems  of  present 
interest,  the  matrix  and  fiber  constituents  are  essentially  the  same, 
i.e.,  mullite/alumina,  and  hence  the  nature  of  bonding  at  the 
fiber-matrix  interface  is  similar  to  that  between  particles  in  the 
matrix.  Consequently,  the  interface  toughness  is  expected  to  be 
proportional  to  the  matrix  toughness:  17  =  col^  where  to  is  a 
nondimensional  parameter  of  order  unity.  Because  of  the  reduced 
efficiency  of  particle  packing  along  the  fiber  surfaces,8  it  is 
expected  that  17  <  Tm  such  that  u>  <  I.  For  conservative  design, 
to  is  taken  to  be  1 . 

The  critical  combination  of  matrix  and  fiber  properties  are 
obtained  in  the  following  way:  For  a  2  0,  the  energy  release  rate 
ratio,  Gj/Gp,  can  be  represented  by  an  approximate  empirical 
equation: 

1 

7TP  =  4(1  -  a)09  (2) 

This  equation  is  an  exact  representation  of  the  He  and  Hutchinson7 
criterion  at  a  =  0  and  has  an  error  <4%  over  the  range  0  £  a  s 
0.95.  On  setting  GJGp  =  I7/r,  in  Eq.  (2)  and  combining  the  result 
with  Eq.  (1),  the  deflection  condition  is  re-expressed  as  follows 
(Here  the  Poisson’s  ratios  of  the  fiber  and  the  matrix  are  assumed 
to  be  the  same.  Consequently,  the  plane  strain  modulus  ratio  ErIEm 
can  be  replaced  with  Young’s  modulus  ratio  £,/Em). 

s-°-,340(|+t.r>“  <3» 

Here  2  is  a  nondimensional  parameter  that  characterizes  the 
propensity  for  crack  deflection.  It  identifies  the  critical  combina¬ 
tion  of  r7rr  and  £m/£,  needed  to  cause  deflection. 

In  light  of  this  prediction,  the  subsequent  experimental  study 
focuses  on  the  matrix  properties  Tm  and  £m,  with  a  view  to 
assessing  whether  the  crack-deflection  condition  in  Eq.  (3)  is 
satisfied  when  the  matrix  is  combined  with  the  specified  fiber, 
notably  Nextel  720™  (3M  Corp.,  St.  Paul.  MN). 


III.  Materials  and  Test  Procedures 

Two  generic  matrix  topologies  have  previously  been  explored.6 
The  first  involves  the  use  of  particle  mixtures  of  muiiite  and 
alumina  and  relies  on  the  more  rapid  sintering  kinetics  of  alumina 
to  bond  the  network  together.  The  second  is  based  on  the  use  of 


alumina  precursor  liquids  in  conjunction  with  either  pure  mullitc 
particles  or  mullite/alumina  particle  mixtures.  The  latter  approach 
has  the  advantage  that  it  does  not  compromise  the  contiguity  of  the 
particle  network  and  should  yield  a  more  stable  pore  structure.  But 
it  comes  at  the  expense  of  additional  steps  in  material  processing, 
beyond  that  associated  with  particle  mixtures  alone.  The  present 
paper  explores  the  former  topology  (based  on  particle  mixtures): 
the  latter  is  developed  elsewhere.6 

Porous  mullite-alumina  compacts  were  fabricated  using  a  pres¬ 
sure  infiltration  technique.  The  compositions  ranged  from  pure 
mullite  to  60%  mullite  +  40%  alumina  by  volume.  The  compo¬ 
sitions  are  henceforth  designated  by  the  percentages  of  mullitc  and 
alumina,  i.e.,  80M/20A  denotes  80%  mullite  +  20%  alumina.  A 
small  number  of  pure  alumina  specimens  were  fabricated  also,  for 
the  purpose  of  assessing  the  efficacy  of  the  mullite  network  in 
inhibiting  shrinkage  associated  with  sintering  of  alumina  in  the 
mixed  compacts.  The  mullite  particulates  were  mullite  (MU- 1 07. 
Showa  Denko  K.K.,  Tokyo.  Japan).  The  size  distribution  was 
rather  broad,  characterized  by  a  mean  diameter  d  =  —  I  pan  and 
diameters  at  the  1 0th  and  90th  percentiles  of  the  distribution  of  d,„ 
=  0.4  p.m  and  d()0  =  2.3  pan,  respectively.  The  alumina  particu¬ 
lates  (AKP-50.  Sumitomo  Chemical.  Tokyo.  Japan),  had  size 
distribution  characterized  by  d  =  —0.2  p.m,  d,()  =  0.1  pun,  and 
dqo  =  0.3  p.m.  An  aqueous  slurry  containing  the  powder  mixture 
was  prepared  and  vacuum-infiltrated  into  a  mold  with  dimensions 
10  mm  X  50  mm  X  90  mm.  The  green  compacts  were  dried  and 
sintered  at  900°C  for  2  h  to  promote  partial  sintering  at  the  particle 
junctions.  The  compacts  were  machined  into  mechanical  test 
specimens  and  the  specimens  given  a  2  It  heat  treatment  at  I200"C. 
To  assess  the  stability  of  the  microstructure  under  conditions  that 
mimic  the  targeted  service  conditions  for  oxide  CFCCs.  some  test 
specimens  were  given  additional  aging  treatments  of  100  h  or 
1000  h  at  I200°C. 

The  porosity.  P.  of  the  compacts  was  determined  from  mea¬ 
surements  of  mass,  both  in  air  and  immersed  in  water.10  The 
pertinent  densities  of  the  constituent  powders  needed  for  calculat¬ 
ing  porosity  were  measured  using  a  pycnometer.  The  measure¬ 
ments  on  the  powders  were  performed  after  the  powders  had  been 
subjected  to  the  same  heat  treatments  as  those  of  the  compacts. 
Even  after  the  longest  aging  treatments,  the  powders  had  not 
sintered  together  and  remained  free  Bowing.  The  measurements 
revealed  a  very  slight  increase  in  the  density  of  the  pure  mullitc 
powder,  from  3.14  to  3.24  g/cm\  as  the  aging  time  increased  from 
2-1000  h.  The  associated  mechanism  is  unknown,  although  the 
values  are  in  broad  agreement  with  those  reported  in  the  litera¬ 
ture."  In  contrast,  the  alumina  powder  density  remained  un¬ 
changed  over  the  same  range  of  aging  time,  at  4.01  ±  0.01  g/cin'. 

The  Young’s  modulus  of  the  fired  compacts  was  measured 
using  strain-gauged  flexure  specimens.  5  mm  x  5  mm  x  50  mm 
long,  loaded  in  four-point  bending.  The  inner  and  outer  loading 
spans  were  20  and  40  mm,  respectively.  The  toughness  was 
obtained  from  work  of  fracture  measurements  on  chevron-notched 
bend  specimens1’  with  the  same  dimensions,  but  loaded  in 
three-point  bending  and  instrumented  with  an  LVDT  to  measure 
load-point  displacement.  In  all  cases,  the  crack  grew  stably  across 
the  specimen  cross  section.  The  fracture  surfaces  were  examined 
subsequently  via  scanning  electron  microscopy  (SEM). 

Henceforth,  the  '  symbol  is  used  to  denote  properties  of  the 
porous  compacts,  i.e.,  E  and  F  for  modulus  and  toughness, 
respectively;  the  subscripts  A  and  M  arc  used  to  denote  pure 
alumina  and  pure  mullite:  and  the  absence  of  a  subscript  implies  a 
particle  mixture. 


IV.  Pure  Mullite 
(1)  Experimental  Results 

The  variation  in  the  porosity  of  the  pure  mullitc  compacts  with 
aging  time  is  plotted  in  Fig.  2.  Despite  the  slight  change  in  powder 
density,  the  porosity  of  the  compacts  remains  essentially  constant 
over  the  entire  range  of  aging  times,  within  0.2%  of  its  initial 
value.  £„. 
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Fig.  2.  (a)  Effects  of  composition  and  aging  time  on  compact  porosity. 
The  data  for  0  h  and  2  h  are  virtually  identical  to  one  another  and  are 
indistinguishable  on  the  graph,  (li)  Results  in  (a),  normalized  by  the  initial 
(green)  porosity  at  the  same  composition.  The  solid  line  for  alumina  in  the 
inset  was  obtained  by  dilatometry.  the  symbols  are  based  on  the 
Archimedes  measurements,  and  all  other  lines  are  simply  curve  fits  through 
the  discrete  data  points. 


Both  the  modulus  and  the  toughness  of  the  mullite  compacts 
increase  appreciably  with  aging  time,  by  a  factor  of  —3-4  (Fig.  3). 
These  property  changes  indicate  sintering  at  the  particle  junctions, 
but  the  absence  of  porosity  change  implies  a  mechanism  that  does 
not  involve  movement  of  the  particle  centers.  Possible  mechanisms 
include  transport  of  matter  from  surface  sources  via  either  surface 
diffusion,  lattice  diffusion,  or  vapor  transport.1''  The  dominant 
transport  mechanism  is  inferred  from  sintering  models  and  the 
mechanical  measurements  in  the  manner  described  in  the  follow¬ 
ing  section. 

SEM  observations  of  the  fracture  surfaces  reveal  only  a  small 
number  of  well-defined  broken  junctions  (Figs.  4(a)  and  (b)). 
These  are  somewhat  more  prevalent  after  the  longest  aging  lime 
because  of  the  increase  injunction  size.  The  small  number  of  such 
junctions  indicates  that  fracture  normally  occurs  “cleanly”  without 
appreciable  crack  meandering,  suggesting  a  low  fracture  energy.  It 
also  precludes  direct  measurement  of  junction  size. 

(2)  Modeling 

When  the  junction  radius,  a.  is  much  smaller  than  the  particle 
radius,  R,  junction  growth  owing  to  sintering  follows  a  power  law 
of  the  form:1-'14 


Fig.  3.  Effects  of  aging  time  on  Young's  modulus  and  toughness  of  pure 
mullite.  Power  law  fits  of  the  experimental  data  yield  exponents  of  0.18  ; 
0.03  and  0.28  ±  0.03  for  modulus  and  toughness,  respectively,  in  accord 
with  those  predicted  for  surface  diffusion  controlled  sintering,  notably 
0.143  and  0.286. 


where  /  is  sintering  time,  /R  is  a  reference  time  and  n  is  a  constant; 
both  tR  and  it  depend  on  the  transport  mechanism.  Sintering 
models  yield  values  of  n  =  3  for  vapor  transport,  n  =  5  for  lattice 
diffusion,  and  ;t  =  7  for  surface  diffusion. 

In  turn,  for  monosized  spherical  particles  with  small  junctions 
(a/R  <<  I ),  the  Young’s  modulus  of  a  bonded  particle  aggregate 
scales  linearly  with  junction  radius  in  accordance  with:1516 


where  Ep  and  v  are  the  Young's  modulus  and  Poisson’s  ratio  of  the 
particles;  z  is  the  particle  coordination  number  (-6  for  random 
packing);  D  is  the  relative  packing  density;  £  is  a  numerical 
parameter;  and  g(v)  is  given  by 

5  -  4v 

giv)  =  (I  —  v5)(5  —  3v)  (6) 

For  values  of  v  in  the  range  0£rs  0.25.  g(v)  =  1  i  0.004;  the 
dependence  on  v  is  extremely  small  and  is  subsequently  neglected. 
This  result,  with  £  =  I,  was  derived  by  Walton."  using  Hertzian 
contact  mechanics  to  describe  junction  stiffness  and  assuming  a 
uniform  aggregate  strain.  A  more  rigorous  model  that  accounts  for 
finite  junction  size  and  multiparticle  interactions  yields  results  in 
the  domain  of  small  junctions  (a/R  £  0.3)  that  are  consistent  with 
Eq.  (5)  when  £  is  taken  to  be  —0.76  (simulation  of  modulus  has 
been  performed  for  two  limiting  cases;  assuming  that  the  torsional 
stiffness  of  the  junctions  is  either  zero  or  infinite.16  The  two  sets 
of  results  are  fit  well  by  Eq.  (5).  using  £  =  0.65  and  0.88. 
respectively.  The  average  of  these,  £  =  0.76.  is  used  in  the  present 
work).  Combining  Eqs.  (4)  and  (5)  and  taking  g(  v)  =  I  yields  the 
time-dependence  of  the  modulus: 


The  relationship  between  toughness  and  junction  radius  for  a 
bonded  particle  aggregate  has  been  obtained  from  numerical 
simulations  of  fracture,  using  a  technique  based  on  ihc  discrete 
element  method  (DEM).  Details  of  the  method  and  ils  implemen¬ 
tation  are  presented  elsewhere.17  Briefly,  a  computer-simulated 
sintering  algorithm  is  used  to  generate  a  random  three-dimensional 
array  of  spherical  particles  with  a  prescribed  junction  radius.  The 
simulated  aggregates  consist  of  —1000  particles.  A  crack  is 
defined  by  a  plane  separating  particles  that  have  had  the  junctions 
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Fig.  4.  Fracture  surfaces  of  compacts  that  cover  the  range  of  alumina  content,  from  0%  (a,  b)  to  40%  (c.  d),  following  a  1000  h  aging  at  I200"C.  Examples 
of  fractured  particle  junctions  are  highlighted  by  the  arrows. 


between  them  broken.  The  simulation  proceeds  by  incrementally 
increasing  the  remote  displacement  while  allowing  the  junctions  at 
the  crack  tip  to  fail  when  the  junction  stress  reaches  a  critical 
value,  ac. 


0.05  0.1  0.2  0.3  0.4  0.5 

Relative  junction  size,  a/R 


Fig.  5.  Results  of  fracture  simulations  based  on  the  discrete  element 
method. 


where  17  is  the  intrinsic  junction  toughness.  The  results  of  the 
simulations  (Fig.  5)  are  well  described  by 

f  /a\2 

Fp  l2\^)  (9) 

Combining  Eqs.  (4)  and  (9)  yields  the  time  dependence  of  the 
toughness: 


The  predictions  of  modulus  and  toughness  are  compared  with 
the  experimental  measurements  in  Fig.  4.  The  best  correlations  are 
obtained  by  taking  n  =  7  in  Eqs.  (7)  and  (10).  The  inference  is  that 
junction  growth  occurs  predominantly  by  surface  diffusion.  For 
this  mechanism,  the  reference  time  is  given  by13 

ATR1 

,R  =  m&Ayft  ( ' 1  * 

where  k  is  Bollzman’s  constant.  T  is  absolute  temperature.  5S  is  the 
effective  surface  thickness,  Ds  is  the  surface  diffusion  coefficient, 
y  is  the  surface  energy  (®“1  J/m2).  and  ft  is  the  molecular  volume 
(2.2  X  I0-28  m3  for  mullite). 

A  further  assessment  of  the  models  was  made  in  the  following 
way:  The  numerical  values  of  the  various  parameters  for  mullite 
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are  known,  with  the  exception  of  the  junction  toughness  17  and  the 
effective  surface  diffusivity  5 SDS.  The  diffusivity  was  inferred 
from  the  modulus  data  by  equating  the  proportionality  constant 
obtained  from  a  1/7  power  law  fit  of  the  experimental  measure¬ 
ments  (Fig.  3)  with  the  value  predicted  by  Eq.  (7).  This  yields 
hsDs  =  3.8  X  10_3°  m3/s.  The  junction  toughness  was  then 
obtained  by  equating  the  proportionality  constant  obtained  from  a 
2/7  power  law  fit  of  the  toughness  measurements  with  the  value 
predicted  by  Eq.  (9),  making  use  of  the  diffusivity  inferred  from 
the  modulus  measurements.  The  result  is  17  ~  3  J/m2.  This  value 
is  only  slightly  higher  than  the  surface  energy  contribution,  T  = 
2y  »»  2  J/m2,  and  well  below  the  toughness  of  bulk  polycrystalline 
mullite,  F  «=  20-30  J/m2."  The  latter  toughness  has  significant 
contributions  from  crack  tilting  and  twisting  and  associated  step 
formation  when  cracks  propagate  by  transgranular  cleavage.18 
This  mechanism  is  not  operative  at  the  particle  junctions  because 
the  particles  themselves  are  single  crystals,  and  the  junctions 
between  them  are  essentially  grain  boundaries. 

V.  Mullite/AIumina  Mixtures 

(1)  Experimental  Results 

The  physical  properties  of  the  mullite/alumina  mixtures  are 
summarized  in  Fig.  2.  The  initial  porosity  exhibits  a  shallow 
minimum  at  XA  *«  0.2-0.3,  a  consequence  of  the  slight  size 
difference  of  the  two  particle  types  (Fig.  2(a)).  However,  the 
minimum  porosity  is  only  —2%  lower  than  that  of  pure  mullite. 
For  moderate  alumina  content,  XA  £  0.2,  the  porosity  remains 
extremely  stable  during  aging,  within  0.5%  of  E,,.  For  higher 
values  (XA  >  0.3),  the  porosity  diminishes  during  aging.  The 
transition  occurs  near  the  expected  percolation  threshold,  i.e..  for 
monosized  spheres  with  random  dense  packing,  the  threshold  is 
predicted  to  be  at  —0.25. 19  Interestingly,  even  beyond  the  thresh¬ 
old  (XA  >  0.25),  the  porosity  change  occurs  remarkably  slowly. 
For  example,  for  XA  =  0.4,  the  porosity  reduction  is  —6%  of  P„ 
after  a  1000  h  aging.  By  comparison,  pure  alumina  attains  this 
porosity  reduction  after  only  10  min  at  the  same  temperature  (inset 
of  Fig.  2(b)).  The  results  demonstrate  that  the  mullite  network  is 
effective  in  inhibiting  densification  even  for  relatively  large 
amounts  of  sinterable  phase,  significantly  beyond  the  percolation 
threshold. 

For  each  aging  treatment,  the  modulus  increases  approximately 
linearly  with  XA  (Fig.  6(a))  at  a  rate  characterized  empirically  by 

£  =  EM(1+PXA)  (12) 

where  p  2-3.  The  increase  is  caused  in  part  to  the  difference  in 
particle  moduli:  £M  =  200  GPa  vs.  EA  =  400  GPa.  But  a 
rule-of-mixtures  prediction  based  on  modulus  difference  alone 
yields  a  lower  modulus  increase,  notably: 

Ea 

p=  —  -1  =  1  (13) 

E  M 

The  disparity  is  attributed  to  the  larger  junction  sizes  at  the 
alumina-alumina  and  alumina-mullite  contacts,  relative  to  mullile- 
mullite,  because  of  the  more  rapid  sintering  of  alumina.  Both 
features  are  incorporated  in  the  model  presented  in  the  following 
section.  Similar  empirical  trends  are  obtained  in  toughness,  char¬ 
acterized  by 

f  =  f„(l  +  yXA)  (14) 

where  7  =  2. 

The  fracture  surfaces  of  the  mixed  particle  compacts  are  similar 
to  those  of  pure  mulllite  in  that  there  is  little  evidence  of  broken 
junctions  even  after  the  longest  aging  time.  Some  examples  are 
shown  in  Figs.  4(c)  and  (d). 

(2)  Modeling 

The  preceding  models  for  junction  growth  and  associated 
property  changes  in  monophase  aggregates  have  been  extended  to 


Fig.  6.  Effects  of  composition  and  aging  time  on  (a)  Young's  modulus 
and  (b)  toughness.  Solid  lines  are  model  predictions  (Eqs.  (16)  and  (18)). 
with  X  =  2,  t)  =  3  and  'P  =  1. 


mullile-alumina  mixtures.  Here,  three  junction  types  are  present: 
alumina-alumina  (AA),  alumina-mullite  (AM),  and  mullite-mullite 
(MM).  The  modeling  approach  involves  an  averaging  of  the 
pertinent  junction  property,  weighted  by  the  number  fraction  of  the 
associated  junction  type,  to  obtain  the  corresponding  aggregate 
property.  For  simplicity,  the  alumina  and  the  mullite  particles  are 
assumed  to  be  the  same  size  and  arranged  randomly  in  the  mixture. 
A  rudimentary  statistical  analysis  yields  the  junction  fractions./: 


> 

> 

II 

>< 

( 1 5«) 

/mm  —  (1  XxY 

(156) 

/am  =  2Xa(  1  —  XA) 

(15c) 

A  model  for  the  Young’s  modulus  of  the  two-phase  mixtures 
has  been  derived  following  an  approach  similar  to  that  of  Walton1’ 
for  monophase  aggregates,  making  appropriate  changes  to  reflect 
differences  in  junction  characteristics:  (i)  The  description  of 
junction  stiffness  is  modified  to  account  for  the  moduli  of  the 
particles  on  either  side  of  the  junction.  For  dissimilar  particle 
junctions,  the  junction  stiffness  in  the  Hertzian  limit  is  proportional 
to  2A/(  I  +  X)  where  X  =  EA/EM.20  (ii)  The  junction  stiffness  is 
further  proportional  to  the  junction  radius.  In  light  of  the  relative 
sintering  rates  of  alumina  and  mullite,  the  radii  of  the  three 
junction  types  are  expected  to  rank  as  r;AA  >  aAM  >  aMM.  (iii)  The 
area  of  the  alumina-mullite  junction  is  assumed  to  be  the  average 
of  the  areas  of  the  alumina-alumina  and  mullite-mullite  junctions. 
This  result  can  be  expressed  in  nondimensional  form  as  «AM/«M\i 
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=  V(l  +  t))/2  where  -r)  is  a  junction  area  ratio,  defined  by  t)  = 
(aAA/nMM)2.  (iv)  The  modulus  is  determined  from  the  mean 
junction  stiffness,  weighted  by  the  fraction  of  the  three  junction 
types,  given  by  Eq.  (15).  Both  the  arithmetic  mean  and  the 
harmonic  mean  are  considered. 

Using  the  arithmetic  mean  of  the  junction  stiffness,  the  resulting 
aggregate  modulus  is  predicted  to  be: 

E  ,  1(1  +  q)  /  2X  \ 

Ym  =  ***  N/n  +  2XA(1  -  xA)  J— J—  +  (1  -  xA)2 

(16) 


The  three  additive  terms  on  the  right  side  represent  the  contribu¬ 
tions  from  alumina-alumina,  alumina-mullite,  and  mullite-mullite 
junctions,  respectively.  Analogously,  the  modulus  prediction  ob¬ 
tained  from  the  harmonic  mean  of  the  junction  stiffness  is: 


E 


^  +  2Xa(1  —  XA) 


1  +  X 


I  +  q  \  2X 


+  0  -*a)2 


(17) 


A  preliminary  assessment  of  these  results  has  been  made 
through  comparisons  with  numerical  simulations  based  on  the 
discrete  element  method.  The  nature  of  the  simulations  is  essen¬ 
tially  the  same  as  those  used  for  determining  toughness  of  particle 
aggregates,  with  the  exception  that  the  aggregate  now  consists  of 
two  particle  types  with  elastic  properties  given  by  those  of  alumina 
and  mullite  (X  =  2).  In  the  present  simulations,  all  junctions  are 
assumed  to  be  of  the  same  size  (ri  =  I ).  Figure  7  shows  the  results 
for  two  sets  of  simulations,  for  junction  radii  a/R  =  0.053  and 
0.313.  The  results  are  presented  in  terms  of  a  modulus  parameter 
that  varies  from  0  for  XA  =  0  to  1  for  XA  =  1 .  Also  shown  are  the 
predictions  from  the  analytical  models  (Eqs.  (16)  and  (17)).  Over 
the  entire  range  of  XA,  the  predictions  from  the  two  analytical 
models  bound  the  simulation  results.  In  the  limit  where  the 
concentration  of  the  stiffer  particles  approaches  zero  (XA  — »  0  in 
this  case),  the  arithmetic  mean  method  yields  more  accurate 
results.  By  contrast,  in  the  limit  where  XA  — »  1,  the  harmonic  mean 
method  is  more  accurate.  For  the  range  of  compositions  considered 
in  this  study  (XA  =  0-0.4).  the  arithmetic  mean  method  is 
preferred  and  is  used  in  subsequent  calculations. 

An  analogous  analytical  model  has  been  developed  for  the 
toughness  of  a  mixed  aggregate.  Following  Eq.  (9),  the  contribu¬ 
tion  from  each  junction  type  is  assumed  to  be  proportional  to 
Tfa/R)2,  weighted  by  the  corresponding  number  fraction  (Eq.  ( 15). 
The  junction  toughnesses  are  denoted  FAA,  FMM  and  TAM;  the 
alumina  to  mullite  toughness  ratio  is  defined  as  'F  =  TAA/rMM 


Alumina  content,  XA 


Fig.  7.  Numerical  simulations17  and  analytical  results  (Eqs.  (16)  and 
(17))  for  Young’s  modulus  of  a  mixture  of  spherical  monosized  alumina 
and  mullite  particles  (X  =  2  and  q  =  1 ). 


and  the  alumina-mullite  junction  toughness  is  taken  to  be  the 
average  of  the  alumina-alumina  and  mullite-mullite  junctions. 
Using  the  arithmetic  averaging  method,  the  toughness  is  predicted 
to  be: 

f  ,  /l+q\/l+'F\ 

Y  =  +  2*a(  I  -  **)(— J  +  ( 1  "  X'y 

(18) 

The  models  were  implemented  in  the  following  way.  The  modulus 
Em  of  the  pure  mullite  aggregate  was  obtained  from  Eq.  (7). 
calibrated  according  to  the  procedures  outlined  above.  Based  on 
the  known  elastic  moduli  of  the  panicles,  X  =  2.  In  comparing  Eq. 
(15)  to  the  experimental  data,  the  remaining  unknown  variable,  q. 
was  used  as  a  fitting  parameter.  The  results  of  this  fitting  arc 
shown  in  Fig.  6(a).  The  procedure  yields  q  =  3  i  I  for  all  aging 
times.  This  value  is  qualitatively  consistent  with  the  expectation 
that  the  alumina-containing  junctions  sinter  more  rapidly  than 
those  with  only  mullite.  and  hence  their  junctions  arc  larger.  A 
similar  procedure  was  used  to  compare  the  toughness  model 
predictions  with  the  experimental  measurements.  For  consistency, 
the  value  of  q  inferred  from  the  modulus  measurements  was  used 
in  Eq.  (18).  The  remaining  unknown  variable.  M7,  was  used  for 
fitting  (Fig.  6(b)).  The  procedure  yields  a  junction  toughness  ratio 
■'F  =  1.0  ±  0.3.  This  result  suggests  that  the  toughnesses  of  the 
three  junction  types  are  similar  to  one  another  and  that  the  increase 
in  aggregate  toughness  with  alumina  content  is  mainly  caused  by 
the  increase  in  the  average  junction  area.  The  quality  of  the  fits  to 
both  the  modulus  and  the  toughness  data  for  all  three  aging  times 
appears  satisfactory,  especially  in  light  of  the  scatter  in  the 
experimental  data.  Furthermore,  it  is  recognized  that  the  scatter 
precludes  determination  of  the  fitting  parameters  q  and  ^  with 
high  accuracy. 


VI.  Implications  for  Crack  Deflection 

Both  the  experimental  results  and  the  model  predictions  have 
been  used  to  calculate  the  variation  in  the  crack  deflection 
parameter.  1.  with  aging  lime,  assuming  the  mixtures  to  be 
combined  with  Nextel  720  fibers  (3M  Corp.)  to  form  a  composite. 
The  results  are  plotted  on  Fig.  8.  Over  the  entire  range  of 
parameter  values  considered  in  this  study,  the  crack  deflection 
condition  is  predicted  to  be  satisfied:  that  is,  I  >  I.  The 
experimental  results  are  also  plotted  on  Fig.  I.  demonstrating  that 


Fig.  8.  Effects  of  aging  time  and  composition  on  the  crack  deflection 
parameter  S.  using  fiber  properties  Ij  =  15  J/nr  and  =  260  (iPa.  Solid 
lines  are  model  predictions.  Extrapolations  of  the  predictions  to  X  -  u> 

1  yield  the  critical  limes,  q,  at  which  crack  penetration  is  predicted  to 
occur. 
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the  property  combinations  lie  within  the  domain  of  crack  deflec¬ 
tion  in  the  He  and  Hutchinson7  diagram.  These  results  reaffirm 
mullite-alumina  mixtures  as  good  candidate  matrixes  for  use  in 
oxide  composites,  even  with  rather  high  alumina  content. 

On  extrapolation  of  the  predictions  in  Fig.  8  to  2,  =  o>  =  I,  the 
critical  aging  times  tc  for  crack  penetration  are  obtained.  The  times 
range  from  4000  h  for  60M/40A  to  60  000  h  for  pure  mullite: 
comparable  to  the  targeted  service  times  for  CFCC  components. 
Based  on  the  differences  in  /c,  matrixes  with  lower  alumina 
contents  would  be  preferred. 


IV.  Conclusion 

Mullite-alumina  mixtures  exhibit  physical  and  mechanical  char¬ 
acteristics  that  make  them  attractive  for  use  as  porous  matrixes  in 
ceramic  composites.  The  key  attributes  include  stable  porosity  at 
the  targeted  service  temperatures  of  oxide  composites,  especially 
for  alumina  concentrations  <20%,  and  combinations  of  modulus 
and  toughness  that  should  lead  to  deflection  of  a  matrix  crack  at  a 
fiber-matrix  interface.  The  models  of  mechanical  properties  of 
particle  aggregates  presented  here  coupled  with  the  crack  deflec¬ 
tion  parameter  provide  a  means  of  extrapolating  the  experimental 
data  to  longer  aging  times,  for  the  purpose  of  determining  the 
critical  times  at  which  damage  tolerance  will  be  lost.  With 
knowledge  of  the  activation  energies  of  the  sintering  mechanisms, 
the  models  could  be  extended  to  make  such  predictions  for  other 
temperatures. 

Design  of  porous  matrixes  requires  consideration  of  other 
properties,  including  shear  and  interlaminar  strength,  as  well  as 
compressive  strength.  These  generally  increase  as  the  matrix  is 
strengthened.  As  a  result,  among  the  mixtures  considered  in  this 
study,  those  with  higher  alumina  contents  would  be  preferred,  in 
contrast  to  the  mullite-rich  mixtures  selected  based  on  crack 
deflection  considerations.  These  issues  are  currently  being 
addressed.2 
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An  analytical  model  to  predict  the  effective  conductivity  of  a  partially  sintered  solid  is  presented.  It 
is  based  on  an  analysis  of  heat  flow  through  a  spherical  particle  embedded  in  a  self-consistent 
matrix.  Determination  of  the  effective  conductivity  reduces  to  solving  an  integral  equation,  and  a 
numerical  solution  of  this  equation  is  presented.  In  the  limiting  cases  wherein  the  contact  size  is 
either  small  or  comparable  to  the  particle  size,  approximate  analytical  results  are  obtained  which 
agree  very  well  with  the  numerical  solution.  Furthermore,  when  pieced  together,  the  two  solutions 
are  found  to  describe  the  effective  conductivity  over  the  entire  range  of  contact  sizes.  ©  2004 
American  Institute  of  Physics.  [DOl:  10. 1 063/ 1.1630360] 


I.  INTRODUCTION 

There  has  been  a  considerable  amount  of  research  di¬ 
rected  toward  predicting  the  effective  thermal  conductivity 
of  porous  and  granular  solids.1"10  In  most  of  these 
treatments,3-7  it  has  been  assumed  that  contacts  between  par¬ 
ticles  are  sufficiently  small  so  that  the  constriction  resistance 
to  heat  flow  corresponds  to  a  contact  spot  separating  two 
infinite  half  spaces.11  This  treatment  is  reasonable  in  the 
analysis  of  granular  solids  where  the  contact  areas  depend 
upon  Hertzian  mechanics.  However,  in  the  case  of  partially 
sintered  solids,  the  finite  contact  size  must  be  taken  into  ac¬ 
count.  Although  there  have  been  two  numerical  studies,910 
this  problem  has  received  less  attention. 

Here,  a  calculation  of  the  self-consistent  type,  where  the 
heat  flow  through  a  sphere  embedded  in  a  self-consistent 
matrix,  is  presented.  No  other  shape  for  equiaxed  particles 
seems  to  offer  any  chance  of  an  analytical  solution.  Similar 
to  the  calculation  done  by  Argento  and  Bouvard  9  we  take 
two  contacts  per  sphere.  These  two  caps  are  held  at  a  self- 
consistent  temperature,  and  the  remaining  portion  of  the 
sphere  surface  is  adiabatic. 

II.  PRELIMINARY  COMMENTS 

From  a  macroscopic  viewpoint,  the  porous  solid  is  iso¬ 
tropic  with  respect  to  heat  flow  and  its  effective  conductivity 
is  ke .  Consequently,  only  heat  flow  in  one  direction  needs  to 
be  analyzed.  Here,  heat  flow  in  the  Z  direction  is  considered, 
and  the  apparent  temperature  is  written  as  %—CZ.  On  a 
finer  scale,  the  microstructure  is  nonuniform,  and  flow  pat¬ 
tern  is  tesellated.  In  the  simplest  case,  there  are  two  phases: 
Solid  particles  with  thermal  conductivity  k ,  and  a  gas  phase 
with  thermal  conductivity  kg .  A  typical  particle  has  a  portion 
of  its  surface  surrounded  by  gas  with  the  remainder  of  its 
surface  being  in  contact  with  its  nearest  neighbors.  Let  p  be 
equal  to  the  ratio  of  the  contact  radius  to  solid  particle  radius. 
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For  sufficiently  small  values  of  p,  the  flow  constriction 
through  the  solid  contacts  becomes  severe  and  the  heat  How 
is  predominantly  through  the  gas  phase.  On  the  other  hand,  it 
has  been  shown  by  Batchelor,4  that  when  kp>kK ,  the  heat 
flow  through  the  gas  phase  is  negligible.  This  condition  is 
satisfied  in  most  partially  sintered  solids. 

Initially,  the  individual  solid  particles  are  approximately 
spherical.  However,  as  sintering  proceeds,  the  shape  of  each 
particle  becomes  irregular.  If  the  particles  remain  equiaxed,  it 
is  useful  to  consider  a  “test”  region  that  contains  an  average 
representative  particle.  The  analysis  presented  here  is  con¬ 
cerned  with  analyzing  the  temperature  field  in  this  represen¬ 
tative  test  region,  and  then  relating  this  field  to  obtain  an 
estimate  for  the  effective  conductivity.  In  a  calculation  of  this 
type,  neither  the  particle  radius  a  nor  the  magnitude  C  of  the 
thermal  gradient,  affects  the  estimate  for  k,. .  Consequently,  it 
is  convenient  to  define  dimensionless  variables  z  —  Z/a.  and 
Te=T(,/Ca.  Using  these  definitions, 


From  this  point  onward  in  the  article,  all  spatial  variables  and 
temperatures  are  dimensionless. 


III.  TEMPERATURE  WITHIN  THE  TEST  PARTICLE 

To  make  the  problem  tractable,  the  shape  of  the  test  par¬ 
ticle  is  taken  as  being  a  sphere  having  a  radius  equal  to  a.  In 
view  of  the  spatial  scaling,  the  sphere  has  a  dimensionless 
unit  radius,  and  the  dimensionless  temperature  in  the  particle 
is  designated  by  T.  For  simplicity,  the  origin  is  placed  at  the 
center  of  the  particle.  Spherical  coordinates  (/-,  0)  are  used 
where  r  is  the  distance  from  the  center  of  the  sphere  and  0  is 
the  angle  measured  from  the  r  axis.  The  temperature  T  must 
satisfy  Laplace’s  equation 

V3T=  0.  (2) 
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(9> 

The  relations  in  Eqs.  (6)  and  (7)  are  referred  to  as  triple 
series  relations.12'13  The  adjective  “triple”  arises  from  the 
fact  that  there  are  three  different  zones:  (i)  a>0s=O,  (ii) 
j. 3>0>a ,  and  (iii)  tt^0>/3.  These  mixed  boundary  value 
problems  are  generally  more  difficult  than  the  standard  prob¬ 
lems  where  there  is  only  a  single  condition  over  the  entire 
surface.  The  aim  is  to  find  the  coefficients  C„  that  satisfy  the 
triple  relations. 

Owing  to  the  symmetry  of  the  present  problem,  the  tem¬ 
perature  is  antisymmetric  with  respect  to  the  plane  z  =  0  so 
that  the  Legendre  coefficients  C0.  C2>  C4,--C2n  must  van¬ 
ish.  As  described  by  Collins,12  in  order  to  determine  the  Leg¬ 
endre  coefficients,  it  is  necessary  to  write  Cn=A„  +  Bn  and  it 
follows  that 


FIG.  1.  Spherical  test  particle  where  the  polar  caps  are  symmetrical,  i.e., 
f)-Tr-a  so  that  the  temperature  is  an  odd  function  of  z. 

The  test  particle  is  shown  in  Fig.  1,  and  on  the  surface,  r 
=  1,  the  following  boundary  conditions  are  imposed.  Within 
the  equatorial  region,  i.e.,  a<0</3 ,  there  is  no  radial  flow 
and  the  boundary  condition  is  then 

0T 

—  =0  U3>0>a).  (3) 

The  apparent  overall  temperature  T(.  =  z. ,  and  the  surface  of 
the  sphere  corresponds  to  c  =  cos  6.  Therefore,  on  the  north 
and  south  polar  caps  subtending  the  cone  of  angle  2a,  the 
boundary  condition 

T=cos0  (a>  05=0, 7r  3=  0> y3),  (4) 

is  imposed.  The  temperature  in  the  test  sphere  is  written  in 
the  form 


C„=An  +  Bn  and  B,,=  -(-  1  )"An .  (10) 

The  triple  relations  in  Eqs.  (6)  and  (7)  are  written  as 

x 

E  (A  „  +  /?„)(  1  +H„)P„(cos  0)  =  O  (a>0s 0),  (II) 

W  =  0 

X 

E  A„(  1  +  2/i)/>„(cos  0)  —  SXT’ifcos  0)  =  O  {ttS*  0>  a), 

n  =  0 

(12) 

X 

E  B„(  1  +  2n)P„(co$  <?)  — 3\P|(cos  0)  =  O  (/3>03= 0), 

/i  =  0 

(13) 

X 

E  (A,l  +  B„)(l+H„)Pn(cos0)  =  O  (t t^0>P).  (14) 

77  =0 


l+2« 

7'=rP,(cos  0)- E  - C„P„(cos  0)/-",  (5) 

n=  I  ft 

where  the  terms  P„( cos  0)  are  Legendre’s  polynomials.  The 
temperature  is  an  odd  function  of  r.  so  that  only  the  terms, 
n=  1,3,5,- ••  are  taken.  The  expression  for  T  satisfies 
Laplace’s  equation,  and  the  coefficients  C„  are  found  using 
the  two  boundary  conditions  that  are  prescribed  above.  Using 
these  two  boundary  conditions,  it  is  found  that  the  coeffi¬ 
cients  C„  must  satisfy  the  following  conditions: 

X 

E  C„(  I  +H„)P„(cos  0)  =  O  (a>0 5=0,773=  0>/3) 

11=  I 

(6) 

X 

E  C„(  1  +2n)P„{cos  0)  =6\/>i(cos  0)  (/ 3>0>a ), 

n  =  I 

(7) 

where  the  various  coefficients  are  defined  by 


The  coefficients  A0,  B0  are  not  equal  to  zero  so  that  the 
lower  limit  on  the  summation  has  been  changed.  To  avoid 
objection  over  the  possibility  of  division  by  zero,  the  term 
H0=  1;  this  definition  is  inconsequential  because  A0+#o 
=  0.  Owing  to  the  symmetry  of  the  problem,  only  the  first 
two  equations  [Eqs.  (II)  and  (12)]  need  to  be  considered. 

It  is  known  that  the  Legendre  polynomials  can  be  repre¬ 
sented  by  the  Mehler-Dirichlet  integrals 


\/2  f# 

P „(cos  0)=  — — 

17  J i=0 


cos[(n+  1/2  )t\dt 
V cos r  —  cos  0 

sin[(n  +  1/2  )r]dt 
ycos  0-cost 


(15) 


(16) 


The  representation  for  P„( cos  0)  in  Eq.  (15)  is  used  in  Eq. 
(II)  and,  after  interchanging  the  order  of  integration  and 
summation,  it  follows  that 
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/, 


®  dt  ^ 

7 .  J  2  (An  +  B„)(\+Hn) 

o  Vcos/  — cos  6  [  «=o 


Xcos (n+  1/2 )t  f  =0. 


(17) 


Similarly,  the  representation  for  P„(cos  0)  in  Eq.  (16)  is  used 
in  Eq.  (12)  and,  after  interchanging  the  order  of  integration 
and  summation,  it  follows  that 


[ 


dt 


=0sjcos  8— cost 
d 


X 


dt 


2  Ancos(n  +  1/2)/- \ cos3//2 

n  =  0 


=  0. 


(18) 


Equations  (17)  and  (18)  can  be  solved  using  Eqs.  (2.3.5)  and 
(2.3.6)  in  Ref.  13.  In  both  cases,  the  quantities  within  the 
{•••}  must  vanish.  In  Eq.  (18),  the  portion  within  [•••]  =  £> 
and  the  constant  D  =  0;  this  follows  by  noting  that  the  ex¬ 
pression  within  [•••]  vanishes  at  /  =  7r.  Thus,  the  two  results 

oc 

2  (An  +  Bn)(  1  +  W„)cos[(rt+  I/2)r]  =  0  («>/s=  0), 

n  =  0 

(19) 


2  /4„cos[(n+  l/2)/]-\ cos(3//2)=0  (-n-s=/>a), 

n  =  0 

(20) 

are  obtained.  Here,  the  coefficient  A„  is  found  and  once  this 
is  known  the  other  two  coefficients,  C„ ,  B„  can  be  found  by 
use  of  Eq.  (10).  To  find  A„,  an  unknown  function  j(t )  is 
defined  as  follows 


2  4„cos| 
11  =  0 


/  n 

| 

[\  +  2/  J 

=1 

[  \  cos(3//2)  -j(t)/2  (a>/s=0) 

[\cos(3//2)  {tr^t>a) 

(21) 

The  functions  cos(n+l/2)f  are  orthogonal  over  the  interval 
(0=sr«7r)  and  thus  the  coefficient  A„  is 


f" 

7  l  \ 

j(u) cos 

K ”  +  2/  ] 

J  11  =  0 

du. 


(22) 


where  <5 ,  „  =  0  when  n  =/=  1.  Expressions  (21)  and  (10)  are  then 
substituted  into  Eq.  (19)  to  obtain 


m  +  2'Z  A„{(-  1)"  — ff „[!-(- 1)"]}C0S 

n  =  0 


f  l 
n+-z\t 
l  1 


—  2  \  cos 


3/' 

T.  ' 


(23) 


Using  the  expression  for  A,,,  and  interchanging  the  order  of 
summation  and  integration,  Eq.  (24)  is  found  for  the  deter¬ 
mination  of  j(t): 

1  f  a 

j(t)  =  cos(3//2)-l —  j(u)J(uj)du 

tt  Jo 

where  (a>t5?0).  (24) 


The  kernel  J(ti,t)  =  K(t  —  u)  +  K(t  +  u)  and  K(^)  is 


^(f)=E  {( —  I H„[  I  -  ( -  1  )"]}cos 

n  =  0 


2  sec[  ~  j  +  7T  sin: 


+  cos|  —  I  In 


tan- 1  - 


(25) 


The  closed  forms  for  the  infinite  summations  in  Eq.  (25)  can 
be  found  from  the  results  given  in  standard  tables.14  For 
general  values  of  a ,  the  function  j(t )  can  only  be  found  by  a 
numerical  treatment. 

To  summarize  this  section,  the  function  j(t)  is  a  solution 
of  the  integral  Eq.  (24).  The  coefficients  An  are  then  found 
from  Eq.  (22)  and  the  odd  Legendre- Fourier  series  coeffi¬ 
cients  C2n+ 1  =2A2ii+  j  .  Finally,  these  coefficients  are  used 
in  Eq.  (5)  to  find  the  temperature. 


IV.  POWER  EQUIVALENT  CONDUCTIVITY 


On  the  macroscopic  level,  the  dimensionless  apparent 
temperature,  T,.=  z,  is  given  in  Eq.  (I).  The  thermal 
“power”  P  (which  here  has  the  same  units  as  the  conductiv¬ 
ity),  expended  within  the  test  volume  V  is  expressed  by  the 
integral 

=  gradrr-grad7'(,£/F  =  /:),V/. 

The  test  volume  V  is  comprised  of  two  phases  and,  if  there  is 
negligible  flow  through  the  gas  phase,  then  P  is  due  to  the 
power  expended  in  the  solid  sphere.  The  power  P  can  be 
alternatively  written  as 

P  =  k\  grad  T  ■  grad  Td  V, 

J  v, 

where  Vv  =  4 7r/3  is  the  dimensionless  volume  of  the  sphere. 
Upon  equating  these  two  expressions  for  P,  the  power 
equivalent  thermal  conductivity  is 

kt,=  kfa,  (26) 

where  /=  Vs/V  is  the  volume  fraction  of  solid,  and  the  con¬ 
ductance  (j  is  defined  as 


grad  T ■  grad  TdV. 


(27) 


The  conductance  a  is  equal  to  the  average  squared  magni¬ 
tude  of  the  temperature  gradient  within  the  sphere.  In  this 
expression,  the  volume  integral  can  be  converted  to  an  inte¬ 
gral  over  the  surface  of  the  sphere  by  the  use  of  the  diver¬ 
gence  theorem.  In  view  of  the  boundary  conditions  in  Eqs. 
(3)  and  (4),  along  with  the  orthogonal  nature  of  the  Legendre 
polynomials,  it  follows  after  substituting  the  series  in  Eq.  (5) 
into  the  surface  integral  that  cr=  1  —  3C).  The  coefficient 
C|  =  2/4]  and  using  the  expression  for  4,  in  Eq.  (22),  the 
expression  for  a  is  found  as 

6  f“  /3k\ 

a--—  y(«)cos  —\du.  (28) 

tt  J  n  =  o  \  2  / 


Downloaded  20  Feb  2004  to  128.1 11.74.93.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://jap.aip.org/jap/copyright.jsp 


J.  Appl.  Phys.,  Vol.  95,  No.  1,  1  January  2004 


J.  R.  Dryden  and  F.  W.  Zok 


159 


As  discussed  later,  this  represents  the  half-space  approxima¬ 
tion,  and  is  shown  in  Fig.  2. 

B.  Approximate  solution  for  cr  when  a>1 

For  larger  values  of  a,  the  integral  equation  is  not  ame¬ 
nable  to  an  iterative  solution  of  the  type  described  above.  On 
the  other  hand,  when  a  becomes  large,  the  perturbation  of 
the  heat  flow  from  the  uniform  flux  is  not  very  severe  and  it 
is  possible  to  find  fairly  accurate  approximate  solutions. 

To  begin,  it  is  known  from  the  theory  of  variational  cal¬ 
culus  that  a  lower  bound  for  cr  can  be  obtained  by  suitable 
mathematical  manipulation  of  the  functional  given  in  Eq. 
(27).  To  obtain  this  bound  it  is  necessary  to  find  a  vector,  say 
s,  that  has  its  divergence  equal  to  zero,  and  also  has  no  com¬ 
ponent  crossing  r—  I  in  the  region  /3>0>a.  This  vector  can 
be  obtained  as  the  gradient  of  a  function  F  and  the  function 
is  harmonic  so  that  V2F=  0.  Such  a  function  can  be  ex¬ 
pressed  as  a  Legendre-Fourier  series 


FIG.  2.  The  conductance  a  as  a  function  of  the  cap  angle  a.  The  symbols  O 
are  obtained  by  numerical  integration.  [The  solid  line  labeled  “small  a”  is 
obtained  from  Eq.  (30).  The  results  of  Eq.  (31)  represent  the  half-space 
approximation.  Finally,  the  dashed  line  labeled  “large  or”,  represents  the 
prediction  of  Eq.  (36)]. 


To  evaluate  this  conductance,  the  function  j(r)  is  found  by 
numerically  solving  the  integral  Eq.  (24)  for  the  following 
ten  values  of  a:  7r/20,  277/20,  377/20,-  ■  ■ ir/2 .  For  each  value  of 
a,  the  region  0<u<a  is  subdivided  into  60  segments  (this  is 
found  to  give  sufficiently  accurate  results),  and  a  set  of  60 
linear  equations  is  generated  using  Eq.  (24).  This  set  of  60 
equations  is  then  solved  to  find  y'(r)  at  the  discrete  incre¬ 
ments,  r  =  /| ,r2,- within  the  interval  0<t<a.  The  val¬ 
ues,  jOOiKh)’'  ”7(f6o)>  are  then  used t0  evaluate  crand  the 
results  are  shown  by  the  symbols  in  Fig.  2. 


4n  +  3 
2/2+  1 


hn+\P  2n+  1 


(cos  8)r 


2/i  +  I 


(32) 


The  coefficients  / 2„+i  are  chosen  so  that  F  satisfies  the 
boundary  conditions 


cos  8 

a>8>  0 

0 

/3>8>a. 

(33) 

cos  8 

77  >  8>/3 

To  satisfy  this  boundary  condition,  the  coefficients  /?„+  ,  are 
expressed  by  the  integrals 


'2/?+  I 


'a 

=  ^2/i+l 

JO 


(cos  #)cos  6  sin  del 8. 


(34) 


A.  Approximate  solution  for  a  when  o<1 

If  a  is  small,  the  integral  equation  is  amenable  to  an 
iterative  solution.  Using  the  expression  for  j(t)  given  in  Eq. 
(24),  the  conductance  can  be  rewritten  as 


cr= 


3 

17 


sin  3  a 
aT-y-J 


+  —  7(w,r)/(w)cosfyj  r/nrfr.  (29) 

When  a  becomes  sufficiently  small,  the  function  j(u) 
=  cos(3k/2),  and  the  kernel  can  be  approximated  by  its  sin¬ 
gular  portion  7(u,f)~In|ir-f2|/2.  Neglecting  the  higher- 
order  terms,  the  conductance  is 


o 


3 
7 r 


a+ 


sin  3  a 
3 


+ 


(30) 


As  shown  in  Fig.  2,  this  estimate  for  cr  is  reasonably  accurate 
for  values  of  a<  1.  In  the  limit  when  a<  1,  the  estimate  for  cr 
reduces  to  linear  equation 


6a 

a - .  (31) 

7 7 


According  to  standard  treatment,15  the  lower  bound  for  cr  is 
found  by  suitable  mathematical  manipulation  upon  the  func¬ 
tional  given  in  Eq.  (27).  After  doing  these  manipulations,  and 
using  the  divergence  theorem,  the  bound  for  cr  is  given  by 
the  following  integral  which  is  evaluated  on  the  spherical 
surface  r  =  1 , 


-/■ 


( 6  T(,  —  3  F) cos  6  sin  Od  0 , 


=  6/,-3E 

h  =  0 


4/i  +  3 


2/2+1 


^2<l+f 


(35) 


When  a— ’77/2,  the  summation  can  be  lepresented  by  its  first 
term  and  the  estimate  for  cr  reduces  to  the  simple  expiession 

cr~  1  —  cos6  a.  (36) 


This  bound  can  be  improved  by  multiplying  F  by  a  constant 
say  D  and  then  adjusting  this  constant  so  as  to  maximize  cr. 
However,  when  a=»7r/2,  this  refinement  does  not  make  much 
improvement  to  the  estimate  for  cr.  The  behavior  of  this 
bound  is  shown  in  Fig.  2  and,  for  large  values  of  a,  it  is  in 
very  good  agreement  with  the  numerical  results. 
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V.  DISCUSSION  AND  CONCLUSIONS 

In  a  real  material,  each  particle  is  surrounded  on  average 
by  N  nearest  neighbors  and  a  is  the  average  cap  angle  of  the 
contacts.  According  to  their  outward  normal,  these  contacts 
can  be  partitioned  into  six  different  directions  ±(i,j,k).  In 
the  analysis  presented  here,  the  cap  angle  a  is  expressed  by 
a=5/V/6.  When  the  contacts  are  very  small  in  comparison 
with  the  size  of  the  particles,  the  conductance  is  given  by  Eq. 
(31),  and  cr=  aNI tr.  The  effective  conductivity  is  then  writ¬ 
ten  as 


This  is  the  same  result  for  the  effective  conductivity  that  has 
previously  been  deduced  by  Batchelor  and  O’Brien.4  As  the 
cap  angle  a  increases,  the  conductance  becomes  a  nonlinear 
function  of  a  and  cr  is  given  approximately  by  either  Eq.  (30) 
or  (36). 

The  present  analysis  can  be  related  to  a  unit-cell  treat¬ 
ment  where  the  particles  are  arranged  in  a  simple  cubic  array. 
When  a<€  1,  the  spherical  particles  are  barely  touching  each 
other  and  the  volume  fraction,  f=ir/6,  then  corresponds  to 
simple  cubic  packing  of  spheres.  Furthermore,  according  to 
Eq.  (31),  the  limiting  behavior  of  the  conductance  is  given 
by  o — 6a/w  so  that  the  effective  conductivity  is  ke  =  kfa 
~kea.  Now,  consider  the  vertical  flow  of  heat  through  a 
simple  cubic  array.  Taking  symmetry  into  account,  the  four 
lateral  faces  are  adiabatic  while  the  top  and  bottom  cube 
faces  are  isothermal.  Thus,  when  1,  each  isothermal  con¬ 
tact  has  a  radius  approximately  equal  to  a,  and,  neglecting 
the  curvature  of  the  sphere,  the  total  constriction  resistance  is 
/?,.=  1/2 ka.  The  unit-cell  cube  has  dimensions  2X2X2,  and 
if  its  conductivty  is  ke  then  its  resistance  is  Re=  1/2 k,. .  Upon 
equating  these  two  resistances,  it  follows  that  kv  =  ka  and 
this  is  in  agreement  with  the  treatment  given  here. 

As  sintering  progresses  the  size  of  the  necks  between 
neighboring  particles  increases.  In  the  unit-cell  model  the 
particle  shape  has  cubic  symmetry;  the  initially  spherical  par¬ 
ticle  acquires  six  facets  corresponding  to  the  cube  faces. 
There  is  no  analytical  method  for  deducing  the  thermal  field 
in  such  a  truncated  shape,  and  it  is  necessary  to  resort  to  a 
numerical  solution.  Selecting  a  spherical  test  particle  has  the 
advantage  that  it  allows  a  semianalytic  solution  for  the  ther¬ 
mal  field. 

The  effective  conductivity  in  Eq.  (26)  has  the  simple 
form  ke  =  kfcr,  and  up  to  this  point,  only  the  behavior  of  cr 
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has  been  considered.  The  conductance  cr  depends  on  the  neck 
growth,  whereas  the  volume  fraction /depends  on  the  degree 
of  densification.  Ashby16  points  out  that  there  are  six  distin¬ 
guishable  sintering  mechanisms,  all  causing  neck  growth.  In 
contrast,  only  some  of  mechanisms  lead  to  densification.  Al¬ 
though  obtaining  a  general  relationship  between  a  and  /  is 
beyond  the  scope  of  this  article,  a  few  remarks  are  possible. 
First,  during  the  early  stages  of  sintering  there  is  only  neck 
growth  while  the  /  remains  constant.  Second,  at  later  stages 
of  sintering,  corresponding  in  the  analysis  to  a— >ir/2,  the 
conductance  a~\.  The  effective  conductivity  is  then  ap¬ 
proximately  given  as  kt.  =  kf  and  the  model  reduces  to  the 
law  of  mixtures.  In  almost  fully  dense  sintered  solids,  the 
porosity  exists  as  isolated  spheres  and  the  effective  thermal 
conductivity  is  better  described  using  Maxwell’s17  estimate 
instead  of  the  law  of  mixtures. 
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A  theoretical  framework  is  developed  for  the  notched  strength 
of  ceramic  composites  that  exhibit  rising  fracture  resistance.  It 
is  based  on  established  concepts  of  crack  stability  under 
stress-controlled  loadings.  On  using  a  linear  representation  of 
the  resistance  curve  (expressed  in  terms  of  an  energy  release 
rate),  straightforward  analytical  solutions  are  obtained  for  the 
strength  as  well  the  amount  of  stable  crack  growth  preceding 
fracture  and  the  associated  fracture  resistance.  Calculations 
are  performed  for  several  test  configurations  commonly  used 
for  material  characterization,  including  single-  and  double¬ 
edge-notched  tension,  center-notched  tension,  and  single-edge- 
notched  bending.  The  results  reveal  salient  trends  in  strength 
with  notch  length  and  specimen  geometry.  An  assessment  of 
the  theory  is  made  through  comparison  with  experimental 
measurements  on  an  all-oxide  fiber  composite.  Transitions  in 
the  degree  of  notch  sensitivity  with  notch  length  are  identified 
and  explored.  The  utility  of  the  theoretical  results  both  for 
rationalizing  the  trends  in  measured  notched  strength  and 
for  guiding  experimental  studies  of  notch  sensitivity  is 
demonstrated. 


I.  Introduction 

Crack  growth  in  continuous  fiber  ceramic  composites  (CFCCs) 
is  characterized  by  a  rising  fracture  resistance  beyond  the 
initiation  toughness.  G0,  a  manifestation  of  fiber  bridging  in  the 
crack  wake.  When  the  bridging  fibers  furthest  from  the  crack  tip 
disengage,  the  resistance  reaches  a  steady-state  value,  Gss,  inde¬ 
pendent  of  further  crack  growth.  The  steady-state  bridging  tough¬ 
ness  is  AGh  =  where  u.  is  the  critical  crack  opening 

displacement  to  disengage  the  bridging  fibers,  o»  is  a  measure  of 
the  pullout  stress,  and  (3  is  a  nondimensional  parameter  of  order 
unity.'  The  amount  of  crack  growth  Aoss  needed  to  achieve  steady 
state  is  dictated  by  a  characteristic  bridging  length  scale:  t?ch  = 
EuJa „,  where  E  is  the  plane  strain  Young’s  modulus.12  For 
typical  values  of  the  material  parameters  ( E  =  100  GPa.  u.  = 
30-100  p,m,  and  o0  =  30-100  MPa),  the  bridging  toughness  is 
lO’-IO4  J/m2  and  the  characteristic  length  scale  is  10-100  mm. 
The  implications  are  2-fold,  (i)  Significant  toughening  can  be 
obtained  at  steady  state,  relative  to  the  initiation  toughness,  (ii) 
Because  of  the  large  amount  of  crack  growth  needed  to  achieve 
steady  state,  the  benefits  of  the  rising  resistance  curve  in  composite 
strength  are  not  fully  realized.  That  is.  crack  growth  becomes 
unstable  with  respect  to  the  applied  stress  long  before  steady  state 
is  achieved. 

The  objective  of  this  paper  is  to  address  the  related  issues  of 
crack  stability  and  notch  sensitivity  of  strength  in  ceramic  com¬ 
posites  that  exhibit  a  rising  fracture  resistance,  both  theoretically 
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and  experimentally.  On  the  theoretical  front,  an  analytical  frame¬ 
work  is  developed  for  describing  crack  stability  in  an  orthotropic 
material.  It  is  based  on  comparisons  of  the  energy  release  rates  of 
common  notched  test  configurations  with  the  material’s  fracture 
resistance.  Analogous  treatments  of  notch  sensitivity  have  been 
used  previously  for  other  materials  that  exhibit  rising  fracture 
resistance,  including  metal-,  whisker-  and  fiber-reinforced  ceram¬ 
ics'1'4  and  ceramics  toughened  by  grain  pullout.'1'6  For  our  purpose, 
the  resistance  curve  GR(Azf)  is  assumed  to  be  linear  over  the  range 
of  crack  growth  preceding  instability.  Furthermore,  the  size  of  the 
inelastic  deformation  zone  at  the  crack  lip  is  considered  to  be  small 
in  relation  to  the  geometric  and  the  bridging  length  scales,  i.c.. 
small-scale  yielding  conditions  are  obtained.  Effects  of  inelasticity 
are  thus  embodied  in  the  initiation  toughness.  Elsewhere  the 
material  is  assumed  to  be  elastically  orthotropic.  To  assess  the 
theory,  comparisons  arc  made  with  experimental  measurements  of 
fracture-resistance  curves  and  notch  sensitivity  of  strength  in  a 
porous-matrix  all-oxide  CFCC. 

It  is  emphasized  that  the  theoretical  approach  is  based  on  an 
assumed  form  of  the  resistance  curve,  notably  linear;  no  attempt  is 
made  to  model  the  resistance  curve  itself.  The  latter  problem  falls 
in  the  domain  of  bridging  micromechanics,  a  field  that  has  been 
explored  exhaustively  over  the  past  decade.  Although  micromc- 
chanics  provides  a  more  fundamental  link  between  the  bridging 
phenomenon  and  the  macroscopic  response,  it  requires  numerical 
solution  of  the  integral  equations  governing  the  crack-opening 
displacement  profile,  the  crack-lip  stress  intensity  factor,  and  the 
bridging  traction  law.  As  a  consequence,  it  is  likely  to  find  limited 
use  in  the  engineering  community.  Arguably,  for  the  purpose  of 
predicting  the  strength  of  notched  or  Hawed  structures,  microme¬ 
chanics  is  not  essential.  Instead,  only  a  convenient  functional  form 
for  the  resistance  curve  is  required.  This  phenomenological  ap¬ 
proach  has  proved  to  be  useful  in  addressing  strength  variability  in 
ceramics.5  Furthermore,  as  illustrated  in  the  present  paper,  in  cases 
where  the  resistance  curve  is  linear  (as  it  is  in  the  oxide  composite 
of  interest),  straightforward  analytical  solutions  for  notched 
strength  are  Obtained.  These  solutions  allow  probing  of  the  effects 
of  specimen  type  and  size  on  notch  sensitivity  and  identifying 
transitions  in  behavior. 


II.  Theory  of  Crack  Stability 

To  begin,  it  is  instructive  to  consider  the  case  of  a  specimen  of 
width  W,  containing  a  small  center  notch  of  length  2 «„  and  loaded 
in  tension  perpendicular  to  the  notch  plane  along  one  of  the 
principal  material  axes.  In  the  limit  where  au  <sc  W,  the  energy 
release  rate  is  given  by 

/  ct2tt«\ 

G  =  g(p,X)^  — £-J  (1) 

where  o  is  the  applied  tensile  stress.  g(p,X)  characterizes  the  effect 
of  orthotropy.7,8  and  p  and  X  are  defined  by 


X  = 


£2 

Et 


(2) 
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P  = 


(E,E2yn- 

2(t,2 


(v,2V2I) 


1/2 


(3) 


with  the  subscripts  1  and  2  referring  to  the  principal  material  axes, 
p.  the  shear  modulus,  and  v  the  Poisson’s  ratio.  Composite 
materials  with  balanced  symmetric  fiber  lay-ups  exhibit  cubic 
in-plane  symmetry,  whereupon  X  =  £,/E2  =  1  and  v,2  =  v21. 
Moreover,  for  0  <  p  <  3,  the  function  g(p)  can  be  well- 
approximated  by8 

g(p)  =  1  +0.44(p-  1)  (4) 

Unstable  crack  growth  (at  the  stress  maximum)  occurs  when9 
G  =  Gr  (5a) 


d  G  dGR 
da  ~  dfl 


(5b) 


These  conditions  are  represented  by  the  tangent  construction  plot 
shown  in  Fig.  I .  Here,  GR  is  plotted  against  the  amount  of  crack 
growth,  A  a,  and  the  energy  release  rate,  G,  is  plotted  against  the 
total  crack  length,  a0  +  A  a.  In  this  context,  the  right  side  of  the 
abscissa  is  interpreted  as  A  a  and  the  left  side  as  a0.  The  resistance 
initiates  at  G0  and  subsequently  increases  with  Art  at  a  rate 
characterized  by  the  tearing  modulus:  T  =  dCR/dn.  The  conditions 
in  Eq.  (5)  are  satisfied  when  the  curve  G(Aa+a0)  is  tangent  to  the 
curve  GR(An).  If  the  resistance  curve  GR(An)  is  linear,  the  stability 
of  the  crack  on  initial  growth  is  dictated  by  the  slopes  of  the  two 
curves  at  G0  (Fig.  1(a)).  That  is.  if  dG/drr  >  T  at  G  =  G0,  crack 
growth  occurs  unstably  under  a  decreasing  stress.  Consequently, 
the  strength  is  controlled  by  the  initiation  toughness  via  the 
Griffith  equation;  no  strength  benefit  is  derived  from  the  increas¬ 
ing  resistance.  Conversely,  if  dG/da  <  T  at  G  =  G0,  continued 
crack  growth  requires  an  increase  in  stress  beyond  that  for 
initiation.  The  transition  from  unstable  to  stable  growth  occurs  at 
a  critical  tearing  modulus,  Tc ,  at  which  dG/da  =  T  at  G  =  G0. 
From  Eq.  (5)  or  inspection  of  Fig.  1,  the  critical  tearing  modulus 
in  nondimensional  form  is 


TtW  _  W 
G0  o0 


(6) 


Interestingly,  since  both  Tc  and  G0  scale  with  g(p,\),  the  nondi¬ 
mensional  result  in  Eq.  (6)  is  independent  of  material  orthotropy. 

The  same  approach  is  used  to  address  the  stability  of  cracks  in 
finite  bodies.  The  effects  of  a/W  on  the  energy  release  rate  can  be 
described  generically  by 


G  =  H(a/W )  g(p,\) 


(7) 


where  H(a/W)  is  a  nondimensional  function  that  depends  on  both 
the  specimen  geometry  and  the  loading  configuration.  Combining 
Eqs.  (5)  and  (7).  the  critical  tearing  modulus  becomes 


T,W  W  (  I  \  dll(arJW) 
~G^~^o  +  d(a/W) 


In  the  limit  where  a</W  — »  0,  the  first  term  on  the  right  side 
dominates  and  the  critical  tearing  modulus  reduces  to  the  I  in  Eq. 
(6),  as  required. 

If  T  <  Tt.  the  instability  occurs  at  the  initiation  toughness  at  the 
Griffith  stress,  cr0,  written  in  nondimensional  form  as 

vlirW  1 

-1 -  _  -  (Q) 

EG0  (flo/VF)  H(aoIW)  g(p,\) 

Conversely,  if  T  >  Tc,  stable  crack  growth  occurs  until  the 
conditions  in  Eq.  (5)  are  satisfied.  This  scenario  is  represented  by 
the  tangent  construction  in  Fig.  1(b).  Provided  that  instability 
occurs  within  the  linear  region  of  the  resistance  curve,  the  critical 
values  of  crack  length,  af,  resistance,  G,,  and  stress,  a,,  are  found 
(from  Eq.  (5))  to  be 


(a) 


(b) 


(c) 


Fig.  1.  Tangent  construction  plots,  showing  the  following:  (a)  the 
transition  from  stable  to  unstable  crack  growth  for  =  0;  (b)  stable 
crack  growth  (starting  at  o0)  and  onset  of  instability  (at  o,)  in  the  linear 
portion  of  the  resistance  curve  for  finite  ti0 IW\  and  (c)  unstable  crack 
growth  at  the  steady-state  toughness.  (Here  the  material  is  assumed  to  be 
elastically  isotropic  such  that  ,t;(p.X)  =  I.) 


fl0  a ,  (W  /  I  \dW(fl,/VF)y  '  (  C0 
IF  ~  IF  U  +  \H(a,/W))  dUdW)  )  +  7TF 
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G r  WT (W  (  1  \ 

Go  _  ~Go  \«r  +  VHaJW)) 


dH(a/W) 
d(  a/W) 


aj-nW  TW  (  1  \ 

EG0  ~  ~G0  \(«/W)  H{am  g(p,X)j 
X /w  /  1  NdH(«,/W) 


\as  \H(ajW)J  d (a/W)  )  v  ' 

If  the  predicted  critical  resistance  G,  exceeds  Gss  (Fig.  1(c)),  the 
strength  is  obtained  at  G  =  Gss  at  a  total  crack  length. 


«»  =  (to  +  =  n0  +  - 

The  corresponding  stress  is 


G„  -  G« 


arir  W  f  r«o/C0\  G0/  CoMl 


III.  Results  Common  Test  Configurations 

Results  for  a  number  of  notched  configurations  commonly  used 
for  material  characterization  have  been  computed  using  the  anal¬ 
ysis  in  Section  II  along  with  established  solutions  for  H{a/W).  The 
configurations  of  interest  are  as  follows:  (i)  center-notched  tension 
(CNT),  (ii)  single-edge-notched  tension  (SENT),  (iii) double-edge- 
notched  tension  (DENT),  (iv)  single-edge-notched  (pure)  bending 
(SENB),  and  (v)  single-edge-notched  beam  in  three-point  bending 
(SENTPB)  with  a  span  to  width  ratio  of  S/W  =  4.  Pertinent 
solutions  for  H{a/W)  are  summarized  in  Table  1; 10  each  is  accurate 
within  0.5%  over  the  entire  range  of  a/W. 

Straightforward  analytical  results  for  the  critical  tearing  modu¬ 
lus  are  obtained  for  the  limiting  case  in  which  the  uncracked 
ligament  becomes  small;  i.e.,  a0/W  — >  1 .  For  SENT,  SENTPB.  and 
SENB,  the  limiting  solutions  for  H(aJW)  take  the  form10 

<l5> 

where  (3,  depends  on  the  loading  configuration  but  is  independent 
of  ag/W.  Combining  this  result  with  Eq.  (8)  yields 


1\W  2  +  W/a0 

~Go  =  I  ~  "a/W  "4) 

Similarly,  for  both  CNT  and  DENT,  the  limiting  solutions  for 
H(aJW)  take  the  form10 

HU‘JW)  -  (.  -  vw  1151 

where  p2  is  another  constant,  independent  of  aJW.  Combining 
with  Eq.  (8)  yields 

1\W  1 

”C7  =  («o/W)(l  -ag/W)  tl6) 

Figure  2  shows  the  variation  in  the  critical  tearing  modulus. 
T^.W/Gg,  with  notch  length  for  the  five  loading  configurations1  (Eq. 
(8)  and  Table  1)  as  well  as  for  the  limiting  cases  of  short  notches 
and  short  ligaments.  In  each  of  the  former  cases,  TVW/G„  exhibits 
a  minimum  at  a  notch  length  in  the  range  a0 /W  =  0.3-0.5  and 
becomes  large  in  the  limits  ag/W  — »  0  and  itg/W  —*  I.  Thus,  for  a 
prescribed  value  of  TW/G0.  crack  stability  is  only  obtained  over  an 
intermediate  range  of  aJW ;  both  short  notches  and  long  notches 
promote  unstable  growth.  Furthermore,  7/ W/G0  approaches  the 
short  notch  limit  (Eq.  (6))  when  aJW  <  0. 1  and  the  pertinent  short 
ligament  limit  (Eq.  (14)  or  ( 16))  when  aJW  >  0.9. 

The  figure  also  shows  the  sensitivity  of  7 \.W/G0  to  the  specimen 
type.  DENT  exhibits  the  lowest  value  over  the  entire  range  of  aJW 
and  hence  the  greatest  propensity  for  stable  crack  growth.  SENT  is 
at  the  other  extreme  with  the  highest  value  of  l\.W/G„.  by  as  much 
as  a  factor  of  2-3  relative  to  that  for  DENT.  CNT.  SENTPB.  and 
SENB  give  intermediate  results. 

The  effects  of  the  normalized  tearing  modulus,  TW/Ga,  on  crack 
stability  in  the  SENB  specimen  are  illustrated  in  Fig.  3.  For  the 
purpose  of  identifying  the  range  of  notch  lengths  over  which  crack 
stability  is  obtained  at  initiation,  the  critical  tearing  modulus  is 
replotted  in  Fig.  3(a)  and  compared  with  selected  values  of  TWIG0: 
10,  20,  and  30.  The  critical  points  arc  indicated  by  the  solid 
symbols  and  labeled  a-f.  The  amount  of  crack  growth  (Eq.  ( I  Orr ))- 


*  Results  lor  SKNTPB  and  SF.NB  are  almost  indistinguishable  from  one  another 
For  clarity,  the  former  results  are  not  included  in  the  figure. 


Table  I.  Summary  of  the  Functions  H(a/W)  for  Common  Loading  Configurations 

(from  Ref.  8) 


Loading  configuration 


IHifl/w)!1'-' 


[I  -  0.025(n/W)2  +  0.06 (a/W) 


I  l  TW  \ 

1  ^  ysec  \2vv7 


l2W  /tt<7  \ 

I  mt  ‘an  \2ivj 


0.752  +  2.02 Ui/W)  +  0.37  I  -  sin 


1,122  —  0.561  (a/W)  -  0.205 (a/W)2  +  0.47 1  (a/W}'  -  0. 1900//W)1 
,  1  -  a/W 


I2W  /it  n) 

I - tan  377/ 

I  it  a  \2W ) 


0.923  +  0.199  I  -  sin 


1.99  -  (a/W){  1  -  o/lV)(2.l5  -  3.93 (a/W)  +  2.1(alW)2) 
v%(l  +  2a/W){]  -  a/W)'12 


SENTPB  (S/W  =  4) 
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Fig.  2.  Effects  of  notch  length  and  specimen  type  on  the  critical  tearing 
modulus  above  which  stable  crack  growth  occurs. 


the  fracture  resistance  (Eq.  ( 10/>)),  and  the  stress  at  instability  (Eq. 
(10c))  are  plotted  in  Figs.  3(b-d),  respectively,  along  with  the 
critical  points,  a-f.  Outside  the  range  that  leads  to  stable  crack 


(a) 


growth  (e.g.,  to  the  left  of  b  or  to  the  right  of  e  for  TWI(J0  =  20). 
unstable  growth  occurs  at  G',/G0  =  1  at  the  Griffith  stress  (Eq.  (9). 
Fig.  3(c)).  In  the  intermediate  range  of  Oq/W,  some  stable  crack 
growth  is  obtained  (Fig.  3(b))  and  the  strength  is  elevated  relative 
to  the  Griffith  stress.  The  amount  of  strength  elevation  and  the 
range  of  notch  lengths  over  which  this  benefit  is  realized  increases 
with  increasing  TW/G0. 

It  is  noteworthy  that  when  viewed  over  an  intermediate  range  of 
afW,  the  shapes  of  the  strength  curves  (Fig.  3(d))  are  not 
strikingly  different  from  one  another.  Consequently,  experimental 
measurements  of  strength  for  several  notch  sizes  could  be  misin¬ 
terpreted  using  linear  elastic  fracture  mechanics  (LEFM),  assum¬ 
ing  a  uniquely  valued  toughness,  and  thus  the  effects  of  the 
resistance  curve  could  be  overlooked.  A  more  critical  distinction 
can  be  made  by  also  performing  tests  on  specimens  of  different 
widths  to  vary  TW/G0,  or  on  other  specimen  types,  especially  ones 
that  yield  substantially  different  values  of  TCW/G0. 

To  illustrate  the  effects  of  specimen  type,  results  are  presented 
in  Fig.  4(a)  for  the  two  extreme  cases:  SENT  and  DENT.  To 
distinguish  the  results  from  those  obtained  from  LEFM,  they  are 
replotted  in  Fig.  4(b)  in  terms  of  or,/<70.  The  latter  ratio  represents 
the  strength  elevation  above  the  Griffith  stress  due  to  the  resistance 
curve.  The  results  demonstrate  that  the  beneficial  effects  of  the 
rising  fracture  resistance  are  more  pronounced  in  DENT  than  in 
SENT.  For  instance,  for  TW/G0  =  20,  the  maximum  strength 
elevation  above  the  initiation  stress  is  only  about  18%  at  a^W  =» 


(b) 


<d) 


Fig.  3.  Effects  of  tearing  modulus  on  crack  stability  and  strength  for  SENB. 
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(a) 


Fig.  4.  Effects  of  specimen  type  on  strength.  Note  that  the  benefits  of  the 
resistance  curve  are  significantly  greater  for  DENT  than  for  SENT  in  the 
intermediate  range  of  a<JW,  consistent  with  the  differences  in  the  critical 
tearing  moduli  shown  in  Fig.  2. 


0.3;  the  corresponding  value  for  DENT  at  the  same  notch  length  is 
70%.  It  is  expected  that  performing  comparative  experiments  such 
as  these  would  allow  a  critical  evaluation  of  the  contributions  to 
strength  from  the  resistance  curve.  Indeed,  results  of  this  type  are 
used  in  designing  some  of  the  experiments  described  below. 

The  preceding  results  are  based  on  the  assumption  that  the 
resistance  at  fracture  is  below  the  steady-state  toughness  and  hence 
the  strength  is  given  by  Eq.  (10c).  Otherwise,  the  strength  is  given 
by  Eq.  (12).  One  set  of  illustrative  results  for  the  latter  case,  where 
Gt  >  Gss,  is  shown  in  Fig.  5.  They  correspond  to  the  SENB 
geometry  for  several  values  of  GJG0  at  a  fixed  value  of  TW/G0. 
The  results  from  Eqs.  (9)  and  (10c)  bound  those  from  Eq.  (12) 
when  GJG0  falls  below  the  maximum  value  of  G,/G0  (Fig.  3(c)). 


IV.  Experimental  Procedure 

The  CFCCs  used  for  experimental  measurements  comprise  seven 
layers  of  3000  denier  Nextel  720  fiber  cloth  in  an  8-harness  satin 
weave  and  a  porous  mullite-alumina  matrix.  The  matrix  was  intro¬ 
duced  in  two  steps:  (i)  pressure-assisted  infiltration  of  a  mullite- 
alumina  slurry  through  the  fiber  preform  and  (ii)  impregnation  and 


*S> 

05 


JZ 

O) 

c: 

£ 

55 


Fig.  5.  Effects  of  steady-state  toughness  on  notched  strength,  when  (7,,  < 
Gr.  When  Gss  >  G,  over  the  entire  range  of  aJW.  the  instability  occurs  in 
the  rising  part  of  the  resistance  curve  and  hence  the  strenelh  is  independent 
of  G.s. 


pyrolysis  of  an  AI2CI(OH),  solution  (an  alumina  precursor).  Process¬ 
ing  details  were  essentially  identical  to  those  described  earlier"  with 
two  exceptions.  First,  the  concentration  of  the  precursor  solution  was 
increased  to  give  a  volumetric  yield  of  6.4%,  compared  with  3%  in 
previous  studies.  The  more  concentrated  solution  allows  a  prescribed 
matrix  density  to  be  achieved  with  fewer  impregnation  and  pyrolysis 
cycles.  In  the  present  study,  four  such  cycles  were  used,  yielding  a 
final  matrix  porosity  of  32%.  Second,  to  prevent  segregation  of  the 
precursor  during  drying,  the  panels  were  gelled  immediately  after 
precursor  impregnation  by  immersing  them  in  an  NH40H  environ¬ 
ment  for  6  h.  The  absence  of  segregation  was  confirmed  through  a 
series  of  matrix  hardness  measurements  across  the  plate  thickness. 
The  composite  panel  thickness  was  3.3  mm. 

As  a  parenthetical  note,  the  matrix  in  the  present  composite  had 
been  designed  to  be  considerably  harder  than  that  in  composites 
described  in  earlier  studies. 121:1  This  selection  was  motivated  by 
the  attendant  improvements  in  the  off-axis  properties,  e.g..  0'790" 
shear  strength.  Although  not  reported  here,  significant  elevations 
in  shear  properties  are  indeed  obtained  with  increasing  matrix 
hardness.  A  comparative  study  on  the  effects  of  matrix  hardness  on 
the  composite  performance  in  both  matrix-  and  fiber-dominated 
loadings  will  be  reported  elsewhere. 

In-plane  elastic  properties  were  obtained  from  unnotched  ten¬ 
sion  and  losipescu  shear  tests,  performed  parallel  to  the  fiber 
directions.  The  tensile  specimens  were  of  a  dog-bone  configuration 
with  a  gauge  section  8  mm  wide  and  60  mm  long.  The  losipescu 
specimens  were  19  mm  wide  and  contained  a  pair  of  V-notchcs 
with  an  included  angle  of  1 05",  to  a  depth  of  3.8  mm.  Both 
specimen  types  were  instrumented  with  strain  gauges.  The  key 
results  from  these  tests  are  as  follows:  E,  =  64  GPa.  v,2  =  0.13. 
and  p, 1 2  =  14  GPa.  Combining  with  Eqs.  (3)  and  (4)  yields  p  = 
2.15  and  g(p)  =  1.51. 

Additional  mechanical  tests  were  performed  in  the  0790'' 
orientation  on  three  types  of  notched  specimens:  SENTPB,  SENT, 
and  DENT.  The  geometric  parameters  were  selected  to  probe  the 
effects  of  W  at  fixed  aJW  and  the  effects  of  af\V  at  fixed  W.  The 
pertinent  values  are  listed  in  Table  II.  For  SENTPB,  the  span 
length  was  SIW  =  4.  The  length  of  the  SENT  specimens  was 
HIW  =  1.0:  for  the  DENT  specimens.  H/W  a  3.  Crack-growth 
measurements  were  made  using  crack  gauges  affixed  to  the 
specimen  surface  (Krak-Gauge,  model  A05-CE,  Hartrun  Coqt.,  St. 
Augustine,  FL).  These  measurements  were  checked  against  direct 
optical  measurements  on  the  opposing  surface.  Image  contrast  was 
enhanced  by  coating  the  viewing  surfaces  with  Au-Pd.  The  cracks 
were  observed  to  be  essentially  planar,  with  no  evidence  of 
inelastic  deformation  or  damage  away  from  the  crack  tip.  The 
nominal  fracture  resistance  was  calculated  via  Eq.  (7),  using  the 
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Table  II.  Notched  Loading  Configurations  Used  in  the 
Experiments 


a,JW 

W  (mni) 

S/W 

w\v 

SENTPB 

0.5 

6,  12,24 

4 

_ 

SENT 

0.5 

6,  12.  24 

- 

1.0 

0.35,  0.5,  0.65 

24 

- 

1.0 

DENT3 

0.5 

6.  12 

- 

>3 

following  the  usual  convention  for  DENT  specimens,  the  full  specimen  width  is 
denoted  2 W.  its  length  is  2 H,  and  the  length  of  each  of  the  two  notches  is  nti. 


measured  load  and  crack  length,  coupled  with  the  computed  value 
of  H(a/W)  and  g(p)  =  1.51.  For  SENT,  the  effects  of  specimen 
rotation  in  the  gripped  sections  were  incorporated  using  the 
approach  described  in  the  Appendix.  The  SENTPB  specimens 
were  also  instrumented  with  a  spring-loaded  one-armed  extensom- 
eter  to  monitor  load-line  displacement.  The  work  of  fracture, 
obtained  from  the  area  under  the  load-displacement  curve,  was 
used  as  a  measure  of  the  steady-state  toughness.14  Base-line  values 
of  unnotched  strength  were  obtained  from  uniaxial  tension  tests 
(described  above)  and  four-point  bend  tests.  The  bend  specimens 
were  either  6  or  12  mm  wide  and  loaded  with  a  span  S/W  «*  9. 

The  experimental  results  are  summarized  in  Figs.  6  and  7.  In  all 
cases,  stable  crack  growth  had  been  observed  before  the  stress 
maximum  and  the  resistance  curves  were  essentially  linear.  All  the 
tensile  geometries  yielded  common  resistance  parameters:*  nota¬ 
bly,  G0  *»  900  ±  200  J/m2  and  T  =  400  ±  100  kJ/m3.  In  contrast, 
the  SENTPB  specimens  exhibited  significantly  higher  fracture 
resistance,  characterized  by  G0  ~  1500  ±  150  J/m2  and  a  tearing 
modulus  T  =  950  ±  100  kJ/m  .  It  is  surmised  that  the  differences 
are  associated  with  the  7-stresses  and  their  effect  on  inelasticity  in 
the  notch  tip  region.  That  is,  for  SENTPB,  the  T-stress  is  positive 
over  the  range  of  af\V  used  in  the  present  experiments.9  When 
positive,  the  T-stress  promotes  matrix  microcracking  ahead  of  the 
notch  tip  which,  in  turn,  reduces  the  normal  stress  along  the 
incipient  fracture  plane  and  hence  increases  fracture  resistance.13 
In  contrast,  for  the  tensile  geometries,  the  T-stress  is  negative  over 
the  pertinent  range  of  aJW,  suppressing  the  inelasticty  that  leads 
to  increased  resistance. 

The  one  notable  exception  to  the  behavior  described  above 
pertains  to  the  narrowest  SENTPB  specimens.  In  this  case,  the 
fracture  resistance  was  anomalously  low  (relative  to  that 
measured  on  the  wider  specimens)  and  the  notched  strength  was 
essentially  equal  to  the  unnotched  strength.  Evidently,  fracture 
is  controlled  by  the  net-section  stress  rather  than  the  A-field 
associated  with  the  notch.  Consequently,  the  analysis  presented 
here  does  not  apply. 

V.  Comparisons  between  Experiment  and  Theory 

Notched  strength  predictions  were  made  using  the  analysis 
presented  in  Section  II  and  the  range  of  resistance  parameters 
shown  in  Fig.  7:  G0  =  1500  ±  150  J/m2  and  T  =  950  kJ/m3  for 
flexure  and  G0  =  900  ±  200  J/m2  and  T  =  400  kJ/m3  for 
tension.  The  results  are  plotted  in  Fig.  6.  Also  shown  for 
comparison  are  the  values  obtained  from  LEFM:  notably,  the 
Griffith  stress  (Eq.  (9)),  calculated  using  the  initiation  tough¬ 
ness.  The  transition  between  the  two  corresponds  to  T,.W/G0. 
Broadly,  the  predictions  correlate  well  with  the  experimental 
measurements 

Viewed  collectively,  the  results  provide  insights  into  the  tran¬ 
sitions  in  behavior  with  increasing  notch  length.  For  instance,  for 
the  SENT  specimens  with  a<JW  =  0.5,  three  domains  exist:  (i) 
net-section-stress-dominated,  denoted  1,  for  a0  s  I  mm,  (ii) 


*The  tearing  modulus  was  obtained  through  a  linear  regression  analysis  of  the 
resistance  curve  for  crack  growth  up  to  An/lV  0.25.  The  curves  exhibited 
nonlinearity  beyond  this  point,  a  manifestation  of  large-scale  bridging.  In  all  pertinent 
cases,  the  peak  stress  was  attained  when  the  crack  was  within  the  linear  portion  of  the 
resistance  curve. 


(a) 


(b) 


Notch  length,  a0/W 

(c) 

Fig.  6.  Notched  strength  of  oxide  composite:  fa)  SENTPB.  ufW  =  0.5: 
(b)  DENT  and  SENT.  aJW  =  0.5.  and  (c)  SENT.  W  =  24  mm.  Solid 
symbols  represent  individual  experimental  measurements,  solid  lines  are 
values  predicted  from  the  model,  and  dashed  lines  are  from  LEFM. 
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Crack  extension,  Aa  (mm) 


Fig.  7.  Typical  fracture  resistance  curves  obtained  from  SENTPB  and 
SENT  tests.  Lines  represent  envelope  of  fit  values,  with  G0  =  1500  ±  150 
i/m2  and  T  —  950  kj/m’  for  flexure  and  G0  =  900  ±  200  and  T  =  400  for 
tension.  The  steady-state  toughness  obtained  from  SENTPB  tests  is  Gss  = 
5300  ±  400  J/m2,  well  beyond  the  values  relevant  to  the  regime  where 
instability  occurs. 


initiation-dominated,  denoted  II,  for  1  mm  <  a0  <  5  mm,  and  (iii) 
T-dominated,  denoted  Ill,  for  a0  s  5  mm.  For  the  SENTPB 
specimens,  the  transition  from  1  to  II  occurs  at  stresses  that  are 
essentially  the  same  as  the  unnotched  strength.  The  implication  is 
that  the  behavior  transitions  directly  from  I  to  III  without  passing 
through  II. 

To  further  explore  these  transitions,  calculations  have  been  per¬ 
formed  for  other  values  of  aJW.  Representative  results  for  the 
SENTPB  geometry  with  aJW  in  the  range  0.25-0.75  are  plotted  in 
Fig.  8;  similar  trends  are  obtained  for  the  other  geometries.  For  clarity, 
results  are  presented  only  for  the  average  values  of  the  initiation 
toughness  (G0  =  1500  J/m2),  the  tearing  modulus  (T  =  950  kJ/m3), 
and  the  steady-state  toughness  (Gss  =  5300  J/m2).  Furthermore,  they 
are  presented  over  a  broad  range  of  notch  lengths  to  capture  the  full 
range  of  behavior.  Two  new  features  are  notable.  First,  there  exists  an 


Fig.  8.  Effects  of  afW  on  notched  strength  of  SENTPB  specimen.  The 
open  circles  show  the  behavioral  transitions.  For  the  transition  II  — >  III  (top 
row  of  circles),  the  critical  notch  length  is  given  by  Eq.  (8);  for  the 
transition  III  — *  IV  (lower  row  of  circles),  it  is  given  by  Eq.  (11). 


additional  domain,  denoted  IV,  for  the  longest  notches,  where  the 
strength  is  dictated  by  the  steady-state  toughness.  In  this  domain,  the 
notched  suength  asymptotically  approaches  an  inverse  square  rtxii 
dependence  on  notch  length,  analogous  to  that  of  the  initiation- 
dominated  domain.  Second,  the  notch  lengths  corresponding  to  the 
transitions  increase  substantially  with  increasing  Uq/W.  For  instance, 
for  a  very  deep  notch,  characterized  by  aJW  =  0.75,  the  transition 
from  I  to  II  occurs  for  a0  **=  14  mm.  in  contrast,  with  notches  of  this 
length  in  a  specimen  with  aJW  =  0.25,  the  behavior  is  well  within 
domain  IV.  As  a  consequence,  striking  differences  in  behavior 
would  be  obtained  over  the  practical  range  of  notch  lengths, 
from  2  to  14  mm.  In  the  former  case  (<70/W  =  0.75),  the  strength 
would  be  indistinguishable  from  the  unnotehed  strength;  in  the 
latter  (n0/W  =  0.25),  significant  notch  sensitivity  would  be 
apparent.  The  latter  tests  are  recommended  to  most  critically 
probe  notch  effects. 


VI.  Concluding  Remarks 

A  theoretical  framework  for  predicting  the  notch  sensitivity  of 
strength  of  ceramic  composites  with  rising  fracture  resistance  has 
been  presented.  An  attractive  feature  of  the  resulting  solutions  is 
their  relatively  simple  analytical  form:  a  consequence  of  the 
assumed  linear  form  of  the  resistance  curve.  The  framework 
provides  a  consistent  description  of  the  strength  of  the  oxide 
composite  studied  here  when  the  notch  length  is  large  and  the 
corresponding  strength  is  low  in  relation  to  the  unnotched  strength. 
Parametric  studies  illustrate  the  salient  features  in  the  transitions  in 
fracture  behavior  and  their  dependence  on  both  the  notch  length 
and  its  normalized  value  laJW). 

In  addition  to  its  use  in  rationalizing  experimenial  results,  the 
theoretical  framework  should  find  utility  in  designing  critical 
experiments  to  probe  notch  sensitivity  of  other  composite  materi¬ 
als.  For  instance,  with  preliminary  estimates  of  G0  and  7',  calcu¬ 
lations  could  be  performed  to  identify  the  specimen  sizes  and  types 
for  which  stable  crack  growth  should  occur  and  experiments 
performed  accordingly.  This  combined  theoretical/experimental 
approach  is  expected  to  yield  significantly  improved  understand¬ 
ing  of  notch  sensitivity  relative  to  that  obtained  through  an  ad  hoc 
experimental  program. 

On  a  cautionary  note,  the  framework  is  not  expected  to  apply  to 
composites  that  exhibit  large-scale  inelasticity  in  the  notch-tip 
region.  In  composites  with  weak  (porous)  matrixes,  such  inelas¬ 
ticity  can  occur  through  a  shear  mode  parallel  to  the  loading 
direction.15  Alternatively,  in  systems  with  strong  (dense)  matrixes, 
it  can  occur  in  a  tensile  mode  through  multiple  matrix  crack¬ 
ing.1516  When  the  spatial  extent  of  the  inelasticity  becomes 
comparable  to  the  structural  dimensions,  the  inelasticity  must  to  be 
taken  into  account  explicitly1316  and  not  simply  embodied  in  the 
initiation  toughness. 

Finally,  the  degree  of  notch  sensitivity  displayed  by  the  present 
oxide  composite  is  significantly  greater  than  that  in  systems  with 
weaker  matrixes.12 13  However,  the  lower  notch  sensitivity  asso¬ 
ciated  with  weaker  matrixes  is  accompanied  by  reductions  in 
off-axis  properties,  notably  the  0°/90f'  shear  strength.  Trade-offs 
between  fiber-dominated  and  matrix-dominated  properties  are  the 
subject  of  a  current  investigation. 


Appendix 

Experience  with  the  SENT  geometry  indicates  that  small 
amounts  of  specimen  rotation  can  occur  under  nominally  fixed 
grip  conditions,  thereby  altering  the  associated  stress  intensity 
factor.  The  magnitude  of  this  effect  is  determined  from  measure¬ 
ments  of  the  crack-mouth-opening  displacement  (CMOD)  coupled 
with  the  following  analysis.17 

The  stresses  at  the  specimen  ends  can  be  partitioned  into  two 
components:  (i)  a  uniform  tensile  stress,  ox,  equal  to  the  average 
measured  stress  and  (ii)  a  bending  stress,  oB.  proportional  to  crT. 
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Fig.  9.  Effects  of  specimen  rotation  on  the  normalized  bending  stress. 
Solid  symbols  represent  individual  experimental  measurements,  solid  lines 
are  values  predicted  from  the  two  limiting  end  conditions,  and  the  dashed 
line  is  an  empirical  fit  to  the  experimental  data. 


but  of  unknown  magnitude.  The  stress  intensity  factor  associated 
with  these  two  stress  components  is 

K  =  Fr(a/W)  Y(p)crx  fut  -  Fb(ciIW)  K(p)crB  nMyi 

=  Ffa/W)  n(.K  frn  (A-l) 

where 

T(p)  =  1  +  0.1(p  -  1)  -  0.01 6(p  -  l)3  +  0.002(p  -  1)J 

(A-2) 

Here.  FT(a/W)  and  FB(a/W)  are  the  nondimensional  functions  for 
SENT  and  SENB  given  in  handbooks,10  FB(a/W)  is  the  function 
that  characterizes  the  stress  intensity  factor  in  the  present  test 
configuration,  and  K(p)  accounts  for  material  orthotropy.7'8  The 
corresponding  CMOD  is 

8  =  W  [  V-{alW)^f)  ~  (A-3) 

where  VT(alW)  and  VB(a/W)  are  the  nondimensional  functions  for 
CMOD  obtained  from  handbooks,  and  g(p,\)  characterizes  the 
material  orthotropy  as  described  in  Section  II.  Rearranging  Eq. 
A-3  yields  the  stress  ratio: 


crB  W  Hp) 

o>  Vb  Vggfp.X) 


)  /  8  E 

X)  \4ocrT 


Combining  Eqs.  A-l  and  A-4  gives  the  normalized  stress  intensity 
factor  as 


crT  \ftta 


/  crT\ 


Fb  K(p)  /  5  E  \1 

=FT-^h+irate)J  (A-5) 

The  stress  intensity  factor  had  been  obtained  from  Eq.  A-5  coupled 
with  measurements  of  8  in  the  elastic  regime  and  the  known  values 
of  the  other  parameters. 

The  degree  to  which  the  bending  stress  had  been  reduced 
relative  to  the  uniform  displacement  condition  was  ascertained 
from  a  comparison  of  the  experimentally  determined  stress  ratio  in 


Eq.  A-4  with  the  predicted  value.  The  results  are  plotted  in  Figure 
9.  Evidently  the  experimental  results  lie  close  to,  but  slightly 
below,  those  for  the  uniform  displacement  condition.18  The  dashed 
line  in  the  figure  shows  an  empirical  fit  through  the  data,  given  by 
the  product  of  the  uniform  displacement  solution  and  m{a/w)  = 
1  -  0.5(fl/n')(l  -  <?/ IT).  This  fit  accurately  represents  the  experi¬ 
mental  data  and  provides  the  correct  limits  to  the  stress  ratio  in  the 
limits  of  a/W  —>  0  and  a/W  —*  I . 


Nomenclature 

Crack  length 

Initial  crack  or  notch  length 
Total  crack  length  at  unstable  crack  growth 
Total  crack  length  at  onset  of  steady  state 
Characteristic  bridging  length  scale 
Young’s  modulus 
Plane  strain  Young’s  modulus 
Orthotropy  function  (Eqs.  I  and  4) 

Energy  release  rale 

Fracture  resistance  at  unstable  crack  growth 
Initiation  toughness 
Fracture  resistance 
Steady-slate  toughness 
Nondimensional  function  (Table  1) 

Loading  span 
Tearing  modulus 
Critical  tearing  modulus 
Critical  crack  opening  displacement 
Specimen  width 
Orthotropy  function  (Eq.  A-2) 

Crack  growth 

Crack  growth  at  onset  of  steady  state 
Elastic  orthotropy  parameter  (Eq.  2) 

Shear  modulus 
Poisson’s  ratio 

Elastic  orthotropy  parameter  (Eq.  3) 

Stress 

Stress  at  onset  of  unstable  crack  growth 
Stress  at  onset  of  crack  growth 
Unnotched  strength 
Characteristic  bridging  stress 
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Controlling  Mechanical  Properties  of  Porous  Mullite/Alumina  Mixtures 
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This  study  focuses  broadly  on  synthesis  and  characterization  of 
porous  mullite/alumina  mixtures  for  use  as  matrices  in  oxide 
fiber  composites.  Specifically,  an  assessment  is  made  of  the  util¬ 
ity  of  a  precursor-derived  alumina  (PDA)  in  controlling  both  the 
modulus  and  the  toughness  of  mullite-rich  particle  mixtures. 
Property  changes  are  probed  through  models  of  the  mechanical 
behavior  of  bonded  particle  aggregates.  Consideration  of  the 
conditions  needed  to  cause  crack  deflection  at  a  fiber-matrix 
interface  yields  an  upper  allowable  limit  on  the  concentration  of 
PDA  to  ensure  damage  tolerance  in  a  fiber  composite.  The  pre¬ 
dicted  critical  concentration  lies  in  the  range  of  about  7-9%, 
depending  on  the  mullite/alumina  ratio  in  the  particle  slurry  and 
the  subsequent  aging  treatment.  Values  slightly  below  this  limit 
should  yield  composites  that  exhibit  a  desirable  balance  between 
fiber-  and  matrix-dominated  properties. 


I.  Introduction 

Damage  tolerance  can  be  enabled  in  oxide  ceramic  compos¬ 
ites  through  the  use  of  a  controlled  amount  of  fine-scale 
matrix  porosity,  obviating  the  need  for  a  fiber  coating.1^*  For 
successful  implementation,  the  matrix  must  satisfy  three  criteria: 
(i)  pore  retention  at  the  targeted  use  temperature,  (ii)  phase 
compatibility  with  the  fibers  (because  of  the  intimate  fiber-ma¬ 
trix  contact  in  the  absence  of  a  coating),  and  (iii)  adequate  me¬ 
chanical  integrity  to  ensure  acceptable  off-axis  composite 
properties.  When  the  reinforcements  are  either  alumina  or  mul¬ 
lite/alumina  fibers,  the  requisite  matrix  characteristics  can  be 
achieved  through  the  use  of  mullite/alumina  mixtures.  Typically, 
the  mullite  is  in  the  form  of  particles,  introduced  by  infiltration 
of  a  slurry  into  the  fiber  preform.  Its  main  role  is  to  form  a 
contiguous  nondensifying  network.  In  contrast,  the  alumina 
serves  as  a  binder  to  impart  mechanical  strength.  It  can  be  in¬ 
troduced  in  one  of  two  ways:  by  mixing  alumina  particles  into 
the  slurry,  or  by  subsequent  impregnation  and  pyrolysis  of  a 
liquid  alumina  precursor.5  These  lead  to  two  distinct  matrix 
topologies,  illustrated  in  Fig.  1.  The  primary  objective  of  the 
present  paper  is  to  compare  the  efficacy  of  the  two  approaches  in 
controlling  matrix  properties.  This  is  achieved  largely  through  a 
series  of  experiments  on  porous  mullite/alumina  compacts  with 
varying  amounts  of  precursor-derived  alumina  (PDA)  and  com¬ 
parisons  with  previously  reported  measurements  on  mullite/alu¬ 
mina  particle  mixtures  (without  precursor  additions).6  A  second 
objective  is  to  develop  and  assess  models  of  bonded  particle  ag- 
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gregates  for  use  in  describing  the  mechanical  properties  of  com¬ 
pacts  containing  a  precursor-derived  phase. 

The  two  routes  for  introducing  alumina  have  three  attributes: 
(i)  The  use  of  a  mixed  mullite/alumina  particle  slurry  allows 
both  matrix  phases  to  be  introduced  simultaneously.  By  con¬ 
trast,  the  precursor  route  requires  additional  steps,  beyond  that 
of  slurry  infiltration.  From  a  cost  viewpoint,  the  former  method 
is  preferred,  (ii)  The  presence  of  alumina  particles  in  the  slurry 
can  compromise  the  stability  of  the  mullite  network,  especially  if 
its  proportion  exceeds  the  percolation  threshold.6  Conversely,  if 
the  slurry  is  comprised  of  only  mullite  and  the  alumina  is  intro¬ 
duced  subsequently  via  the  precursor  route,  the  contiguity  of  the 
mullite  network  is  ensured,  (iii)  The  slurry  processing  route  re¬ 
sults  in  an  essentially  fixed  (high)  porosity.  Consequently,  a 
modest  upper  limit  on  the  strength  of  the  network  is  expected.  In 
contrast,  the  precursor  route  allows  for  filling  of  the  void  space 
between  the  particles  in  the  network,  with  attendant  elevations 
in  both  the  relative  density  and  the  mechanical  properties.  A 
further  benefit  of  the  precursor  route  might  be  obtained,  if, 
because  of  capillarity,  the  PDA  forms  preferentially  at  the 
junctions  of  the  mullite  particles.  The  present  experiments 
were  designed  to  address  some  of  these  issues. 

To  motivate  the  mechanical  measurements,  a  crack  deflection 
parameter,  Z,  is  used,  defined  by6'7 


1  =  0.134 


(£)('+:£) 


0.9 


0) 


where  T  is  the  toughness,  E  is  Young’s  modulus,  and  subscripts 
m  and  f  refer  to  matrix  and  fiber,  respectively.  A  matrix  crack  is 
predicted  to  arrest  or  deflect  at  the  fiber-matrix  interface  when  Z 
exceeds  a  critical  value,  defined  by  to  =  ryrm,  where  T;  is  the 
interface  toughness.  Since  the  nature  of  bonding  at  the  fiber- 
matrix  interface  is  similar  to  that  between  particles  in  the  matrix, 
the  interface  toughness  is  expected  to  be  comparable  with  that  of 
the  matrix,  i.e.  co  «  i .  Thus,  upon  setting  Z  =  to  =  1 ,  the  critical 
combination  of  rm/rr  and  Em/Er  is  obtained.  In  light  of  this 
prediction,  the  subsequent  experimental  study  focuses  on  the 
matrix  properties  rnl  and  Em.  An  assessment  is  then  made  of 
whether  the  crack  deflection  condition  is  satisfied  when  a  can¬ 
didate  matrix  is  combined  with  the  targeted  fiber,  namely  Nextel 
720. 

Despite  an  extensive  literature  on  the  mechanical  properties 
of  porous  ceramics,  only  a  small  fraction  of  the  effort  has  been 
devoted  to  the  development  and  assessment  of  mechanistic  mod¬ 
els  based  on  realistic  representations  of  pore  topology.  Perhaps 
the  greatest  deficiency  lies  in  the  low  porosity  domain,  wherein 
the  materials  consist  of  particle  aggregates  bonded  to  one  an¬ 
other  at  discrete  contact  points.  In  this  case,  macroscopic  prop¬ 
erties  are  extremely  sensitive  to  changes  in  contact  size,  even 
when  unaccompanied  by  changes  in  macroscopic  porosity.6'*'9 
Elastic  properties  of  such  aggregates  have  been  successfully 
modeled  by  treating  the  particle  junctions  as  Hertzian  contacts 
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Nonsinterabte  Sinterabte  phase: 

phase:  Mulile  Alumina  particles 


(a) 


Nonsinter  able  Sinterabte  phase: 

phase:  Muttte  Precursor-derrved  alumina 


Fig.  1.  Schematics  showing  the  two  matrix  topologies  produced  by  (a) 
mullite/alumina  particle  mixtures  and  (b)  mullite  particles  bonded  by 
precursor-derived  alumina. 


and  analyzing  the  aggregate  response  accordingly.10-13  A  similar 
approach  has  been  used  in  conjunction  with  the  discrete  element 
method  (DEM)  for  simulating  the  fracture  process.6  Otherwise, 
most  previous  studies  have  relied  on  purely  empirical  relations 
to  rationalize  experimental  measurements,  usually  couched  sole¬ 
ly  in  terms  of  the  porosity  level.8,14-16  Others  have  implemented 
models  based  on  rather  simplistic  (arguably  unrealistic)  micro- 
structural  representations,  e.g.  a  collection  of  aligned  cylinders 
of  varying  diameters,  each  represented  mechanically  as  a  one¬ 
dimensional  spring  with  stiffness  dictated  by  cylinder  diame¬ 
ter.17,18  Although  models  of  this  type  can  be  readily  analyzed, 
they  lack  the  realism  needed  for  establishing  connections  be¬ 
tween  microstructure  and  properties.  Analogous  phenomeno¬ 
logical  approaches19  have  proved  to  be  of  similarly  limited 
utility  in  modeling  fracture  energy. 

The  present  article  is  organized  as  follows.  Section  II  de¬ 
scribes  the  nature  of  the  materials  systems,  the  processing  route, 
and  the  measurement  procedures.  The  resulting  physical  prop¬ 
erties  and  microstructures  are  presented  in  Section  III.  The  key 
mechanical  properties — modulus  and  toughness — are  addressed 
in  Sections  IV  and  V.  Section  VI  considers  the  implications  on 
crack  deflection  in  fiber  composites.  For  this  purpose,  the  cal¬ 
ibrated  models  are  used  to  establish  the  critical  concentration  of 
the  precursor-derived  phase  at  which  the  crack  deflection  con¬ 
dition  is  no  longer  satisfied.  Conclusions  are  summarized  in 
Section  VII.  A  comprehensive  list  of  key  symbols  and  their  def¬ 
initions  is  presented  in  Table  I. 


II.  Materials  and  Test  Procedures 

Porous  mullite  compacts,  10  mm  x  50  mm  x  90  mm  in  size,  were 
fabricated  by  vacuum  filtration,  following  procedures  detailed 
elsewhere.5  The  particles  were  MU-107  mullite  (Showa  Denko 
K..K.,  Tokyo,  Japan),  with  an  average  diameter  of  1  pm,  dis¬ 
persed  in  aqueous  HNOj  at  a  pH  of  3.  Following  filtration,  the 
compacts  were  dried,  sintered  at  900°C  for  2  h  to  promote  par¬ 
tial  sintering  at  the  particle  junctions,  and  machined  into  me¬ 
chanical  test  specimens.  Some  specimens  were  subsequently 
impregnated  with  an  alumina  precursor  solution,  following  the 
procedures  described  below,  and  pyrolyzed  at  900°C  for  2  h.  All 
were  subjected  to  a  final  aging  treatment  of  either  2  h  at  1200°C 
(representing  the  as-processed  condition)  or  1000  h  at  1200°C 
(to  assess  thermal  stability). 

The  precursor  solution  was  made  by  slowly  digesting  Al  pow¬ 
der  in  a  solution  of  AlClj  •  6H20  in  deionized  water  at  85°C.  The 
fractional  volumetric  yield  of  alumina  was  determined  by  firing 
known  quantities  of  the  precursor  solution  at  either  900°  or 
1200°C  and  subsequently  measuring  both  the  mass  and  the  den¬ 
sity  of  the  pyrolyzed  product.  The  corresponding  densities,  3.48 
and  3.99  g/cm,3  are  consistent  with  y-  and  a-alumina,  respec¬ 
tively.  The  subsequent  yield  measurements  are  reported  on  the 
basis  of  the  y-alumina;  the  two  yields  are  related  by  yjy1  =  py/ 
pa  ss  0.87  where  y  is  the  volumetric  yield,  p  is  the  mass  density, 
and  the  subscripts  refer  to  the  two  phases. 


Table  I.  Symbols  and  Definitions 

a  Junction  radius 

oq  Initial  junction  radius  (before  precursor  impregnation) 
A  Surface  area  per  particle 

Ao  Initial  surface  area  per  particle  (before  precursor 
impregnation) 

E(  Fiber  Young’s  modulus 

Em  Matrix  Young’s  modulus 

Es  Young’s  modulus  of  solid  phase 

h0  Height  of  truncated  spherical  cap  (Fig.  Al) 
n  Sintering  exponent 

p0  Initial  compact  porosity  (before  precursor 
impregnation) 
p  Final  porosity 

t  Time 

/R  Reference  time 

R  Particle  radius 

Em  Effective  particle  radius 

Kpda  Volumetric  concentration  of  precursor-derived 
alumina 

Ko  Particle  volume 

ya  Fractional  volumetric  yield  (ot-phase) 

_yY  Fractional  volumetric  yield  (y- phase) 

a  Normalize  junction  size  (a//?err) 

ao  Normalize  initial  junction  size  ( a0/R ) 
p  Nondimensional  precursor  “segregation”  parameter 

Tf  Fiber  toughness 

rm  Matrix  toughness 

r,  Interface  toughness 

Tj  Junction  toughness 

p  Mass  density 

E  Crack  deflection  parameter  (Eq.  (1)) 

<D  Toughness  ratio  (r i/fj) 


The  precursor  concentration  was  selected  in  the  following 
way.  To  achieve  large  porosity  reductions  with  minimal  impreg¬ 
nation  cycles,  concentrated  solutions  are  preferred.  However,  a 
practical  upper  limit  on  the  allowable  solution  concentration  is 
obtained  because  of  solution  viscosity  and  compact  permeabil¬ 
ity.  That  is,  solution  impregnation  into  the  compact  at  atmos¬ 
pheric  pressure  becomes  impractical  when  the  concentration, 
and  hence  viscosity,  are  excessively  high.  An  experimental  de¬ 
termination  of  this  limit  was  made  by  impregnating  samples  of 
known  porosity  with  solutions  of  varying  concentrations,  firing 
at  900°C,  and  measuring  the  weight  gain  because  of  the  addition 
of  the  PDA.  The  results  are  plotted  in  Fig.  2.  Initially,  the  weight 
increases  linearly  with  solution  concentration.  The  rate  is  con¬ 
sistent  with  a  prediction  based  on  the  measured  precursor  yield, 
assuming  that  all  of  the  pore  space  within  the  compact  is  filled 
by  the  solution.  Beyond  a  critical  concentration  (  «  340  mg 
AbOj/cm3  of  solution)  corresponding  to  a  precursor  yield  of 
about  9.5%,  the  weight  gain  diminishes  rapidly,  implying  min¬ 
imal  impregnation.  Based  in  part  on  these  results,  two  solution 
concentrations  were  selected  for  subsequent  experiments: 
yy  se  3%  and  8%.  The  latter  value  lies  somewhat  below  the 
maximum  dictated  by  permeability.  The  former  value  was  se¬ 
lected  to  be  consistent  with  that  used  in  previous  studies  on  po¬ 
rous  matrix  composites.3,4  The  number  of  impregnation  cycles 
was  varied  from  2  to  10.  As  shown  later,  the  properties  were 
found  to  depend  only  on  the  total  amount  of  PDA,  independent 
of  the  specific  combination  of  precursor  concentration  and 
number  of  impregnation  cycles. 

Previous  work  on  composite  fabrication  had  revealed  some 
benefits  in  using  a  hybrid  approach  involving  both  alumina  par¬ 
ticle  additions  to  the  slurry  and  subsequent  precursor  impreg¬ 
nation  and  pyrolysis.1,2  The  alumina  particles  were  found  to 
improve  the  mechanical  integrity  of  the  green  panels  following 
the  initial  low-temperature  heat  treatment,  relative  to  that  con¬ 
taining  pure  mullite,  and  allowed  handling  during  subsequent 
precursor  impregnation  without  imparting  damage  to  the  panel. 
Furthermore,  density  measurements  on  various  mullite/alumina 
slurries  showed  some  evidence  of  densification  for  alumina  con¬ 
tents  exceeding  20%.  With  this  experience,  the  composites  had 
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Solution  Yield,  yY  (%) 


Solution  Concentration  (mg  Al^Oj/cm3) 


Fig.  2.  Effects  of  precursor  solution  concentration  on  viscosity  and,  in 
turn,  on  the  degree  of  solution  impregnation  into  a  porous  matrix  com¬ 
posite.  Solid  line  shows  the  predicted  trend  in  weight  gain  assuming  all  of 
the  pores  within  the  matrix  are  filled  by  the  solution. 

been  fabricated  using  a  matrix  slurry  of  80%  mullite  and  20% 
alumina  (80M/20A),  followed  by  precursor  impregnation  and 
pyrolysis.  Consequently,  the  effects  of  the  PDA  on  the  proper¬ 
ties  of  compacts  with  the  latter  composition  (80M/20A)  were 
also  assessed  in  the  present  study. 

The  open  and  closed  porosities  were  determined  from  a 
combination  of  Archimedes  measurements  on  the  compacts20 
and  density  measurements  on  the  constituent  particles  using  a 
pycnometer.  In  some  cases,  the  surface  area  was  measured  using 
the  BET  method.21  The  microstructures  were  probed  in  two 
ways:  by  scanning  electron  microscopy  of  fracture  surfaces,  and 
by  energy-dispersive  spectroscopy  (EDS)  of  foils  in  a  transmis¬ 
sion  electron  microscope  (TEM). 

Young’s  modulus  of  each  fired  compact  was  measured  using 
strain-gauged  flexure  specimens,  5  mm  x  5  mm  x  50  mm  long, 
loaded  in  four-point  bending.  The  inner  and  outer  loading  spans 
were  20  and  40  mm,  respectively.  Toughness  was  obtained  from 
three-point  bend  tests  on  chevron-notched  specimens  with  the 
same  dimensions  and  instrumented  with  an  LVDT  to  measure 
load-point  displacement.  In  cases  where  the  crack  grew  stably 
across  the  specimen  cross-section,  the  toughness  T  was  calcu¬ 
lated  from  the  work  of  fracture.22  In  the  few  instances  where 
fracture  occurred  catastrophically  at  the  load  maximum,  the 
fracture  toughness,  Kc,  was  calculated  from  the  peak  load  and 
the  stress  intensity  solutions  developed  by  He  and  Evans.23  It 
was  then  converted  to  T  through  the  Irwin  relation  using  the 
average  measured  Young’s  modulus  from  the  unnotched  tests. 
Figure  3  shows  comparisons  of  the  results  from  the  two  methods 
for  cases  in  which  stable  growth  had  occurred.  In  the  vast  ma¬ 
jority  of  cases,  the  two  values  are  within  about  30%  of  one  an¬ 
other.  This  difference  coincides  closely  with  the  variability  in  the 
modulus  measurements  for  a  given  material,  indicating  that  the 
two  measurement  techniques  are  consistent  with  one  another. 

III.  Physical  Properties  and  Microstructures 

The  effects  of  precursor  impregnation  and  pyrolysis  on  compact 
porosity  are  summarized  in  Fig.  4.  Also  shown  for  comparison 
is  a  prediction,  assuming  that:  (i)  all  of  the  pore  space  is  filled  by 
the  precursor  solution  on  each  impregnation  cycle;  (ii)  during 
pyrolysis  (900°C),  the  precursor  solution  is  converted  to  y-alu- 
mina;  and  (iii)  during  the  final  heat  treatment  (1200°C),  the 
Y-aiumina  is  converted  to  a-alumina.  The  predicted  porosity  is 
given  by 

0) 


r  (J/m2)  from  work  of  fracture 

Fig.  3.  Comparisons  of  toughness  values  determined  from  work  of 
fracture  and  from  the  peak  load  in  chevron-notched  bend  specimens. 


where  po  is  the  porosity  of  the  compact  before  impregnation  and 
N  is  the  number  of  impregnation  cycles.  Upon  fitting  Eq.  (2)  to 
the  experimental  data,  the  volumetric  yields  of  the  two  precur¬ 
sors  are  inferred  to  be  2.4%  and  6.7%.  These  are  about  20% 
lower  than  the  values  obtained  from  the  precursor  solutions 
alone  (3%  and  8%).  The  discrepancy  may  arise  from  one  or 
more  sources.  Losses  of  precursor  solution  from  the  compact 
may  occur:  notably,  during  the  period  between  the  completion 
of  precursor  impregnation  and  the  onset  of  gelation.  While  no 
significant  closed  porosity  was  detected  in  the  final  compacts,  it 
is  possible  that  closed  porosity  could  have  arisen  during  the 
conversion  of  the  dried  precursor  to  transient  alumina.  In  that 
event,  the  net  amount  of  open  porosity  susceptible  to  impreg¬ 
nation  would  decrease,  and  hence  the  amount  of  alumina 
that  can  be  incorporated  in  each  impregnation  cycle.  This  “tran¬ 
sient”  closed  porosity  could  be  subsequently  reopened  by 
shrinkage  of  the  precursor-derived  oxide  during  the  final  con¬ 
version  to  corundum,  and/or  densified  in  part  as  the  typically 
nanocrystalline  transient  alumina  is  heated  to  higher  tempera¬ 
tures  during  the  conversion  stage. 

Representative  EDS  maps  of  two  prototypical  microstruc¬ 
tures  are  shown  in  Fig.  5.  One  is  of  the  80%  M/20%  A  particle 
mixture  (without  precursor  addition),  and  the  other  comprises 
particles  of  pure  mullite  and  about  1 5%  PDA,  the  latter  coming 
from  10  impregnation  cycles  of  the  concentrated  precursor.  The 
maps  are  produced  by  separately  mapping  out  regions  contain¬ 
ing  A1  and  Si  and  then  superimposing  the  two  images  in  such  a 


Number  of  Impregnations,  N 

Fig.4.  Reduction  in  porosity  following  precursor  impregnation  and 
pyrolysis.  Solid  lines  arc  fits  of  Eq.  (2)  to  the  experimental  data. 
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Fig.  5.  Energy-dispersive  spectroscopy/transmission  electron  micro¬ 
scope  maps  of  thin  foils  of  (a)  80%  mullite/20%  alumina  (without  pre¬ 
cursor  addition)  and  (b)  pure  mullite  compact  impregnated  10  times  with 
the  concentrated  precursor,  yielding  PPDA  =  15.8%. 


way  that  alumina  appears  red  and  mullite  appears  blue.  In  the 
former  case,  the  distinct  nature  of  the  two  particle  types  and 
their  distribution  are  evident;  in  the  latter,  the  alumina  forms  a 
coating  (albeit  somewhat  nonuniform)  on  the  mullite  particles, 
roughly  consistent  with  the  schematic  shown  in  Fig.  1(b). 

The  results  from  a  series  of  surface  area  measurements  show¬ 
ing  the  effect  of  the  PDA  are  plotted  in  Fig.  6.  Here  the  specific 
surface  area  is  normalized  by  the  value.  So,  prior  to  precursor 
impregnation.  The  results  show  a  reduction  in  specific  surface 
area  with  precursor  addition,  but  generally  lower  than  that  ex¬ 
pected  from  simple  geometric  models  based  on  randomly 
packed  spherical  powders.  Two  extreme  scenarios  are  analyzed 
(Appendix  A).  One  assumes  that  the  precursor-derived  oxide  is 
uniformly  distributed  over  the  surfaces  of  the  powders,  building 
up  the  necks  as  its  thickness  increases  (Fig.  A-l(a)).  The  other 
considers  the  precursor  to  segregate  completely  to  the  inter-par¬ 
ticle  contacts,  with  no  material  added  to  the  surface  of  the  pow¬ 
ders  away  from  the  necks  (Fig.  A-l(b)).  The  models  allow  for 
the  presence  of  an  initial  neck  size,  produced  during  the  matrix¬ 
strengthening  treatment  prior  to  the  first  precursor  impregna¬ 
tion.  These  two  scenarios,  albeit  idealized,  provide  useful 
bounds  to  the  reduction  of  pore  surface  area  with  increasing 
densification,  as  well  as  to  the  rate  of  buildup  of  the  necks  be¬ 
tween  particles  (Fig.  6(b)).  Comparison  with  the  experimental 
results  in  Fig.  6(a)  reveals  that,  in  the  early  stages  of  the  process, 
the  normalized  area  falls  within  these  bounds  and  is  closer  to  the 
uniform  distribution  model.  There  is,  however,  evidence  that 
some  amount  of  segregation  occurs,  as  discussed  in  Section  IV. 
The  implication  is  that  the  experimental  points  should  continue 


Relative  Density,  D  =  1-p 


Fig.  6.  Effects  of  precursor-derived  alumina  on  (a)  normalized  specific 
surface  area  and  (b)  neck  size.  Details  of  the  model  predictions  are  given 
in  the  Appendix  A.  (Uncertainty  in  each  measurement  in  (a)  is  less  than 
the  symbol  size.) 

to  fall  within  the  boundaries  as  densification  proceeds,  but  in 
reality  the  pore  area  decreases  at  a  slower  rate  than  that  expect¬ 
ed  from  either  model. 

Further  insight  into  the  sources  of  the  above  discrepancy  can 
be  obtained  from  the  fracture  surfaces  of  the  broken  porous 
compacts  (Fig.  7).  For  moderate  PDA  additions,  the  fracture 
surfaces  are  seemingly  indistinguishable  from  those  of  the  par¬ 
ticle  aggregates  alone;  compare,  for  example,  Fig.  7(a)  and  (b) 
with  (c)  and  (d).  This  result  suggests  that  the  coating  is  thin  and 
uniform,  in  qualitative  agreement  with  the  specific  surface  area 
measurements  for  low  PDA  loadings  in  Fig.  6.  For  larger  PDA 
additions  (Fig.  7(e)  and  (1)),  the  PDA  coating  breaks  up  into 
small  discrete  islands  on  the  larger  mullite  particles.  Such  island 
formation  leads  to  an  increase  in  specific  surface  area  and  hence 
a  lower  rate  of  pore  surface  reduction  with  densification  than 
that  expected  from  the  model  predictions  in  Fig.  6(a). 


IV.  Young’s  Modulus 
(1)  Experimental  Measurements 

The  modulus  measurements  are  plotted  in  Fig.  8  as  a  function 
of  the  percentage  of  total  volume  occupied  by  the  PDA: 
V pda  =  Po~P-  In  this  context,  the  PDA  is  extremely  effective. 
For  instance,  a  10%  addition  of  PDA  to  a  pure  mullite  particle 
aggregate  leads  to  about  a  10-fold  increase  in  modulus.  In  con¬ 
trast,  the  particulate  alumina  in  the  slurry  plays  only  a  minor 
role;  the  results  for  the  two  slurry  compositions — 100%M  and 
80%M/20%A — are  similar  to  one  another. 
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Fig.  7.  Fracture  surfaces  of  three  representative  porous  compacts,  all 
with  pure  mullite  particle  slurry,  but  with  varying  amounts  of  precursor- 
derived  alumina  (PDA):  (a,  b)  FPDA  =  0;  (c,  d)  FPDa  =  8.3%;  (e,  0 
FPda  =  15.8%.  All  had  been  aged  for  2  h  at  1 200"C. 


The  effects  of  the  aging  treatment  are  most  pronounced  in  the 
particle  compacts  that  are  free  of  the  PDA.  For  instance,  for 
Kpda  =  0,  the  modulus  of  the  pure  mullite  compact  increases  by 
a  factor  of  about  3  during  aging,  despite  the  absence  of  meas¬ 
urable  density  change.6  The  change  has  been  attributed  to  sin¬ 
tering  at  the  particle  junctions  by  a  (nondensifying)  surface 
diffusion  mechanism.6  The  large  relative  increase  in  modulus 
results  from  the  very  small  junction  size  (and  hence  modulus)  in 
the  starting  material  such  that  even  modest  amounts  of  junction 
growth  yield  large  increases  in  modulus.  Similar  trends  are  ob¬ 
tained  in  compacts  of  80%M/20%A.  The  relative  magnitude  of 
the  sintering  effect  becomes  progressively  smaller  as  FPDA  in- 


Fig.8.  Summary  of  Young’s  modulus  measurements  and  model  pre¬ 
dictions.  The  predictions  are  based  on  the  results  of  the  discrete  element 
simulations,  shown  in  Fig.  9  and  described  in  the  text.  Strictly,  these 
predictions  pertain  only  to  the  pure  mullite  slurries.  Consequently,  valid 
comparisons  can  only  be  made  between  the  two  solid  lines  and  the  two 
sets  of  filled  symbols. 


creases.  Indeed,  for  FPDA  >  8-10%,  the  effect  is  overshadowed 
by  statistical  variations  in  the  measurements.  This  implies  that 
the  increase  in  modulus  because  of  the  PDA  is  much  greater 
than  that  associated  with  sintering. 

(2)  Simulations  and  Modeling 

Effects  of  the  PDA  on  the  modulus  of  the  compacts  formed 
from  100%  mullite  particles  have  been  computed  using  simula¬ 
tions  based  on  the  DEM.13  In  this  method,  the  material  is 
modeled  as  a  collection  of  particle  elements.  The  junction  re¬ 
sponse  that  defines  the  particle  element  properties  is  derived 
from  a  finite-element  analysis  (FEA)  of  a  single  particle  in  a 
periodic  array.  The  interaction  between  particle  junctions,  char¬ 
acterized  by  the  displacement  of  one  junction  because  of  the 
force  acting  on  another,  is  also  derived  from  the  FEA.  Next,  an 
isotropic  random  aggregate  of  touching  spherical  particles  is 
generated  using  a  computer  algorithm.  To  ensure  equilibrium, 
each  particle  is  required  to  touch  at  least  three  neighbors  upon 
placement  onto  the  existing  aggregate.  The  final  particle  packing 
density  is  about  55%  and  the  average  coordination  number  is 
about  6.  For  monophase  systems,  junction  growth  is  simulated 
by  uniformly  expanding  the  particles  and  redistributing  the 
overlapping  material  uniformly  over  the  free  surface  of  the  par¬ 
ticles.  For  the  two-phase  system  comprising  a  contiguous 
particle  network  and  a  precursor-derived  phase,  the  material  is 
modeled  as  an  aggregate  of  touching  particles,  all  coated  with  a 
uniform  layer,  thickness  h,  of  the  second  phase.  This  idealization 
is  depicted  in  Fig.  A 1 .  The  analyses  described  in  the  Appendix  A 
show  that  the  coating  thickness  can  be  related  to  the  normalized 
junction  size  a  and  the  volume  fraction  of  PDA.  The  predicted 
values  are  plotted  in  Fig.  6(b)  for  a  coordination  number  of  6 
and  two  values  of  the  initial  neck  size  (oq,  =  0,  0. 1).  To  first  or¬ 
der,  the  junction  size  is  given  by 


2KroA  +Q i 


3(1  —  po) 


(3) 


For  the  two-phase  system,  the  elastic  response  of  the  junc¬ 
tions  is  calculated  using  a  periodic  array  FEM  model,  with  the 
two  phases  explicitly  discretized.  For  both  mono-  and  two-phase 
systems,  the  random  particle  network  is  then  subjected  to  a  pre¬ 
scribed  macroscopically  uniform  strain  field  and  the  effective 
elastic  response  is  determined  using  the  DEM.13  In  the  present 
calculations.  Young’s  moduli  E  and  Poisson’s  ratios  v  of  the 
two  phases  are  taken  to  be  £A  =  400  GPa,  EM  =  200  GPa,  and 
vA  =  vM  =  0.2,  where  the  subscripts  A  and  M  refer  to  alumina 
and  mullite,  respectively. 

Figure  9  shows  the  results  of  these  simulations  for  the  two- 
phase  system  (with  oq>  =  0)  as  well  as  for  monophase  aggregates 


Fig.  9.  Discrete  element  method  (DEM)  simulations  of  Young’s  mod¬ 
ulus  of  monophase  particle  aggregates  (either  alumina  or  mullite)  and 
mullite  particle  aggregates  with  a  uniform  coating  of  alumina. 
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of  pure  mullite  and  pure  alumina,  plotted  against  a.  The  latter 
results  are  well  described  by  the  empirical  formula 

■—  =  0.4a(l  +  14  a3)  (4) 

where  E%  is  the  modulus  of  the  solid.  Moreover,  the  modulus  of 
the  two-phase  aggregate  can  be  estimated  by  the  volume  frac¬ 
tion-weighted  average  of  the  moduli  of  the  two  monophase  ag¬ 
gregates,  evaluated  at  the  same  value  of  a.  This  estimate  is 
obtained  by  replacing  Es  in  Eq.  (4)  with  the  volume-weighted 
average  of  the  moduli  of  the  two  solid  phases.  Upon  calculating 
the  volume  ratios  from  geometry,  the  estimated  modulus  of  the 
two-phase  system  becomes 

E  =  0.4a(l  +  14a3)  [£mom+£a(1 -«m)]  (5) 

where 


Fig.  6(a).  Furthermore,  although  this  prediction  is  strictly  ap¬ 
plicable  to  the  compacts  containing  pure  mullite  particles,  it  ap¬ 
pears  to  underestimate  only  slightly  the  results  for  the  80M/20A 
mixtures. 

The  effects  of  the  extended  (1000  h)  aging  treatment  on  the 
modulus  of  compacts  containing  the  PDA  are  treated  in  the 
following  way.  In  an  approximate  sense,  the  addition  of  PDA 
can  be  viewed  as  being  equivalent  to  a  sintering  treatment,  pro¬ 
vided  both  yield  the  same  changes  in  junction  size  without  af¬ 
fecting  interparticle  spacing.  This  equivalence  stems  from  the 
fact  that  the  modulus  is  dictated  by  the  current  junction  size, 
essentially  independent  of  the  route  taken  to  achieve  it  (not¬ 
withstanding  slight  differences  in  the  moduli  of  the  two  phases). 
In  the  sintering  case,  the  junction  size  evolves  according  to  a 
power  law  of  the  form:6 


1  ~Po 

M  1  —  Po  +  k'pDA 


(6) 


Upon  calculating  FPDA  as  a  function  of  a  from  Eq.  (3),  with 
oto  =  0,  and  combining  Eqs.  (4),  (5),  and  (6),  the  estimated  mod¬ 
ulus  of  the  two-phase  system  becomes 


E  =  0.4a(l  +  14  a3)  Em 


( 1  ~  Po)  +  (£a/ Em  )  V pda 
(I  —  po)  +  F PDA 


(7) 


where  t  is  the  sintering  time,  tR  is  a  reference  time,  and  n  is  a 
constant,  the  latter  two  depending  on  the  transport  mechanism. 
(For  mullite,  junction  growth  occurs  by  surface  diffusion, 
whereupon  n  =  7.6)  The  junction  size  of  that  results  from  pre¬ 
cursor  addition  (according  to  Eq.  (8))  can  be  viewed  as  having 
been  developed  through  a  sintering  treatment  of  the  monophase 
particle  aggregate,  for  an  equivalent  sintering  time,  tc q,  given  by 

<cq  =  '*(*')"  (10) 


The  modulus  estimates  for  the  mullite/alumina  system  ob¬ 
tained  by  this  averaging  procedure,  shown  by  the  dashed  line  in 
Fig.  9,  agree  well  with  the  results  of  the  DEM  simulations. 

The  model  results  are  replotted  in  Fig.  8,  for  comparison  with 
the  experimental  measurements.  Several  features  are  notewor¬ 
thy.  Because  the  junctions  between  the  mullite  particles  are 
modeled  as  point  contacts  when  ao  =  0,  the  predicted  modulus 
goes  to  zero  when  FPDA  =  0.  This  behavior  differs  from  that 
obtained  experimentally.  That  is,  the  modulus  of  the  mullite 
particle  aggregate  is  finite,  at  about  10  GPa:  a  consequence  of 
the  finite  (but  small)  junction  size.  The  predicted  modulus  sub¬ 
sequently  increases  with  KPDA,  but  at  a  rate  that  is  significantly 
lower  than  that  of  the  measurements.  The  rather  large  discrep¬ 
ancy  at  large  values  of  FPda  (approaching  60  GPa)  suggests  that 
other  sources  of  error  must  be  present,  apart  from  the  neglect  of 
the  initial  junction  size.  It  is  surmised  that,  because  of  capillarity, 
the  PDA  forms  preferentially  at  the  mullite  particle  contacts, 
rather  than  as  a  uniform  particle  coating,  thereby  stiffening  the 
particle  junctions  and  elevating  the  composite  modulus. 

The  effects  of  both  the  initial  junction  size  and  the  nonuni¬ 
formity  of  the  PDA  can  be  taken  into  account  through  an  ap¬ 
proximate  modification  to  the  DEM  results.  For  this  purpose, 
Eq.  (3)  is  modified  to  read  as 


After  additional  sintering  for  time  t,  the  final  junction  size  is 
evaluated  from  Eq.  (9),  using  the  total  sintering  time,  l+tcq.  The 
result  is 


(11) 


The  model  prediction  (Eqs.  (7)  and  (11))  is  plotted  in  Fig.  8, 
using  a  reference  time  /R  =  2x  106  h.  This  adequately  captures 
the  trend  for  the  pure  mullite  particle  aggregates  with  added 
PDA.  It  is  of  interest  to  note  that  the  corresponding  reference 
time  inferred  from  sintering  treatments  of  pure  mullite  (absent 
PDA)  is  ir  =  8x  106  h:6  the  fourfold  difference  reflecting  the 
higher  sinterability  of  the  alumina  in  the  PDA-containing  mix¬ 
tures.  Indeed,  it  is  for  this  reason  that  the  present  model  pre¬ 
dictions  (using  the  lower  value  of  fR)  slightly  overestimate  the 
modulus  for  the  pure  mullite  (FPDA  =  0). 

In  light  of  the  rather  weak  sensitivity  of  the  modulus  to  the 
slurry  composition  (I00M  vs  80M/20A),  no  attempt  was  made 
to  explicitly  model  a  two-phase  particle  aggregate  containing  a 
precursor-derived  phase.  Models  for  the  behavior  of  the  two- 
phase  aggregate  alone  are  presented  elsewhere.6 


;  a2  2PPda  .  o  to\ 

+  <  <!> 

where  P  is  a  nondimensional  parameter,  >  1,  that  accounts  for 
preferential  alumina  formation  at  the  particle  junctions.  In  this 
context,  the  normalized  junction  size  a  in  both  Fig.  9  and  Eq.  (8) 
is  reinterpreted  as  of.  Combining  Eqs.  (7)  and  (8)  then  yields  the 
variation  in  E  with  KPda-  The  resulting  predictions  for  P  =  1.2 
and  o<o  =  0. 1  are  shown  in  Fig.  8.  This  combination  yields  a  good 
fit  to  the  data  for  the  2  h  aging  over  the  entire  range  of  FPda- 
The  corresponding  buildup  in  neck  size  is  compared  in  Fig.  6(b) 
with  those  resulting  from  the  two  extreme  scenarios  analyzed  in 
the  Appendix  A.  It  is  noted  that  the  inferred  level  of  segregation, 
reflected  in  the  value  of  p,  is  modest.  The  neck  size  curve  in 
Fig.  6(b)  is  closer  to  the  uniform  distribution  scenario  than 
to  the  full  segregation  model,  in  qualitative  agreement  with 
the  specific  area  measurements  discussed  in  the  context  of 


V.  Toughness 

The  toughness  measurements  are  plotted  in  Fig.  10.  For  the  as- 
processed  materials,  the  toughness  increases  approximately  lin¬ 
early  with  Vpda,  almost  independent  of  the  slurry  composition 
(100M  vs  80M/20A).  Furthermore,  the  toughness  increases 
slightly  following  the  additional  1000  h  aging.  These  results 
are  consistent  with  the  associated  increases  in  junction  size  be¬ 
cause  of  both  addition  of  PDA  and  sintering. 

Also  shown  in  Fig.  10  are  the  predictions  of  numerical  sim¬ 
ulations  of  fracture,  again  using  the  DEM.  The  results  of  the 
simulations  for  monophase  aggregates  are  well  described  by:6 

f=12«2  (12) 

1  j 

where  Tj  is  the  intrinsic  junction  toughness.  This  result  is  adapt¬ 
ed  to  the  PDA-containing  aggregates  by  using  the  junction  size 
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Fig.  10.  Summary  of  toughness  measurements  and  model  predictions. 
The  model  predictions  are  based  on  the  parameter  values  used  to  cal¬ 
culate  Young’s  modulus  in  Fig.  8. 

given  in  Eq.  (8),  whereupon  the  toughness  becomes 


2FppA 

3(1  -Po) 


(13) 


For  consistency  with  the  modulus  predictions,  simulations  for 
the  2  h  aging  treatment  were  performed  using  cto  =  0.1  and 
P  =  1.2,  with  Tj  being  a  fitting  parameter.  This  procedure  yields 
Tj  =  4  J/m2:  only  slightly  greater  than  the  value  (3  J/m2)  that  had 
been  inferred  previously  from  measurements  on  pure  mullite.6 

The  effects  of  extended  aging  are  taken  into  account  in  the 
same  manner  described  above  for  the  modulus.  That  is,  junction 
growth  is  assumed  to  follow  in  accordance  with  Eq.  (11),  with  a 
reference  time,  /R  =  2  x  1 06  h.  Combining  this  result  with 
Eq.  (12)  yields  the  predicted  curve  in  Fig.  10.  Although  the  pre¬ 
diction  yields  the  correct  qualitative  trends,  the  scatter  in  the 
experimental  measurements  precludes  a  more  critical  assessment 
of  the  model.  Indeed,  some  discrepancy  is  expected,  since  the 
pertinent  reference  time  should  depend  on  the  amounts  of  the 
two  phases  present  at  the  particle  junctions,  as  dictated  by 
the  PDA  concentration. 


VI.  Crack  Deflection 

The  measured  values  of  modulus  and  toughness  have  been  used 
to  compute  the  crack  deflection  parameter  £  (Eq.  (1)),  assuming 


PDA  concentration,  VPDa  (%) 


Fig.  11.  Effects  of  precursor-derived  alumina  (PDA)  concentration  and 
aging  on  the  crack  deflection  parameter,  X,  calculated  using  fiber  prop¬ 
erties  rr=  15  J/m2  and  Er  =  260  GPa. 


that  the  aggregates  are  combined  with  Nextel  720  fibers  to  form 
a  composite.  The  results  are  plotted  against  Kpda  in  Fig.  1 1 . 
Here,  again,  the  two  slurry  compositions  yield  similar  results. 
For  the  as-processed  materials,  £  decreases  with  increasing 
Vpda.  and  eventually  falls  below  the  critical  value,  aw  I,  at 
Fpda  ~  9%.  This  point  is  expected  to  mark  the  onset  of  crack 
penetration  into  the  fibers  and  essentially  a  complete  loss  in 
damage  tolerance.  Extended  aging  causes  a  reduction  in  £,  es¬ 
pecially  at  low  values  of  KpDA.  Furthermore,  it  reduces  the  crit¬ 
ical  value  associated  with  crack  penetration,  to  Kpda  ~  7-8%. 

The  predicted  results  from  the  models  for  modulus  and 
toughness  (calibrated  accordingly)  are  also  shown  in  Fig.  11. 
A  satisfactory  correlation  with  the  experimental  measure¬ 
ments  is  obtained.  Most  notably,  the  predicted  critical  value, 
Kpda  «  9-10%,  is  broadly  in  line  with  that  inferred  from  the 
measurements. 


VII.  Conclusions 

The  mechanical  properties  of  mullite  particle  aggregates  can  be 
effectively  tailored  by  the  addition  of  PDA.  Tenfold  changes  in 
both  modulus  and  toughness  can  be  obtained  through  the  ad¬ 
dition  of  moderate  amounts  of  alumina,  in  the  range  of  10-15%. 
Significantly  smaller  changes  can  be  effected  through  the  use  of 
alumina  particle  additions.  The  precursor-derived  route  has  the 
additional  advantage  over  the  mixed  mullite/alumina  particle 
route  in  that  the  contiguity  of  the  mullite  particle  network  is  not 
compromised.  Consequently,  the  structure  is  relatively  immune 
from  densification  at  the  targeted  upper  use  temperature  for  the 
corresponding  fiber  composites. 

The  mechanical  properties  of  the  two-phase  aggregates  can  be 
described  using  idealizations  of  the  particle  and  precursor  to¬ 
pology  coupled  with  simulations  based  on  the  DEM.  Compar¬ 
isons  between  predicted  and  measured  properties  suggest  that 
the  precursor-derived  phase  exhibits  some  degree  of  segregation 
to  the  particle  junctions,  yielding  larger  property  elevations  rel¬ 
ative  to  that  for  a  uniform  precursor  distribution. 

Limits  have  been  established  on  the  maximum  acceptable 
concentration  of  the  precursor-derived  phase,  such  that  a  com¬ 
posite  with  the  candidate  matrix  would  satisfy  the  crack  defec¬ 
tion  condition.  In  order  to  produce  composites  that  are  damage 
tolerant  under  fiber-dominated  loadings  yet  retain  adequate  off- 
axis  properties,  precursor  concentrations  slightly  below  the  crit¬ 
ical  value  would  be  preferred.  Experiments  are  in  progress  to 
assess  these  predictions. 


Appendix  A:  Models  for  Precursor  Buildup  on  Particles 
(A.l)  Preliminaries 

Models  representing  two  extreme  scenarios  of  precursor  buildup 
around  the  matrix  particles  are  developed.  The  geometrical  fea¬ 
tures  of  the  models  are  depicted  in  Fig.  Al.  The  goal  is  to  de¬ 
scribe  the  variation  in  specific  surface  area  (defined  as  pore  area 
per  unit  mass  or  particles)  and  contact  area  between  particles  as 
the  relative  density  increases  with  the  addition  of  precursor-de¬ 
rived  ceramic.  The  particles  are  assumed  to  be  initially  spherical 
with  radius  R  and  packed  in  a  random  configuration  with  co¬ 
ordination  number  z.  The  models  include  the  possibility  of  a  fi¬ 
nite  contact  area  of  radius  a0,  developing  during  the  initial  heat 
treatment,  prior  to  precursor  impregnation.  In  that  event,  the 
initial  surface  area  and  volume  per  particle  are 


(A-la) 


(A- lb) 


where  ho  is  the  height  of  the  spherical  cap  that  is  truncated  from 
the  particle  to  create  an  initial  neck  of  radius  a0  (Fig.  Al). 
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(a)  Uniform  (b)  Segregated 

Fig.Al.  Schematics  of  the  two  extreme  scenarios  in  precursor  topolo¬ 
gy:  (a)  uniform  distribution  on  the  particle  surfaces,  and  (b)  complete 
segregation  to  the  particle  necks. 


In  general,  the  relative  (specific)  surface  area  is  given  by  the 
expression 


S__A_ 
So  Aq 


(A-2) 


where  p0  and  pA  are  the  densities  of  the  particles  and  the  PDA, 
respectively,  and  A  and  AV  are  the  pore  area  and  volume  of 
PDA  added  per  particle,  respectively.  The  porosity  is  given  by 


P=Po-{\-Po)-tt  (A-3) 


where  p0  is  the  initial  porosity.  The  normalized  neck  radius  is 
defined  as 


a 

Rett 


(A-4) 


where  is  the  current  radius  of  the  particle  away  from  the 
neck  area.  Its  initial  value,  before  precursor  impregnation,  is 


( A. 2)  Uniform  Distribution 

In  one  scenario  the  PDA  is  assumed  to  be  distributed  uniformly 
over  the  free  surface  of  the  particles  as  a  coating  with  thickness 
h,  except  where  prohibited  by  particle  contact  (see  Fig.  A  1(a)). 
The  analysis  is  based  on  a  single  representative  particle  and  its 
corresponding  coating  and  void  space.  The  pore  surface  area 
and  PDA  volume  per  particle  are  given  by 


f)  + 

1  2  \i 

V 

L  2  \  r) 

(sH'-a®’} 


(A-6a) 
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(A-6b) 

The  porosity  is  related  to  the  normalized  thickness  h/R  through 
Eqs.  (A-3)  and  (A-6b).  For  this  geometry  Rcn=  R+h  and  the 
normalized  neck  radius  is 


(A-7) 


(A-8) 


(A.3)  Precursor  Segregation  to  Necks 
In  a  diverging  scenario,  the  precursor  is  assumed  to  be  fully 
segregated  to  the  necks.  For  simplicity,  the  PDA  buildup  is  tak¬ 
en  to  have  a  cylindrical  geometry  as  illustrated  in  Fig.  A  1(b). 
The  half  thickness  at  the  edge  of  the  neck  is  designated  as  k  and 
is  analogous  but  not  equivalent  to  the  h  used  in  the  earlier  case. 
The  pore  surface  area  and  PDA  volume  per  particle  are  now 
given  by 


A  =  4nR2 


[  2  \R)  2\r) 


(A-9a) 


AV  =  -n/?3 


{jO’K'-S) -’(*)]}  "*> 


As  before,  the  porosity  is  related  to  the  normalized  thickness 
k/R  through  Eqs.  (A-3)  and  (A-9b).  For  this  geometry,  /?efr  = 
R  and  the  normalized  neck  radius  is 


(A- 10) 


Combining  with  Eq.  (A-5)  yields 


a  = 


a®  +  ^\/1  ~°t®~ 


(A- 11) 


Comparison  of  the  models  is  performed  by  setting  a  value  of 
p0  and  then  calculating  the  values  of  h/R  and  k/R  implicitly  from 
the  combination  of  Eq.  (A-3)  and  either  Eq.  (A-6b)  or  (A-9b). 
Once  the  pertinent  value  of  the  normalized  thickness  is  known, 
S/S0  and  a  can  be  readily  calculated.  The  results  of  such  calcu¬ 
lations  are  compared  with  the  experimental  measurements  in 
Fig.  6. 
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The  present  study  explores  the  effects  of  exposure  in  a  laboratory  combustor  on  microstructural  stability  and  property 
retention  of  an  all-oxide  fiber-reinforced  ceramic  composite.  The  material  consists  of  a  porous  mullite-alumina  matrix  and 
Nextel  720  fibers  in  an  eight-harness  satin  weave.  To  assess  the  effects  of  matrix  strength,  two  matrix  conditions  are  used, 
distinguished  from  one  another  by  the  amount  of  alumina  added  through  precursor  impregnation  and  pyrolysis  (1.8%  and 
4.8%).  In  both  cases,  the  dominant  damage  mode  upon  exposure  involves  interply  delamination  along  the  panel  midplane. 
However,  significant  reductions  in  the  rate  and  extent  of  cracking  are  obtained  in  the  material  with  higher  alumina  content:  a 
result  of  the  higher  delamination  resistance.  Mechanical  tests  performed  on  exposed  specimens  reveal  a  slight  (10-20%) 
reduction  in  tensile  strength  along  the  fiber  direction  and  a  comparable  increase  in  shear  strength.  These  trends  suggest  some 
sintering  of  the  matrix  upon  exposure.  Examinations  of  fracture  surfaces  provide  additional  supporting  evidence.  Implications 
for  long-term  performance  and  strategies  for  imparting  improvement  in  microstructural  stability  and  delamination  resistance 
are  discussed. 


Introduction 

The  ever-increasing  demands  on  propulsion  sys¬ 
tems  for  aircraft  engines  have  spurred  the  development 
of  continuous  fiber  ceramic  composites  (CFCCs).  Ma¬ 
terials  systems  based  on  all-oxide  constituents  are  of 
particular  interest  for  applications  requiring  long  service 
lives  in  oxidizing  environments.  Although  a  number  of 
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these  have  been  shown  to  exhibit  good  stability  in  am¬ 
bient  air  at  elevated  temperatures,1^*  their  performance 
in  the  targeted  service  environments  has  yet  to  be  crit¬ 
ically  assessed.  The  purpose  of  the  current  article  is  to 
present  the  results  of  an  experimental  investigation  of 
the  stability  of  one  particular  oxide  CFCC  upon  expo¬ 
sure  in  a  combustion  environment. 

The  material  of  interest  is  a  porous-matrix  com¬ 
posite  without  a  fiber  coating.5-8  It  has  been  shown  to 
exhibit  excellent  damage  tolerance  provided  the  modu¬ 
lus  and  the  toughness  of  the  matrix  are  not  excessively 
high.9,10  Otherwise,  cracks  formed  within  the  matrix 
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penetrate  into  the  fibers  without  deflection  at  the  inter¬ 
vening  interface.  The  main  matrix  constituent  is  mul- 
lite,  chosen  for  its  resistance  to  sintering  at  the  upper  use 
temperature  of  the  fibers  (about  1200°C).  The  minor 
constituent  is  alumina,  present  both  in  the  form  of  fine 
particulates  as  well  as  through  pyrolysis  of  a  precursor 
liquid.  The  composite  has  demonstrated  good  retention 
in  mechanical  properties  following  long  (lOOOh) 
exposures  at  elevated  temperatures  (to  1200°C)  in 
ambient  air. 

In  the  current  study,  annular  test  specimens  of  the 
composite  are  exposed  in  a  laboratory-scale  trapped- 
vortex  combustor  rig  under  conditions  that  simulate 
the  pertinent  heat  loads  and  combustion  environ¬ 
ments.1  1-14  Specifically,  the  combustor  is  operated 
using  liquid  JP-8+100  jet  fuel  and  the  peak  material 
temperature  is  slightly  above  1200°C.  Coupled  with 
complementary  mechanical  tests  and  microstructural 
examinations,  such  tests  yield  insights  into  degradation 
and  damage  mechanisms  associated  with  the  exposure, 
and  thus  provide  guidance  for  future  research  activities 
in  composite  development. 

Materials 

The  reinforcements  were  1500-denier  Nextel  720 
fiber  tows,  woven  into  an  eight-harness  satin-weave 
cloth.  Each  composite  panel  consisted  of  nine  fiber  lay¬ 
ers  and  a  porous  matrix  of  ninemullite  and  alumina. 
The  average  fiber  volume  fraction  was  38%  and  the 
panel  thickness  was  2.6  mm.  A  typical  cross-section 
through  one  such  panel  is  shown  in  Fig.  1. 

The  processing  route  was  essentially  identical  to 
that  described  in  earlier  articles.5'8  Briefly,  the  matrix 


Fig.  1.  Scanning  electron  micrograph  of  polished  cross-section  of 
pristine  f. unexposed)  composite  (4.8%  PDA). 


was  introduced  in  two  steps:  (i)  vacuum  infiltration  of  a 
slurry  containing  80%  of  1  pm  mullite  and  20%  of 
0.2  pm  alumina  particles  into  the  fiber  perform,  and  (ii) 
impregnation  and  pyrolysis  of  an  alumina  precursor  so¬ 
lution  (Al2Cl(OH)5).  The  latter  impregnation  and 
pyrolysis  process  was  performed  twice.  Two  precursor 
solution  concentrations  were  used,  as  described  below. 
To  prevent  segregation  during  drying,  the  panels  were 
gelled  by  placing  them  in  a  chamber  containing  gaseous 
NH4OH  for  6h  immediately  following  solution  im¬ 
pregnation.  After  drying,  the  panels  were  fired  at  900°C 
for  2  h.  Following  completion  of  both  impregnation 
cycles,  the  panels  were  given  a  final  heat  treatment  at 
1200°C  for  2h. 

Selection  of  the  precursor  solution  concentration 
was  based  in  pan  on  the  results  of  concurrent  studies, 
focused  on  the  effects  of  the  volume  fraction  of  precur¬ 
sor-derived  alumina,  VpDA,  on  both  the  matrix  proper¬ 
ties,  and  the  strength  and  notch  sensitivity  of  the 
corresponding  composites.'5  The  results  indicate  that 
the  conditions  for  the  deflection  of  matrix  cracks  at  the 
fiber-matrix  interface  are  satisfied  when  Vp[>A<9%10 
At  higher  concentrations,  matrix  cracks  penetrate  into 
the  fibers,  resulting  in  complete  loss  in  damage  toler¬ 
ance.  This  sets  limits  on  the  allowable  range  of  VPqA. 
Moreover,  within  the  deflection  domain  (0%< 
KPDA<9%),  the  in-plane  notched  strength  decreases 
with  increasing  VpDA:  a  consequence  of  reduced  fiber 
pullout  and  reduced  fracture  energy.15  In  contrast,  the 
delamination  resistance  increases  with  VpoA-7  Because 
of  these  divergent  trends,  the  optimal  precursor  con¬ 
centration  is  somewhat  ill-defined.  Consequently,  in  the 
present  study,  two  different  panels  were  fabricated  and 
tested,  distinguished  from  one  another  by  the  precursor 
solution  concentration.  Following  two  impregnation 
and  pyrolysis  cycles,  the  two  solutions  yielded  PDA 
concentrations  of  1.8%  and  4.8%.  Although  both  values 
fall  within  the  crack  deflection  domain,  the  panel  with 
higher  PDA  concentration  is  expected  to  exhibit  supe¬ 
rior  delamination  resistance  and  inferior  in-plane  dam¬ 
age  tolerance,  and  vice  versa. 

Exposure  and  Testing  Methods 

Composite  specimens  were  exposed  in  a  laboratory- 
scale  ultra-compact  combustor  rig,  shown  in  Fig.  2.  The 
rig  has  been  designed  to  emulate  an  inter-stage  turbine 
burner.  Such  burners  are  being  targeted  for  use  in  future 
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Fig.  2.  (a)  Schematic  of  longitudinal  cross-section  of  ultra¬ 

compact  combustor  rig.  (b)  The  combustor  rig,  disassembled  to 
reveal  the  internal  parts  of  the  combustion  chamber.  (For  scale,  the 
inner  radius  of  the  combustor  cavity  is  59  mm.) 


propulsion  systems,  for  the  purpose  of  improving 
specific  thrust.  The  rig  has  also  been  used  as  a  test-bed 
for  assessment  of  new  materials  for  gas  turbine  engines, 
notably  ceramic  composites.  Its  design,  operation,  and 
performance  are  described  in  detail  elsewhere. 11-1 4 
Briefly,  the  rig  consists  of  an  annular  combustor  cavity 


with  six  fuel  and  24  air  injection  ducts,  equally  spaced 
around  the  circumference.  The  ducts  drive  a  high-g 
swirling  flow  within  the  cavity.  Here,  the  fuel  mixture  is 
partially  burned  before  being  transported  to  the  axial 
main  air  zone.  The  mixture  is  subsequently  diluted  and 
the  reaction  continues  to  completion.  Test  specimens 
for  exposure  in  the  rig  are  in  the  form  of  flat  annular 
rings,  121  mm  outer  diameter  and  83  mm  inner  diam¬ 
eter.  When  inserted  into  the  rig,  they  cover  the  entire  aft 
face  of  the  combustion  cavity.  Attachment  is  accom¬ 
plished  via  a  window-frame  assembly.  The  test  speci¬ 
mens  from  the  two  composite  panels,  with  Vpn>A  = 
1.8%  and  4.8%,  are  henceforth  referred  to  as  ring  #1 
and  #2,  respectively. 

Prior  to  insertion  into  the  rig,  the  test  specimens 
were  coated  with  a  temperature-sensitive  paint  (MC  490- 
10,  Thermographic  Measurements  Co.  Ltd.,  Flintshire, 
UK)  The  paint  was  sprayed  on  and  subjected  to  a  pre¬ 
treatment  of  30  min  at  250°C  to  remove  the  solvent.  The 
paint  undergoes  ten  color  changes  between  490°C  and 
1250°C.  For  calibration,  a  series  of  small  coated  test 
specimens  had  been  subjected  to  prescribed  temperatures. 
Based  on  this  calibration  and  the  paint  colors  following 
exposure,  the  peak  surface  temperature  in  the  combustor 
rig  was  estimated  to  be  1 200-1 225°C  in  localized  hot 
spots  around  each  fuel  nozzle.  These  results  are  in  close 
agreement  with  the  predictions  by  Katta  et  al. 1 6 

The  rig  was  typically  operated  uninterrupted  for 
periods  of  1 — 2  h,  using  liquid  JP-8+100  jet  fuel  with 
essentially  constant  airflow  rate.  To  evaluate  combustor 
performance,  a  variety  of  test  conditions  were  employed, 
including  rich  burn,  lean  burn,  and  flame  blow-out. 
Periodically,  the  test  specimens  were  extracted  from  the 
combustor  rig  and  examined  in  an  optical  microscope. 
These  examinations  were  used  to  guide  decisions  re¬ 
garding  further  combustor  rig  exposure.  For  instance,  in 
the  case  of  ring  #1,  numerous  delamination  cracks  were 
found  emanating  from  the  inner  ring  surface  following 
20  h  of  exposure.  Consequently,  the  exposure  experi¬ 
ments  were  discontinued  and  the  ring  then  sectioned 
for  microstructural  examination.  In  contrast,  ring  #2 
exhibited  minimal  delamination  following  a  similar  ex¬ 
posure  and  hence  the  rig  exposure  test  was  continued. 
After  a  considerably  longer  exposure  time,  small  delam¬ 
ination  cracks  were  found,  including  on  the  inner  ring 
surface.  The  exposure  test  on  this  ring  was  terminated 
after  86  h. 

Following  exposure,  the  rings  were  examined  in  an 
optical  microscope  to  identify  the  dominant  damage 
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modes.  Additionally,  measurements  were  made  of  the 
matrix  porosity  (via  the  Archimedes  method17)  to  assess 
potential  changes  associated  with  sintering.  Finally, 
measurements  of  plate  thickness,  coupled  with  scanning 
electron  microscope  (SEM)  examinations  of  near¬ 
surface  regions  of  polished  cross-sections,  were  used  to 
determine  whether  any  material  had  been  volatized  as  a 
consequence  of  the  interaction  with  the  combustion 
environment. 

To  further  assess  the  effects  of  combustor  rig  expo¬ 
sure,  mechanical  tests  were  performed  to  characterize 
both  composite  and  matrix  properties.  This  part  of  the 
study  was  restricted  to  the  composite  used  to  make  ring 
#2,  i.e.,  the  one  with  the  higher  PDA  concentration, 
because  of  the  extensive  delamination  in  the  other  ring. 
Comparisons  were  made  with  the  corresponding  prop¬ 
erties  of  the  pristine  (as-processed)  material. 

The  composite  strength  was  measured  using 
straight  tensile  bars  cut  from  the  rings,  in  the  manner 
illustrated  in  Fig.  3.  Each  had  a  width  of  10  mm,  a 
gauge  length  of  35  mm,  and  was  oriented  with  one  of 
the  fiber  directions  aligned  with  the  loading  axis.  Lon¬ 
gitudinal  strains  were  measured  through  a  strain  gauge. 
The  fracture  surfaces  were  subsequently  examined  in  an 
SEM.  The  shear  strength  was  measured  using  losipescu 
tests,  including  in  the  0°/90°  orientation  (Fig.  3). 
Strains  were  measured  using  gauges  at  +45°/— 45°. 


Fiber 

directions  Tensile 


Fig.  3.  Schematic  of  the  ring  specimens  used  for  combustor  rig 
exposure  and  the  locations  from  which  the  mechanical  test  specimens 
were  machined. 


The  matrix  strength  was  characterized  through 
Vickers  hardness  measurements.  For  this  purpose,  small 
specimens  were  machined  from  the  panels  (both  pristine 
and  exposed)  and  subsequently  polished  to  a  0.1  pm 
finish.  To  facilitate  polishing,  the  specimens  were  first 
impregnated  with  a  low  viscosity  epoxy.  This  epoxy  was 
burned  out  after  polishing  by  heating  the  specimens 
slowly,  at  5°C/min,  to  550°C,  holding  at  the  temper¬ 
ature  for  2  h,  and  cooling  slowly  to  room  temperature. 
Vickers  indentations  were  made  in  two  location  types: 

(i)  within  matrix-rich  regions  between  fiber  cloths,  and 

(ii)  at  the  centers  of  “end-on”  tows.  In  both  cases,  the 
indent  size  was  restricted  to  one-half  of  the  pertinent 
available  space,  i.e.,  either  the  width  of  the  matrix  region 
or  the  width  of  the  fiber  tow.  This  was  accomplished  by 
selecting  the  indentation  load  accordingly.  Specifically, 
the  loads  were  200  and  500  g  for  the  matrix-rich  regions 
and  the  fiber  tows,  respectively. 

The  same  polishing  procedure  was  used  to  prepare 
specimens  for  examination  of  damage  mechanisms.  The 
latter  were  performed  in  an  SEM. 

Measurements  and  Observations 

The  only  visible  microstructural  damage  following 
combustor  rig  exposure  was  in  the  form  of  delamination 
cracks,  most  emanating  from  the  inner  ring  surface 
(Figs.  4-6);  there  was  no  evidence  of  fiber  damage. 
The  cracks  were  readily  observed  with  the  unaided  eye  in 


Fig.  4.  Macrophotograph  of  the  inner  surface  of  ring  #2  following 
86 h  of  rig  exposure.  Arrows  highlight  short  segments  of  delamina¬ 
tion  cracks  along  the  panel  midplane. 
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Fig.  5.  Typical  delamination  cracks  in  ring  #1  (V pda  =  1.8%, 
20  h  exposure),  at  (a)  inner  and  (b)  outer  surfaces.  Inset  in  (a) 
shows  the  crack  following  an  interply  path  as  well  as  the  extent  of 
crack  opening  (  ~  lOOflm). 

ring  #1.  In  contrast,  despite  the  longer  exposure  time, 
the  delamination  cracks  in  ring  #2  were  considerably 
smaller  and  only  seen  upon  inspection  in  an  optical  mi¬ 
croscope  (Fig.  4).  Micrographs  of  cross-sections  also  re¬ 
vealed  periodic  transverse  matrix  cracks  (Figs.  5  and  6):  a 
consequence  of  drying  during  processing.  This  type  of 
cracking  was  also  seen  in  the  pristine  material  (Fig.  1). 

The  extent  of  delamination  was  significantly  dif¬ 
ferent  in  the  two  test  specimens.  For  ring  #1,  the  cracks 
penetrated  to  a  depth  of  about  10  mm,  with  crack  open¬ 
ing  displacements  at  the  inner  ring  surface  exceeding 
100  pm.  In  contrast,  for  ring  #2,  the  cracks  typically 
penetrated  only  about  1—2  mm,  with  correspondingly 
lower  crack  opening  displacements  (typically  about 
30  pm).  This  difference  is  attributable  to  the  increase 
in  matrix  strength  resulting  from  the  additional  PDA 
content  in  ring  #2  and  its  effect  on  delamination  resist- 


Fig.  6.  Delamination  in  ring  #2  (V pda  =  4.8%,  86 h  exposure). 
In  contrast  to  ring  #1,  the  extent  of  cracking  is  diminished,  as 
characterized  by  crack  length  and  opening  displacement. 


Fig.  7.  Tensile  response  of  pristine  and  exposed  materials. 

ance.  This  assertion  is  supported  by  the  matrix  hardness 
values:  30  +  8  and  101  +  17kg/mm2  for  rings  #1  and 
#2,  respectively.  Similar  increases  in  hardness  and  cor¬ 
responding  increases  in  toughness  have  been  found  in 
neat  matrix  specimens  of  the  same  compositions. 

The  in-plane  tensile  strength  of  ring  #2  decreased 
slightly  following  exposure,  by  about  10—20%  (Fig.  7). 
This  reduction  can  be  rationalized  on  the  basis  of  ob¬ 
servations  of  the  fracture  surfaces  (Fig.  8).  That  is,  the 
extent  of  fiber  pullout  within  individual  tows  was  no¬ 
ticeably  greater  in  the  pristine  material  than  that  in  the 
exposed  material.  This  indicates  an  increase  in  the 
strength  of  the  fiber-matrix  bond  through  sintering. 
Such  sintering  is  also  suggested  by  the  porosity  meas¬ 
urements.  Specifically,  the  matrix  porosity  diminished 
slightly  upon  exposure,  by  about  1%,  from  34.0  +  0.2% 
to  33.0+0.2%.  Similar  trends  in  both  strength  and  fib¬ 
er  pullout  have  been  reported  for  a  comparable  oxide 
CFCC  after  exposure  at  1200°C  in  air. 

Increases  in  the  strengths  of  both  the  interface  bond 
and  the  matrix  are  also  inferred  from  the  shear  tests  (Fig. 
9).  That  is,  the  shear  strength  increased  slightly,  by 
about  15%,  upon  exposure.  Nevertheless,  the  large  in¬ 
elastic  straining  capability  (  >  4%)  was  retained. 

The  hardness  measurements  pertinent  to  ring  #2 
are  summarized  in  Fig.  1 0.  The  values  for  both  the  ma¬ 
trix-rich  regions  and  the  fiber  tows  are  uniform  across 
the  plate  thickness.  More  importantly,  the  average  val¬ 
ues  appear  unchanged  following  exposure.  But,  since  the 
coefficients  of  variation  on  each  measurement  popula- 
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Fig.  8.  Fracture  surfaces  of  (a,  b)  pristine  (as-processed)  material  and  (c,  d)  after  exposure  in  the  combustor  (ring  #2). 


B 

B 


to 

a> 

c 

•o 

l— 

CO 

-C 

y> 

<5 

o 

5 


Normalized  position 

Fig.  10.  Variation  in  hardness  across  plate  thickness.  Inset 
summarizes  average  values  and  standard  deviations. 
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tion  are  about  20%,  relative  changes  in  hardness  of  less 
than  this  value  cannot  be  ascertained.  The  inference  is 
that  changes  in  both  matrix  and  interface  strength  are 
small:  probably  less  than  20%.  Indeed,  the  shear  meas¬ 
urements  suggest  a  change  of  only  about  1 5%,  (assum¬ 
ing  that  the  composite  shear  strength  is  proportional  to 
the  matrix  hardness).  It  should  be  noted  that  such 
changes  are  essentially  negligible  in  relation  to  the 
three-fold  change  in  hardness  associated  with  the  in¬ 
crease  in  precursor-derived  alumina  concentration  over 
the  narrow  range  of  1. 8-4.8%. 

Discussion  and  Concluding  Remarks 

Some  insights  into  the  mechanisms  of  matrix  sin¬ 
tering,  especially  within  the  fiber  tows,  have  been  ob¬ 
tained  from  independent  examinations  of  the  matrix 
microstructure  both  in  the  as-processed  condition  and 
after  exposure  at  high  temperature  in  air.18  Three  ob¬ 
servations  are  particularly  notable.  First,  in  the  mullite/ 
alumina  matrix  systems  used  here,  sintering  of  matrix 
particles  appears  to  be  restricted  to  the  narrow  matrix 
regions  within  fiber  tows.  In  contrast,  the  larger  ma¬ 
trix  regions  between  tows  appear  to  be  largely  unaffect¬ 
ed,  even  after  lOOOh  at  1200°C.  Second,  the  regions 
within  tows  show  higher  Al203/Si02  ratios  than  the 
bulk  matrix  (based  on  energy-dispersive  spectroscopy  in 
a  transmission  electron  microscope).  This  suggests  a 
higher  concentration  of  the  more-sinterable  alumina 
particles  within  the  tows  relative  to  that  in  the  bulk. 
Third,  the  regions  within  the  tows  tend  to  be  denuded 
of  the  larger  particles,  which  are  exclusively  mullite. 
This  observation  provides  additional  evidence  for  ele¬ 
vated  alumina  levels.  A  working  hypothesis,  consistent 
with  these  observations,  is  that  particle  segregation  oc¬ 
curs  during  slurry  infiltration.  Segregation  is  believed  to 
arise  from  two  phenomena.  The  first  one  is  the  physical 
exclusion  of  the  larger  mullite  particles  ( ~  1  pm)  from 
the  narrow  spaces  between  fibers  in  a  tow,  resulting  in  a 
higher  concentration  of  the  smaller  (~0.2  pm)  alumina 
particles.  The  second  is  the  “wall  effect,”  which  con¬ 
strains  the  packing  efficiency  of  larger  particles  around 
fiber  surfaces  and  creates  excess  void  space  that  is  pref¬ 
erentially  filled  by  the  finer  alumina  particles.19  In  com¬ 
bination,  these  phenomena  should  lead  to  increased 
sintering  within  the  tows  and  at  the  interfaces  with 
the  fibers,  consistent  with  the  measurements  and  obser¬ 
vations  presented  here. 


Two  strategies  for  improving  the  stability  of  the 
matrix  microstructure  against  such  sintering  are  envis¬ 
aged:  (i)  To  prevent  segregation,  the  sizes  of  the  two 
particle  types  should  be  more  closely  matched.  This 
could  be  accomplished  by  grinding  the  mullite  particles 
used  in  the  present  composites.  However,  this  size 
matching  has  undesirable  consequences  in  the  efficien¬ 
cy  of  particle  packing,  (ii)  The  alumina  panicles  within 
the  matrix  could  be  eliminated  altogether,  leaving  just 
the  mullite.  Bonding  of  the  particles  would  then  be  ef¬ 
fected  solely  by  the  precursor-derived  alumina.  Alumina 
in  the  latter  form  is  not  subject  to  this  type  of  segrega¬ 
tion  and  should,  in  principle,  yield  a  more  uniform 
microstructure.  Composites  of  this  type  have  been  suc¬ 
cessfully  fabricated,  although  an  assessment  of  their  per¬ 
formance  and  examinations  of  their  microstructure  have 
yet  to  be  performed. 

The  presence  of  delamination  cracks  in  the  exposed 
rings  highlights  the  long-known  problem  of  low  trans¬ 
verse  toughness  of  ceramic  composites  with  2D  fiber  ar¬ 
chitectures.  The  problem  is  exacerbated  in  the  porous 
systems  relative  to  those  with  fully  dense  matrices.  Sig¬ 
nificant  improvements  in  delamination  resistance  of  such 
systems  can  only  be  achieved  through  the  use  of  through 
thickness  reinforcements,  introduced,  for  example,  by 
stitching  2D  fiber  architectures,  or  by  using  3D  weaves. 

Notwithstanding  the  deficiencies,  the  present  com¬ 
posite  system  has  performed  well  in  relation  to  other 
porous  oxide  composites,  especially  those  with  alumino¬ 
silicate  matrices.  In  the  latter  systems,  combustor  rig 
tests  of  the  type  performed  here  have  led  to  the  devel¬ 
opment  of  macroscopic  transverse  cracks  in  which  the 
fibers  are  broken,  even  after  only  a  few  hours  of  expo¬ 
sure.20  This  suggests  rapid  sintering  of  the  matrix  and 
an  attendant  loss  in  in-plane  damage  tolerance.  For  rea¬ 
sons  not  yet  understood,  this  sintering  appears  to  be 
accelerated  in  the  presence  of  water  vapor,  such  as  that 
in  the  combustion  environment.  This  presents  new 
challenges  in  the  design  and  synthesis  of  materials  sys¬ 
tems  that  retain  their  damage  tolerance  characteristics 
over  extended  service  periods. 
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The  effects  of  matrix  strength  on  the  notched  and  unnotched 
tensile  properties  of  a  family  of  porous-matrix  oxide  composites 
are  examined  both  experimentally  and  theoretically.  Experi¬ 
ments  are  performed  on  three  composites,  distinguished  from 
one  another  by  the  amount  of  binding  alumina  within  the  matrix. 
Increases  in  alumina  concentration  produce  elevations  in  un¬ 
notched  tensile  and  shear  strengths,  but  the  benefits  are  offset  by 
an  increase  in  notch  sensitivity.  The  degree  of  notch  sensitivity  is 
rationalized  on  the  basis  of  a  model  that  accounts  for  interac¬ 
tions  between  notch  tip  tensile  and  shear  bands.  The  model  pre¬ 
dictions  are  cast  in  terms  of  the  ratio  of  the  notch  length  to  a 
characteristic  bridging  length  scale.  These  results,  in  turn,  form 
the  basis  for  a  simple  analytical  formula  for  notched  strength, 
accounting  for  effects  of  elastic  anisotropy  and  finite  sample 
size.  The  utility  of  this  formula  in  predicting  notched  strength  is 
assessed.  Issues  associated  with  bridging  law  shapes  and  bridg¬ 
ing  length  scales  are  addressed.  The  effect  of  alumina  concen¬ 
tration  on  notch  sensitivity  is  discussed  in  terms  of  its  influence 
on  the  bridging  length  scale,  dictated  by  the  interplay  between 
the  unnotched  tensile  strength,  the  longitudinal  Young’s  modu¬ 
lus,  the  degree  of  in-plane  elastic  anisotropy,  and  the  fracture 
energy.  The  net  result  is  a  decreasing  bridging  length  scale  and 
hence  increasing  notch  sensitivity  as  the  matrix  is  strengthened 
with  alumina. 


I.  Introduction 

Fiber-reinforced  oxide  composites  are  being  targeted  for  use 
in  future  gas  turbine  engines.  Their  obvious  advantage  over 
monolithic  ceramics  is  their  superior  damage  tolerance  and  low 
notch  sensitivity.  Notwithstanding  progress  in  the  development 
of  fiber  coatings  and  environmental  barrier  coatings  for  use  with 
SiC/SiC  composites,  oxide  systems  remain  top  contenders  for 
applications  requiring  long  service  lives  (>  104  h)  in  oxidizing 
environments,  albeit  at  operating  temperatures  (about  1200°C) 
that  are  some  100-200  K  below  the  targeted  levels  for  SiC/SiC. 
This  has  motivated  numerous  activities  in  material  development 
and  assessment  in  industrial1-5  and  government  laboratories,6-12 
as  well  as  academia.13-19 

The  present  article  focuses  on  tensile  strength  and  notch  sen¬ 
sitivity  of  one  particular  class  of  oxide  composite:  that  in  which 
damage  tolerance  is  enabled  by  a  porous  matrix,  rather  than 
through  the  conventional  route  involving  fiber  coatings.  It 
builds  on  recent  experience  in  composite  fabrication  and  micro- 
structural  control,  through  appropriate  selection  of  matrix  con¬ 
stituents  and  their  topologies.15,20  The  emphasis  on  notch 
sensitivity  stems  from  the  targeted  applications  where  holes 
and  other  stress  concentrators  are  expected  to  be  present,  for 
cooling  or  attachment. 
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In  notched  laminates  with  2D  fiber  architectures  in  a  0790' 
orientation,  failure  occurs  through  two  concurrent  damage 
processes  (Fig.  1).  The  first  involves  shear  deformation  ahead 
of  the  notch,  along  bands  aligned  parallel  to  the  loading  direc¬ 
tion.1216,21  This  deformation  occurs  as  a  result  of  the  shear 
stress  acting  parallel  to  the  fiber  axes  coupled  with  the  inherent 
anisotropy  in  the  inelastic  straining  characteristics  of  2D  lam¬ 
inates.  The  second  involves  a  tensile  process  zone,  oriented  per¬ 
pendicular  to  the  applied  stress.  Here,  the  dominant  damage 
mechanisms  are  matrix  cracking,  interfacial  debonding  and 
sliding,  and  fiber  fracture  and  pull  out.  These  ultimately  lead 
to  the  formation  and  propagation  of  a  macroscopic  crack.  The 
interplay  between  the  shear  bands  and  the  tensile  process  zone 
dictates  the  degree  of  notch  sensitivity.  That  is,  the  shear  de¬ 
formation  reduces  the  crack  tip  stresses,  thereby  increasing  the 
applied  stress  needed  for  fiber  fracture  and  composite  rupture. 

The  subsequent  study  is  predicated  on  the  notion  that  the 
state  of  the  matrix  should  influence  the  notched  properties 
in  two  ways,  (i)  For  uncoated  fiber  systems,  the  key  interface 
properties  (the  debond  toughness  and  the  sliding  resistance)  are 
expected  to  move  in  tandem  with  the  matrix  mechanical 
properties:  increasing  with  decreasing  porosity.  In  the  domain 
of  global  load  sharing,  wherein  debonding  and  sliding  occur 
readily  at  the  fiber-matrix  interface,  elevations  in  interface  prop¬ 
erties  increase  the  fiber  bundle  strength.22  At  the  other  extreme, 
where  the  interface  is  strong  and  hence  strongly  local-load  shar¬ 
ing  conditions  are  obtained,  the  bundle  strength  diminishes  with 
increasing  matrix  strength  and  becomes  inherently  stochastic.23 
This  implies  an  optimum  matrix  condition  at  which  the  un¬ 
notched  composite  strength  achieves  its  maximum  value,  (ii)  El¬ 
evations  in  matrix  strength  increase  the  shear  strength  of  the 
composite  and  hence  reduce  the  spatial  extent  of  the  shear  bands 
in  the  notch  tip  region.  As  a  result,  the  notched  strength  is  re¬ 
duced  relative  to  the  unnotched  value.  One  of  the  objectives  of 
the  present  study  is  to  probe  experimentally  the  interplay  be¬ 
tween  the  shear  strength  and  the  notched  and  unnotched  tensile 
strengths  of  the  composite  as  the  state  of  the  matrix  is  varied, 
from  very  weak  to  strong.  A  second  objective  is  to  assess  the 
utility  of  an  existing  model  for  notch  sensitivity24  in  describing 
the  experimental  results. 

The  article  is  organized  in  the  following  manner.  A  model  for 
the  notched  strength  of  porous  matrix  composites  is  presented 
first,  with  a  view  to  identifying  the  key  governing  material  prop¬ 
erties.  This  is  followed  by  a  description  of  the  procedures  used  in 
composite  fabrication  and  characterization  as  well  as  mechan¬ 
ical  testing.  The  experiments  are  guided  largely  by  the  predic¬ 
tions  of  the  model.  An  assessment  of  the  model  is  then  made 
through  appropriate  comparisons  between  the  measurements 
and  the  predictions.  The  latter  includes  a  proposed  method  of 
partitioning  the  fracture  energy  into  contributions  coming  from 
two  domains:  before  and  after  fiber  bundle  failure.  Because  of 
differences  in  the  characteristic  bridging  length  scales  of  the  two 
bridging  domains,  their  effects  on  notch  sensitivity  vary  with 
notch  length. 


II.  A  Model  of  Notched  Strength 

Numerical  simulations  probing  the  combined  effects  of  shear 
and  tensile  bands  on  notch  sensitivity  of  strength  have  been 
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Fig.  1.  Notch  tip  damage  processes  in  porous  matrix  ceramic  composites. 

performed  by  He  et  a!.24  Here,  the  tensile  band  is  treated  as  a 
cohesive  zone  and  is  represented  mechanically  by  a  continuous 
array  of  linear  springs  with  a  strength,  a0,  equal  to  the  un¬ 
notched  composite  strength.  Analogously,  the  shear  band  is 
treated  as  a  rigid-perfectly  plastic  cohesive  zone,  with  a  shear 
strength,  t0.  Two  types  of  simulations  have  been  performed.  In 
the  first,  an  infinite,  elastically  isotropic  specimen  with  a  semi¬ 
infinite  notch  is  subjected  to  a  remote  displacement  field  pre¬ 
scribed  by  the  square  root  singular  solution  and  scaled  with  the 
stress  intensity  factor,  K.  The  fracture  toughness,  Kc,  is  obtained 
when  the  maximum  stress  within  the  tensile  band  reaches  its 
strength,  cr0.  The  corresponding  fracture  energy,  f,  is  calculated 
from  K:  using  the  Irwin  relation.  The  results  of  these  simulations 
are  plotted  in  Fig.  2(a).  When  xo/ao<0.2,  the  fracture  energy 
increases  rapidly  with  decreasing  t0/ct0,  in  approximate  accord¬ 
ance  with  the  power  law  indicated  in  the  figure.  In  contrast, 
when  To/oo  >  0.2,  the  fracture  energy  rapidly  approaches  the 
value  intrinsic  to  the  tensile  band,  r0;  the  beneficial  effect  of  the 
shear  band  on  the  fracture  energy  is  lost. 

In  the  second,  a  semi-infinite,  elastically  isotropic  specimen 
with  a  finite  notch  of  length  a()  is  subjected  to  a  uniform  remote 
tension.  The  objective  is  to  ascertain  the  degree  of  notch  sensi¬ 
tivity  of  strength  as  the  notch  length  and  the  material  properties 
are  vaned._When  plotted  against  the  normalized  notch  length, 
aoOo/£T  ( E  being  the  plane  strain  modulus),  the  results  for  the 
normalized  notch  strength  fall  essentially  onto  a  single,  mono- 
tonically  decreasing  band,  only  weakly  dependent  on  T0/a0-  Ev¬ 
idently,  the  beneficial  effects  of  the  shear  band  are  embodied  in 
the  fracture  energy,  T,  in  the  sense  that  reductions  in  t0/ct0  lead 
to  an  increase  in  T  and  hence  a  reduction  in  notch  sensitivity. 
Furthermore,  the  trend  can  be  well  described  by  the  empirical 
formula 


This  formula  is  obtained  essentially  by  interpolating  between 
the  strengths  for  the  two  limiting  cases  of  very  short  and  very 
long  notches.  Jn  these  limits,  Eq.  (1)  yields:  (i)  oN/a0 
— *  las  jiCqUo/ET— >0  (the  notch  insensitive  domain)  and 
(ii)  ohJgo  — *  {Er/nolao)^2  (the  Griffith  stress)  when 
nolao/Er  »  1. 

The  same  formula  has  found  use  in  describing  notched 
strength  in  cases  where  fracture  occurs  exclusively  through  a 
tensile  band  (the  shear  band  being  absent).  Here,  again,  Eq.  (1) 
provides  an  accurate  representation  of  numerical  (finite  element) 
results,  within  about  10%,  for  a  variety  of  simple  monotonic 


(a) 


Fig.  2.  Effects  of  strength  ratio,  t0/cto,  on:  (a)  normalized  composite 
toughness  (obtained  in  an  infinite  notched  specimen)  and  (b)  strength  of 
an  edge-notched  semi-infinite  specimen  (adapted  from  He  el  at24).  Also 
shown  in  (a)  arc  the  pertinent  strength  ratios  of  the  materials  tested  in 
this  study,  designated  by  the  concentration  of  precursor-derived  alumina 
in  the  matrix. 


traction  laws,  including  linear  and  parabolic  hardening,  linear 
softening,  and  rectilinear.25  For  more  complex  traction  laws, 
shape  plays  an  important  role.  For  instance,  when  bridging  is 
characterized  by  a  two-level  rectilinear  law,  there  are  two  bridg¬ 
ing  stress  levels  and  two  fracture  energies:  features  not  embodied 
in  Eq.  (1).  The  ramifications  of  this  deficiency  are  addressed 
further  in  Section  VI. 

Inspection  of  Eq.  (1)  reveals  a  characteristic  length  scale  that 
dictates  notch  sensitivity,  defined  by 

£T 

tfch  =  —  (2) 

co 

With  this  definition,  the  notched  strength  formula  can  be  written 
more  compactly  as 


In  this  form,  the  role  of  och  in  the  notch  sensitivity  can  be  as¬ 
certained.  Two  notable  results  emerge,  (i)  When  a0  =  ach,  the 
predicted  strength  ratio  is  very  close  to  0.5:  a  rather  significant 
reduction,  relative  to  unity.  Only  when  a0  <C  ach  is  the  notch 
sensitivity  small,  i.e.  ctn/ct0->  1 .  In  this  context,  ach  represents 
the  pertinent  length  against  which  notch  or  flaw  lengths  must 
be  compared  in  assessing  the  degree  of  notch  sensitivity,  (ii)  A 
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complementary  view  of  ach  is  obtained  by  comparing  the  pre¬ 
dictions  of  Eq.  (3)  with  the  Griffith  stress.  For  instance,  when 
aQ  =  ach,  the  normalized  Griffith  stress  (assuming  all  of  the  frac¬ 
ture  energy  is  fully  utilized  in  notched  strength),  is  about  0.57: 
only  12%  higher  than  the  value  predicted  by  Eq.  (3).  Thereafter, 
for  a0>ach,  the  two  predictions  converge,  allowing  use  of  the 
Griffith  stress  in  failure  prediction.  In  this  context,  ac h  can  be 
viewed  as  the  critical  notch  length  beyond  which  Griffith-like 
behavior  is  obtained,  with  the  associated  inverse  square  root 
sensitivity  of  strength  to  notch  length. 

To  compare  the  preceding  predictions  with  experiment,  Eq. 
(1)  is  modified  to  account  for  both  finite  specimen_dimensions 
and  elastic  anisotropy.  In  the  limit  where  nalao/Er  —>  0,  the 
normalized  notched  strength  (expressed  on  a  net-section  basis) 
approaches  the  sanre  value,  notably  aN/a0  -»  I .  In  the  Griffith 
limit,  where  rta^i0/£r  >  1,  the  energy  release  rate  is  modified 
to  read  as 


With  the  modified  definition,  the  notched  strength  formula  in 
Eq.  (9)  reduces  to  the  one  in  Eq.  (3). 

An  assessment  of  Eq.  (9)  is  made  through  comparison  with 
the  experiments  described  in  the  following  sections.  The  model 
identifies  five  fundamental  material  properties  that  govern 
notched  strength:  the  elastic  constants  £),  v,2,  and  pi ,  2,  the  un¬ 
notched  tensile  strength,  a0,  and  the  fracture  energy,  T.  Addi¬ 
tional  insights  into  the  contributions  from  shear  deformation  to 
notch  sensitivity  are  obtained  through  measurement  of  the  shear 
strength,  r0»  >n  combination  with  the  results  plotted  in  Fig.  2(a). 
The  elastic  constants  and  the  strengths  are  obtained  from  stand¬ 
ard  instrumented  tension  and  shear  tests,  whereas  the  fracture 
energy  is  estimated  from  the  work  of  fracture  in  notched  tension 
tests.  The  latter  tests  are  performed  on  specimens  with  varying 
notch  lengths,  for  the  purpose  of  establishing  the  dependence  of 
strength  on  size. 


G  =  ^g(p,X)F2(9;|)(l -aa/Wf  (4) 

where  W  is  the  specimen  width,  H  is  the  specimen  length,  F[ao/ 
W,H/W)  is  a  non-dimensional  function  that  depends  on  both 
the  specimen  geometry  and  the  loading  configuration,26  g(p,  X) 
characterizes  the  effect  of  orthotropy, 27-28  and  p  and  X  are 
defined  by 


with  the  subscripts  1  and  2  referring  to  the  principal  material 
axes,  p  the  shear  modulus,  and  v  Poisson’s  ratio.  Composites 
with  balanced  symmetric  fiber  lay-ups  exhibit  cubic  in-plane 
symmetry,  whereupon  X  =  EfE2  =  1  and  v12  =  v2i,  and  g( p,  X)  is 
given  by 


=  (-^-£)  '  (1  +0.1(1  -p)  -0.016(1  -p)2. 
+0.002(1  -  p)3) 


(7) 


Upon  setting  G  =  T,  the  notched  strength  becomes 
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Interpolating  between  the  two  limits,  in  a  manner  analogous  to 
that  used  to  obtain  Eq.  (1),  gives  the  notched  strength  as 
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In  light  of  this  result,  the  definition  of  the  characteristic  bridging 
length  ach  is  modified  to  take  into  account  both  the  elastic  an¬ 
isotropy  as  well  as  the  finite  specimen  width: 


_ £T _ 

S( P-X) 


(10) 


III.  Materials 

The  composites  were  comprised  of  seven  layers  of  Nextel  720 
fiber  cloth  in  an  8-hamess  satin  weave  and  a  porous  matrix  of 
mullite  and  alumina.  The  processing  route  was  essentially  iden¬ 
tical  to  that  described  in  earlier  articles.1718  Briefly,  the  matrix 
was  introduced  in  two  steps:  (i)  vacuum  infiltration  of  a  slurry 
containing  80%  mullite  and  20%  alumina  particles  into  the  fiber 
preform  and  (ii)  impregnation  and  pyrolysis  of  an  alumina  pre¬ 
cursor  solution  (Al2Cl(OH)5).  The  effective  volumetric  yield  of 
the  precursor  solution  upon  conversion  to  y-alumina,  deter¬ 
mined  in  the  manner  described  below,  was  6.4%.  To  prevent 
segregation  of  the  precursor  during  drying,  the  panels  were 
gelled  by  placing  them  in  a  chamber  containing  gaseous 
NH4OH  for  6  h  immediately  following  impregnation.  After  dry¬ 
ing,  the  panels  were  fired  at  900°C  for  2  h,  thereby  converting 
the  precursor  to  y-alumina.  After  all  impregnation  cycles  were 
completed,  the  panels  were  given  a  final  heat  treatment  at 
1200°C  for  2  h.  In  the  latter  step,  the  y-alumina  is  converted 
to  a-alumina.  The  average  fiber  volume  fraction  was  38%. 

Three  different  composites  were  fabricated,  distinguished 
from  one  another  by  the  number  of  precursor  impregnation 
and  pyrolysis  cycles.  Selection  of  the  pertinent  range  was  based 
in  part  on  the  results  of  a  separate  study,  focused  on  the  effects 
of  the  volume  fraction  of  precursor-derived  alumina,  KPDA,  on 
matrix  properties.29  The  results  of  the  latter  study,  coupled  with 
the  He  and  Hutchinson30  crack  deflection  criterion,  suggest  that 
the  deflection  condition  is  not  satisfied  once  KPDa  exceeds  a 
critical  value  of  about  9%— 10%.  With  the  current  precursor, 
this  critical  value  is  achieved  after  five  impregnation  and  pyr¬ 
olysis  cycles.  Consequently,  to  explore  the  effects  of  matrix 
strength,  all  within  the  crack  deflection  domain,  panels  were 
fabricated  with  either  zero,  two,  or  four  impregnation  cycles. 

Matrix  porosity  and  concentration  of  PDA  were  measured 
using  the  Archimedes  method.31  The  initial  matrix  porosity  (fol¬ 
lowing  slurry  infiltration  but  prior  to  precursor  impregnation) 
was  p°m  =  39%.  Assuming  that  the  pore  space  remains  open 
during  processing  and  the  precursor  solution  completely  fills  the 
pore  volume  during  each  impregnation,  the  matrix  porosity  af¬ 
ter  N  impregnation  cycles  is  predicted  to  follow  the  relation 

pm=P°m(l-(l-(l-3'P)'V)^)  (ID 

where  pm  is  the  matrix  porosity,  yp  is  the  precursor  yield  (upon 
conversion  to  y-alumina),  and  p  is  the  density  of  the  PDA  in  the 
a  or  y  states.  The  precursor  yield  (yp  =  6.4%)  was  calculated 
from  Eq.  (1 1)  using  the  measured  values  of  pm  and  p°m  for  both 
N  =2  and  4.  The  change  in  porosity  following  precursor  im¬ 
pregnation  was  also  used  to  determine  the  KPDA,  as  summarized 
in  Table  I.  The  latter  spans  the  targeted  range,  from  0%  to  8%. 

The  effect  of  PDA  concentration  on  matrix  strength  was 
characterized  through  hardness  measurements.  For  this  pur¬ 
pose,  Vickers  indentations  were  made  on  polished  cross-sections 
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Table  I.  Summary  of  Elastic  Constants  and  Anisotropy 
Parameters  as  a  Function  of  Precursor  Concentration 


N 

IW  (%) 

E,  (GPa) 

VJ2 

lt,2  (GPa) 

p 

s(p) 

0 

0 

48.3 

0.06 

4.3 

5.5 

3.1 

2 

3.9 

56.1 

0.10 

8.4 

3.3 

2.0 

4 

8.2 

63.0 

0.13 

14.0 

2.2 

1.5 

PDA,  precursor-derived  alumina. 


Fig.  3.  Variation  in  matrix  hardness  with  precursor-derived  alumina 
(PDA)  concentration. 


within  matrix-rich  regions  between  the  fiber  cloths,  with  a  mass 
of  200  g.  The  results  from  such  measurements,  plotted  in  Fig.  3, 
reveal  a  strong  sensitivity  of  matrix  hardness,  Hm,  to  Fpdai  in¬ 
creasing  by  more  than  an  order  of  magnitude  over  the  range 
V pda  =  0%  8%.  Moreover,  they  demonstrate  that  the  alumina 
derived  from  the  precursor  has  a  much  stronger  influence  on  the 
matrix  than  that  associated  with  the  particulates  in  the  slurry. 
For  instance,  the  addition  of  less  than  1  %  PDA  has  the  effect  of 
doubling  the  matrix  hardness,  relative  to  the  unimpregnated 
material  containing  20%  alumina  in  particulate  form. 

Hereafter,  experimental  results  for  various  composite  prop¬ 
erties  are  presented  in  terms  of  either  the  PDA  concentration  or 
the  matrix  hardness,  as  appropriate.  Their  connections  to  the 
number  of  precursor  impregnation  cycles  can  be  established 
through  the  results  plotted  in  Fig.  3  and  summarized  in  Table  1. 


IV.  Mechanical  Testing  Procedures 

A  variety  of  mechanical  tests  were  performed  on  each  of  the 
three  composite  panels,  with  the  goal  of  assessing  the  degree  of 
notch  sensitivity  of  strength  in  the  0°/90°  orientation.  To  this 
end,  measurements  were  made  of  the  key  properties  described 
above,  notably:  the  elastic  constants  £),  vi2,  and  pi2,  the  un¬ 
notched  tensile  strength,  a0,  the  shear  strength,  t0,  and  the  frac¬ 
ture  energy,  T.  A  secondary  objective  was  to  assess  the  effects  of 
matrix  strength  on  the  pertinent  composite  properties. 

The  unnotched  tensile  strengths  were  measured  using  uniaxial 
tensile  bars,  approximately  10  mm  wide,  and  with  gauge  lengths, 
L,  of  either  20  or  64  mm.  (The  rationale  for  probing  gauge 
length  effects  stems  from  the  length  of  the  fracture  process  zone 
in  the  notched  specimens,  as  described  below.)  The  displacement 
rate  was  0.012  mm/min.  In  all  cases,  the  longitudinal  strains 
were  measured  both  through  a  non-contacting  laser  extensome- 
ter  and  a  strain  gauge.  Transverse  strains  were  also  measured 
with  a  strain  gauge. 

To  assess  the  effects  of  the  matrix  on  the  0790°  strength, 
the  tensile  properties  of  dry  fiber  cloth  were  also  measured. 


Gauge  length,  L  (mm) 

Fig.  4.  Effects  of  gauge  length  and  matrix  hardness  on  unnotchcd  ten¬ 
sile  strength. 

using  test  coupons  with  a  gauge  length  of  64  mm  (identical 
to  the  length  of  one  set  of  composite  specimens)  and  a  gauge 
width  of  10  mm.  Test  samples  were  cut  from  the  same  fiber 
cloth  used  in  composite  fabrication  and  fired  at  550°C  for  2  h 
to  remove  the  factory-applied  sizing.  The  sample  ends  were 
tabbed  with  silicone  rubber  and  clamped  between  frictional 
grips  for  testing. 

Notch  sensitivity  of  strength  was  examined  through  single¬ 
edge  notched  tension  tests  (SENT).*  The  ratio  of  notch  length  to 
specimen  width  was  fixed  at  0.5  and  scale  effects  were  probed  by 
varying  the  specimen  width,  from  12  to  24  mm.  The  height  to 
width  ratio  was  fixed  at  2.5.  Crack  mouth  opening  displace¬ 
ments  were  measured  using  a  spring-loaded  gauge  mounted 
between  knife  blades  affixed  to  the  sample.  The  remote 
displacement  was  also  monitored  using  an  LVDT.  Additional¬ 
ly,  the  sequence  of  damage  events  leading  to  fracture  was  ob¬ 
tained  by  optical  imaging  of  the  external  specimen  surface.  To 
enhance  the  contrast  associated  with  such  damage,  the  surfaces 
were  coated  with  a  thin  (30  nm)  layer  of  Au-Pd. 

The  in-plane  shear  response  was  measured  using  the  Iosipescu 
test  in  the  0°/90°  orientation.32  Specimens  were  19  mm  wide  with 
symmetrical  V-notches,  each  3.8  mm  deep,  with  an  included 
angle  of  105°.  Gauges  were  mounted  at  +45°  to  measure  shear 
strains.  The  tests  were  performed  at  a  displacement  rate  of 
0.05  mm/min. 


V.  Experimental  Measurements  and  Observations 
(1)  Unnotched  Tensile  Properties 

Figure  4  demonstrates  the  effects  of  gauge  length  and  PDA 
concentration  on  the  unnotched  tensile  strength.  The  unimpreg¬ 
nated  composite  exhibits  a  significant  sensitivity  of  strength  to 
gauge  length,  diminishing  by  about  1/3,  from  120  to  80  MPa, 
over  the  range  L  =  20-64  mm.  In  contrast,  the  composites  with 
either  4%  or  8%  PDA  show  no  apparent  length  dependence. 
Furthermore,  the  tensile  strength  increases  with  Vpoa  for  con¬ 
stant  gauge  length.  For  L  =  20  mm,  the  strength  improvement  is 
a  modest  10%-15%  over  the  range  KpDA  =  0%-8%;  for  L  -  64 
mm,  the  increase  is  more  substantial,  approaching  60%.  Chang¬ 
es  in  E\  and  v,2  are  also  observed  over  this  range  of  PDA  con¬ 
centrations  (Table  1).  Modulus  rises  30%,  from  48  to  63  GPa. 
while  Poisson’s  ratio  increases  nearly  twofold,  from  0.06  to  0.13. 

Also  shown  for  comparison  in  Fig.  4  is  the  strength  of  the 
fiber  cloth,  taken  to  be  the  product  of  the  fiber  volume  fraction 

!Wiih  damped  ends,  SENT  provides  a  more  critical  assessment  of  notch  sensitivity  than 
double  edge  notched  tension  (DENT),  because  of  the  higher  stress  intensity  at  a  prescribed 
normalized  notch  length,  a^W.  A  secondary  advantage  of  SENT  is  that  it  requires  signif¬ 
icantly  less  materia)  for  testing:  both  the  specimen  width  and  the  specimen  length  arc  one 
half  of  the  values  needed  for  DENT  at  a  fixed  oq/  W. 
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Fig.  5.  Macrophotographs  of  fractured  tensile  specimens,  showing  re¬ 
duction  in  degree  of  fiber  pullout  with  increasing  precursor-derived  alu¬ 
mina  (PDA)  concentration. 


(38%)  and  the  mean  cloth  strength  (186  +  33  MPa).  The  result¬ 
ing  value  (71  ±13  MPa)  is  in  reasonable  agreement  with  the 
mean  strength  (about  80  MPa)  of  the  composite  with  no  PDA. 
This  correlation  indicates  that  the  composite  strength  is  dictated 
entirely  by  the  fiber  cloth,  with  no  beneficial  effect  of  the  matrix. 
Moreover,  it  suggests  that  the  gauge  length  dependence  of  the 
strength  of  this  composite  can  be  rationalized  on  the  basis  of  the 
corresponding  length  dependence  of  the  cloth  itself.  Assuming  a 
Weibull  distribution  for  the  fiber  strength,  the  predicted  length 
dependence  is  given  by 


where  a0  and  Lo  arc  the  reference  strength  and  reference  length, 
respectively.  From  the  slope  of  the  line  in  Fig.  4  (presented  in 
log-log  coordinates),  m  is  estimated  to  be  about  2.5  + 0.4.  (Com¬ 
parisons  of  m  with  previously  reported  values  on  pristine  fibers 
is  problematic,  since  the  latter  values  span  the  range  of  2.8-26, 
depending  on  the  method  33~3S  Further,  these  values  pertain 
only  to  pristine  fibers  and  may  not  accurately  represent  the  fiber 
condition  after  weaving.) 

Accompanying  the  elevations  in  strength  with  FPDA  is  a 
change  in  the  appearance  of  the  fracture  surfaces  (Fig.  5).  In 
the  unimpregnated  composite,  fiber  pullout  is  extensive,  with 
long  exposed  fiber  lengths.  As  the  matrix  strength  is  increased 
with  PDA  addition,  a  significant  reduction  in  pullout  lengths  is 
apparent.  This  observation  is  qualitatively  consistent  with  mod¬ 
els  of  fiber  fragmentation  and  bundle  failure,  where  pull¬ 
out  lengths  are  predicted  to  fall  with  elevations  in  interface 
properties.22 

(2)  Notched  Tensile  Properties 

The  notched  tensile  response  and  associated  damage  mecha¬ 
nisms  are  shown  in  Figs.  6-8.  For  samples  with  FPDA  =  0%,  the 
first  apparent  damage  event  is  a  shear  band  at  the  notch  tip, 
located  at  the  boundary  between  fiber  tows  (Fig.  6(a)),  and  is 
manifested  in  the  onset  of  non-linearity  in  the  stress-displace¬ 
ment  response  (Fig.  7).  As  the  stress  increases,  additional  shear 
bands  form  along  tow  boundaries,  with  small  intervening  tensile 
cracks  (the  latter  only  barely  evident  on  Fig.  6(a)).  The  shear 
bands  extend  to  a  distance  of  about  10  mm  on  either  side  of  the 
notch  plane.  The  sequence  of  shear  banding  and  tensile  cracking 
proceeds  until  all  fiber  bundles  have  fractured.  The  subsequent 
response  exhibits  significant  softening  as  the  broken  fibers  pull 
out,  presumably  resisted  by  friction  as  they  slide  past  the  sur¬ 
rounding  matrix.  The  spatial  extent  of  the  damage  and  the 
pullout  is  evident  in  the  low-magnification  image  of  a  typical 
fractured  specimen,  shown  in  Fig.  8(a). 


For  composites  with  the  intermediate  matrix  (  FPDA  =  4%),  a 
similar  sequence  of  damage  events  accompanies  fracture.  How¬ 
ever,  in  this  case,  the  length  of  the  shear  bands  and  the  extent  of 
fiber  pullout  are  reduced  significantly,  as  seen  in  the  sequence 
shown  in  Fig.  6(b)  as  well  as  the  fractured  specimen  in  Fig.  8(b). 
Additionally,  the  peak  stress  is  lower  and  occurs  at  a  lower  value 
of  remote  displacement,  the  latter  being  consistent  with  the  re¬ 
duced  shear  band  length. 

For  the  composites  with  the  strongest  matrix  ( KPDA  =  8%), 
shear  banding  is  almost  negligible.  Instead,  fracture  occurs  es¬ 
sentially  by  the  propagation  of  a  single  tensile  crack,  with  some 
meandering  from  a  pure  mode  1  path.  This  view  is  consistent 
with  the  rather  narrow  band  of  damage  and  short  pullout 
lengths  seen  on  the  fracture  surface  (Fig.  8(c)).  Accompanying 
this  diminished  shear  banding  is  a  further  reduction  in  the  peak 
stress  and  the  corresponding  displacement. 

Figure  9  summarizes  the  effects  of  notch  length  and  precursor 
addition  on  the  notched  strength.  In  all  materials,  some  notch 
sensitivity  is  observed,  with  strength  diminishing  with  increasing 
notch  length.  Further,  the  degree  of  notch  sensitivity  increases 
with  PDA  concentration.  Over  the  range  ao  =  0-12  mm,  the 
unimpregnated  composite  strength  falls  approximately  25%.  In 
contrast,  the  densest  composite  ( FPDA  =  8%)  exhibits  a  60% 
strength  reduction  over  the  same  notch  length  range.  These 
trends  suggest  strong  dependence  of  notch  sensitivity  on  the 
spatial  extent  of  the  shear  bands,  in  qualitative  agreement  with 
the  theoretical  results  presented  in  Section  II.  A  quantitative 
assessment  is  made  in  Section  VI. 

(3)  Shear  Properties 

Changes  in  shear  response  with  PDA  concentration  are  shown 
in  Fig.  10.  In  each  case,  the  response  is  essentially  elastic-plastic, 
with  only  modest  strain  hardening.  Failure  occurs  through  the 
formation  of  a  localized  damage  band  across  the  specimen  sec¬ 
tion,  leading  to  a  plateau  stress  at  a  strain  of  2%-3%:  well  be¬ 
yond  the  tensile  failure  strain  of  the  fibers.  This  behavior  is 
reminiscent  of  the  shear  response  of  fiber  composites  with  duc¬ 
tile  matrices,  including  both  polymers  and  metals.  Further,  there 
is  a  pronounced  sensitivity  of  the  shear  modulus  and  strength  to 
the  state  of  the  matrix.  Roughly  fourfold  increases  in  both  prop¬ 
erties  are  obtained  as  the  PDA  concentration  is  increased  from 
0%  to  8%,  with  the  shear  modulus  rising  from  4  to  14  GPa 
(Table  I)  and  the  strength  increasing  from  8  to  32  MPa. 

These  elastic-plastic  characteristics  suggest  that  the  compos¬ 
ite  shear  strength  should  scale  with  the  shear  “flow  stress”  of  the 
matrix,  as  obtained  in  metal  and  polymer  matrix  composites. 
Using  the  Tresca  yield  criterion  (applicable  to  materials  that 
deform  by  dislocation  plasticity),  the  matrix  shear  strength  is 
predicted  to  be  half  the  uniaxial  yield  stress.  Furthermore,  the 
uniaxial  yield  stress  is  roughly  three  times  the  matrix  hardness, 
Hm.  Upon  combining  these  results,  the  predicted  composite 
shear  strength  is  about  Hm/ 6. 

Figure  10(b)  shows  the  dependence  of  the  composite  shear 
strength  on  matrix  hardness.  For  the  weakest  matrix,  the  com¬ 
posite  strength  is  consistent  with  the  rudimentary  prediction: 
Hm/6.  For  stronger  matrices,  the  shear  strength  rises  with  matrix 
hardness,  but  at  a  much  lower  rate  than  predicted.  The  discrep¬ 
ancy  is  likely  because  of  the  assumed  form  of  the  yield  criterion. 
Specifically,  the  porous  ceramic  matrix  likely  follows  a  different 
yield  condition  (not  the  Tresca  criterion).  The  discrepancy  high¬ 
lights  the  need  for  the  development  of  appropriate  flow  laws  for 
the  porous  matrices,  for  use  in  predicting  composite  properties. 


VI.  Analysis  of  Notch  Sensitivity 

The  predicted  notched  strengths  from  the  model  presented  in 
Section  II  have  been  calculated  using  Eq.  (9)  in  combination 
with  the  measured  elastic  constants  E t,  p)2,  and  G]2,  the  un¬ 
notched  tensile  strength,  a0,  and  the  work  of  fracture,  T.  The 
elastic  parameters  and  corresponding  anisotropy  function  g(p) 
are  summarized  in  Table  1.  Because  of  the  absence  of  gauge 
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Fig.  6.  Progression  of  damage  in  edge-notched  tensile  specimens:  (a)  VrDA  =  0,  (b)  KPDA  =  4,  (c)  KPDA  =  8.  In  each  sequence,  the  first  picture  represents 
the  pristine  (unloaded)  specimen.  The  subsequent  images  are  at  progressively  increasing  displacements,  indicated  in  Fig.  7.  PDA,  precursor-derived 
alumina. 


length  dependence  in  materials  with  either  4%  or  8%  PDA,  o0 
was  taken  to  be  the  average  of  all  measured  unnotched  tensile 
strengths  (including  both  20  and  64  mm  gauge  lengths).  In  con¬ 
trast,  for  the  unimpregnated  material,  only  the  values  measured 
for  the  20  mm  gauge  length  were  used,  since  this  value  corre¬ 
sponds  more  closely  to  the  length  of  the  damage  zone  at  the 
notch  tip.  T  was  measured  from  the  area  under  the  load-dis¬ 
placement  curves  from  the  notched  tests. 

Comparisons  between  the  predictions  and  the  experiments 
are  presented  in  Fig.  1 1.  In  general,  the  predictions  lie  above  the 
measured  values.  Additionally,  the  magnitude  of  the  discrepan¬ 
cy  increases  from  roughly  10%  to  over  40%  as  the  PDA  con¬ 
centration  is  increased  from  0%  to  8%.  The  discrepancy  is 
believed  to  be  associated  mainly  with  the  shape  of  the  traction 
law,  as  described  below. 


The  stress-displacement  curves  in  Fig.  7,  as  well  as  previous 
theoretical  models  of  fiber  failure  and  pullout,  suggest  that  the 
bridging  traction  law  can  be  partitioned  into  two  parts5 
(Fig.  12).  The  first  is  characterized  by  rapid  hardening,  a  high 
peak  stress  (comparable  with  the  unnotched  strength),  and  a 
small  critical  displacement,  with  a  fracture  energy  denoted  rP  In 
this  domain,  the  inelastic  displacement  occurs  through  matrix 
cracking  and  associated  interfacial  debonding  and  sliding,  with 
minimal  fiber  failure.  The  second  part  is  characterized  by  expo¬ 
nential  softening,  a  low  strength  (that  because  of  pullout  of 
broken  fibers),  and  a  large  critical  displacement  (dictated  by  the 
pullout  length).  In  this  domain,  virtually  all  fibers  are  broken. 


SA  similar  two-part  traction  model  has  been  used  to  rationalize  notch  sensitivity  in  fiber- 
reinforced  Ti  matrix  composites.36 
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0  0.1  0.2  0.3  0.4  0.5 

Remote  Displacement  (mm) 

Fig.  7.  Stress-displacement  response  of  single-edge  notched  tension 
tests  (SENT)  specimens.  Symbols  represent  the  coordinates  correspond¬ 
ing  to  the  micrographs  in  Fig.  6. 


Fig.  8.  Macrophotographs  of  fractured  notched  tensile  specimens. 


0  5  10  15 

Notch  Length,  a0  (mm) 

Fig.  9.  Notched  strength  as  a  function  of  notch  length  and  precursor- 
derived  alumina  (PDA)  concentration. 


The  pullout  process  yields  additional  fracture  energy,  Tp,  yield¬ 
ing  a  total,  T  =  Tj+rp.  Despite  differences  in  their  shapes,  each 
of  the  two  bridging  laws  contributes  a  comparable  amount  to 
the  total  fracture  energy.  The  notch  length  range  in  which  each 


Shear  Strain  (%) 


Fig.  10.  Related  effects  of  precursor-derived  alumina  (PDA)  concen¬ 
tration  and  matrix  hardness  on  composite  shear  properties:  (a)  stress- 
strain  response  and  (b)  shear  strength.  Dashed  lines  in  (a)  denote 
limit  stress  values  at  larger  strains.  Also  shown  in  (b)  is  the  expected 
trend  for  a  ductile  matrix  composite  (dashed  line),  subject  to  the  Trcsca 
yield  criterion. 


a09  (P.  X)[o0F(a</W)(1-at/W)]2/[rE] 

Fig.  11.  Variation  in  notched  strength  with  normalized  notch  length. 
Solid  line  is  prediction  of  Eq.  (9).  Open  symbols  are  experimental  values, 
calculated  using  the  total  measured  fracture  energy;  solid  symbols  are 
experimental  values,  calculated  using  the  fracture  energy  without  the 
pullout  contribution. 
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Fig.  12.  Partitioning  of  the  bridging  traction  law.  Total  stress  is  the 
number-weighted  average  of  contributions  from  intact  and  broken  fib¬ 
ers.  Tj  and  Tp  represent  the  energies  associated  with  deformation  of  in¬ 
tact  fibers  and  pullout  of  broken  fibers,  respectively. 


has  a  significant  effect  on  notch  sensitivity  depends  on  the  char¬ 
acteristic  bridging  length  scale,  calculated  using  Eq.  (2),  with  the 
appropriate  value  of  fracture  energy.37  Specifically,  the  length  at 
pertinent  to  the  first  part  of  the  bridging  law  is  obtained  using  Tj 
in  Eq.  (2),  whereas  the  length  ap  where  pullout  plays  a  significant 
role  is  found  using  the  total  fracture  energy:  Fj+rp.  When  aQ< 
Oj,  notch  sensitivity  is  dictated  by  the  first  part  of  the  bridging 
law;  pullout  contributes  significantly  to  notched  strength  only  at 
larger  values  of  notch  length:  a0>ap.36 

This  energy  partitioning  scheme  was  implemented  in  the  fol¬ 
lowing  way.  The  pullout  contribution  was  obtained  from  the 
stress-displacement  curves  of  the  notched  specimens  at  large 
displacements  where  all  fibers  are  clearly  broken.  Assuming  that 
the  pullout  lengths  follow  an  exponential  distribution  and  that 
the  fibers  are  resisted  by  friction,  with  a  sliding  stress,  ts,  the 
predicted  response  for  a  2D  fiber  composite  in  a  0°/90°  orien¬ 
tation  is3839 


where  u  is  the  inelastic  remote  displacement,  R  is  the  fiber  ra¬ 
dius,  and  ho  is  the  mean  pullout  length.  Representative  curves 
from  each  of  the  three  materials,  shown  in  Fig.  13,  confirm  the 


Fig.  14.  Effects  of  precursor-derived  alumina  (PDA)  concentration  on 
both  interface  sliding  stress  and  mean  fiber  pullout,  as  determined  from 
the  pullout  response  (Fig.  1 3). 


Fig.  15.  Effects  of  matrix  hardness  on  the  total  fracture  energy.  T,  and 
the  fracture  energy  without  the  pullout  contribution,  f,. 


Inelastic  Displacement,  u  (mm) 


Fig.  13.  Fiber  pullout  response  following  bundle  fracture.  Solid  lines 
are  fits  of  Eq.  (13)  to  the  experimental  data. 


validity  of  Eq.  (13)  in  this  domain.  The  inferred  values  of  h0  and 
xs  are  plotted  in  Fig.  14.  Strictly,  the  contribution  from  pullout 
at  smaller  displacements  is  obtained  by  weighting  the  pullout 
stress  by  the  number  fraction  of  broken  fibers,  as  shown  sche¬ 
matically  by  the  dotted  curve  in  Fig.  12.  The  corresponding 
pullout  energy,  Tp,  is  the  integral  of  this  weighted  traction  law 
over  the  displacement  interval  O-oo .  Information  on  the  fraction 
of  broken  fibers  being  absent,  the  pullout  energy  is  estimated  by 
extrapolating  the  exponential  law  (assuming  all  fibers  are  bro¬ 
ken)  back  to  the  peak  stress  and  then  integrating  over  the  ap¬ 
propriate  limits.  Although  the  details  of  the  traction  law  are 
likely  in  error  in  the  transition  where  the  fiber  failure  rate  is  high 
(overestimating  the  bridging  stress  in  some  parts  and  underes¬ 
timating  it  in  others),  the  fracture  energy  is  likely  captured  ad¬ 
equately,  since  most  of  it  is  dissipated  at  larger  displacements. 
The  fracture  energy  pertinent  to  notch  sensitivity  in  the  short 
notch  domain,  T*,  is  then  taken  as  the  difference:  r-Tp.  The  re¬ 
sults  of  this  energy  partitioning  are  summarized  in  Fig.  1 5.  They 
demonstrate  that  Tp  represents  a  significant  fraction  of  T,  es¬ 
pecially  for  large  values  of  FPDA. 

To  justify  use  of  Tj  in  notched  strength  predictions,  the  values 
of  a,  and  ap  have  been  calculated  from  Eq.  (10)  and  the  fracture 
energies  in  Fig.  15,  and  compared  with  a0  in  Fig.  16.  In  the  intact 
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Fig.  16.  Comparisons  of  bridging  length  scales  a,  and  ap  with  the  notch 
lengths  used  in  the  present  study. 


fiber  domain,  a,  falls  in  the  range  of  10-40  mm,  decreasing  with 
increasing  FPDA.  This  range  is  comparable  with  (but  generally 
somewhat  greater  than)  the  notch  lengths  used  in  the  present 
experiments  (a0=  6-12  mm).  The  values  in  the  pullout  domain 
range  from  20  to  60  mm,  considerably  larger  than  a0.  The  in¬ 
ference  is  that  bridging  from  pullout  contributes  little  to  the 
notched  strength  over  the  current  notch  length  range,  but  would 
become  significant  for  notch  lengths  on  the  order  of  several 
centimeters. 

Revised  strength  predictions  utilizing  Tj  as  the  fracture  energy 
in  Eq.  (9)  are  shown  in  Fig.  11.  In  this  case,  the  data  for  ail 
materials  collapse  onto  essentially  a  single  band,  independent  of 
Epda-  Furthermore,  the  model  predictions  are  considerably 
closer  to  the  experimental  measurements  (within  about  I0%- 
20%).  These  results  provide  confidence  in  the  proposed  energy 
partitioning  method  for  ascertaining  the  contribution  pertinent 
to  short  notches  and  provides  preliminary  validation  of  the 
modified  strength  formula  in  Eq.  (9). 

Insights  into  the  effects  of  PDA  concentration  on  notch  sen¬ 
sitivity  are  obtained  by  considering  changes  in  the  parameter 
values  that  control  the  bridging  length  scale, j?ch  (Eq.  (10)).  Over 
the  PDA  range  probed  in  this  study,  a0  and  E  increase  by  about 
20%  and  30%,  respectively.  Being  fiber-dominated  properties, 
they  exhibit  the  weakest  sensitivity  to  the  state  of  the  matrix. 
Furthermore,  their  changes  have  essentially  offsetting  effects  on 
aCh.  In  contrast,  g(p,  X)  depends  on  a  matrix-sensitive  property, 
notably  the  shear  modulus.  Consequently,  it  decreases  by  a 
greater  amount:  a  factor  of  about  2.  This  change  produces  a 
proportionate  twofold  increase  in  ach.  Tj  is  by  far  the  most  ma¬ 
trix  sensitive,  because  of  its  strong  dependence  on  the  composite 
shear  strength,  as  illustrated  in  Fig.  2.  Indeed,  it  decreases  by 
almost  an  order  of  magnitude,  resulting  in  a  proportionate 
(nearly  10-fold)  reduction  in  ach.  The  latter  reduction  in  ach  is 
significantly  greater  than  the  corresponding  increase  because  of 
g(p,X).  The  net  result  of  these  dependencies  is  a  decrease  in  ach, 
by  a  factor  of  about  4,  and  hence  an  increase  in  notch  sensitivity, 
in  accordance  with  Eq.  (3). 


VII.  Summary 

The  strength  and  notch  sensitivity  of  a  family  of  porous-matrix 
composites  have  been  examined  experimentally.  The  unnotched 
tensile  and  shear  strengths  can  be  increased  by  strengthening  the 
matrix,  via  addition  of  a  PDA,  but  the  benefits  are  offset  by  an 
increase  in  notch  sensitivity. 

A  model,  based  originally  upon  finite  element  simulations, 
has  been  used  to  rationalize  the  degree  of  notch  sensitivity.  The 


simulations  account  for  the  interactions  between  notch  tip  ten¬ 
sile  and  shear  bands  in  semi-infinite,  notched  specimens  of 
elastically  isotropic  composites.  By  interpolating  between  the 
limiting  cases  of  very  short  and  very  long  notches,  a  simple  for¬ 
mula  for  notched  strength  is  obtained.  This  formula  has  been 
modified  to  account  for  the  effects  of  elastic  anisotropy  and  fi¬ 
nite  sample  size,  allowing  comparison  with  the  experimental 
measurements. 

Predictions  of  the  notched  strength  overestimate  the  experi¬ 
mental  measurements  when  the  total  (measured)  fracture  energy 
is  used.  This  discrepancy  arises  because  of  the  two  dissimilar 
parts  of  the  bridging  law.  To  remedy  the  problem,  a  method  for 
partitioning  the  bridging  law  and  the  associated  fracture  ener¬ 
gies  has  been  proposed  and  implemented.  For  the  range  of  notch 
lengths  used  in  the  present  study,  notch  sensitivity  is  dictated  by 
the  initial  rising  part  of  the  bridging  law,  where  the  fibers  are 
largely  intact.  Use  of  the  fracture  energy  from  this  part  of  the 
bridging  law  in  combination  with  the  notched  strength  formula 
yields  predictions  that  are  in  reasonable  agreement  with  the 
measurements.  The  choice  of  this  fracture  energy  can  be  justified 
on  the  basis  of  the  characteristic  bridging  length  scales  of  the 
two  parts  of  the  bridging  law  relative  to  the  notch  length. 

Effects  of  PDA  concentration  on  notch  sensitivity  can  be  ra¬ 
tionalized  in  terms  of  the  characteristic  bridging  length  scale, 
which,  in  turn,  is  dictated  by  the  interplay  between  a0,  E,  g( p,  X), 
and  Tj  in  Eq.  (3).  The  net  result  is  a  decreasing  ach  and  hence 
heightened  notch  sensitivity  with  increasing  PDA  concentration. 

Finally,  the  model  of  notch  sensitivity  presented  here  would 
be  applicable  to  other  ceramic  composites  with  weak  matrices, 
such  as  C/C  and  SiC/C.40  In  these  systems,  fracture  of  notched 
specimens  is  accompanied  by  formation  of  shear  bands  perpen¬ 
dicular  to  the  notch  plane,  in  a  manner  similar  to  that  of  the 
present  oxide  composites. 
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